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A  HIGH  POWER  PHASED  ARRAY  ANTENNA  FOR  SUB-SCALE 
AERIAL  TARGET  APPLICATIONS 

Junes  De  Vries  end  Kax  Dot 
Pacific  Missile  Test  Center 
Point  Mugu,  California  93042 

Fred  Lauriente 

H jrowave  Applications  Croup 
3030  Industrial  Parkway 
Santa  Maria,  California  93455 

ABSTRACT 

The  simulation  of  airborne  radars  place  special  requirements  on 
antenna  systems  used  in  sub-scale  targets.  Typical  airborne  radar 
antennas  are  steerable  in  azimuth  and  elevation  and  provide  gains 
on  the  order  of  3u  dB.  The  simulation  of  these  antennas  on  a  sub¬ 
scale  target  has  been  a  difficult  task.  Low  gain  non-steerable 
antennas  provide  good  coverage  but  suffer  from  marginal  ERP. 
Gimballed  antenna  systems  are  steer? Me  and  provide  nigh  gain  but 
not  without  the  attendant  problems  of  excessive  weight,  large  size 
and  high  cost.  In  an  effort  to  find  a  practical  solution,  a 
limited  scan  prototype  phased  array  antenna  has  been  developed. 

The  prototype  phased  array  effectively  utilizes  the  available 
aperture  to  provide  a  high  gain  antenna  that  is  steerable  over 
small  angles  in  azimuth  and  elevation. 

The  phased  array  antenna  described  herein  is  iaeal  for  this 
application  due  to  its  compact  size,  light  weight,  and  low  cost. 
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1.0  Introduction 


A  steerable  phased  array  antenna  system  was  designed  and  developed 
for  aerial  subscale  target  testing.  The  principal  purpose  of  this 
design  effort  was  to  increase  the  effective  radiated  power  (ERP)  by 
a  factor  of  d  dB  over  fixed  beam  antennas.  The  antenna  system 
fully  utilizes  the  limited  available  space  in  subscale  targets.  A 
microprocessor  control  subsystem  derives  steering  commands  from 
analog  navigational  and  altimeter  data,  and  provides  drive  signals 
to  steer  the  antenna.  The  operational  requirements  for  the 
antenna  system  of  ±20  degrees  in  azimuth  and  ±10  degrees  in  eleva¬ 
tion  were  satisfied  by  using  only  eight  elements.  This  minimal 
number  of  elements  reduces  system  complexity  while  providing  a  low 
cost  antenna  system. 

The  phased  array  system  is  comprised  of  the  following  components: 
a  two  bv  four  array  of  horn  elements,  eight  high  accuracy  rotary 
field  phase  shifters,  an  eight-wav  equal  line  length  waveguide 
power  divider  and  a  beam  steering  computer.  Each  of  those 
components  can  he  seen  in  the  photograph  of  Figure  1-1  and  is 
discussed  individually. 
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2.0  Horn  Radiator 


! 

I 

j  The  horn  radiator  is  an  E-plane  sectoral  horn  fed  by  a  six  element 

i 

i  »  slotted  waveguide  standing  wave  array.  An  exploded  view  of  the 

horn  is  shown  in  Figure  2-1.  The  slotted  waveguide  array  is  center 
fed  through  a  series  tee  junction  with  three  radiating  slots  on 

I 

j  either  side  of  the  junction.  With  an  impedance  matching  structure, 

I 

consisting  of  a  resonant  iris  and  dielectric  plug  in  the  input  arm, 

J  an  impedance  match  of  better  than  17  dB  was  achieved  over  the  14.0 

to  15.2  GHz  bandwidth  for  each  horn,  as  shown  in  Figure  2-2. 

i 

The  sectoral  horn  contains  dielectric  mode  launching  obstacles 

i 

attached  to  the  broad  walls  for  the  purpose  of  correcting  the  phase 
curvature  in  the  E-plane,  A  similar  set  on  the  side  walls  provides 

i 

I 

j  control  in  the  H-plane.  The  use  of  this  technique  has  been  described 

in  the  literature  (1,2).  The  radiation  patterns  of  the  horn  are 
shown  in  Figure  2-3  at  14.4  GHz.  The  measured  gain  was  20.7  dB 
.  corresponding  to  a,i  aperture  efficiency  of  approximately  80.2  per¬ 

cent. 

I 

i 

The  horn  radiator  was  subjected  to  high  power  testing  at  15  KW  peak 
and  15  watts  average  power  at  an  altitude  of  45,000  feet  without 
evidence  of  corona  or  voltage  breakdown. 
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3.0  Power  Divider  Network 


The  power  divider  network  consists  of  commercially  available  wave¬ 
guide  components.  Seven  folded  H-plane  WR-62  waveguide  magic  tees 
are  used  as  shown  in  Figure  3-1.  The  series  port  of  each  tee  is 
terminated  in  a  short  profile  slab  load,  capable  of  dissipating  one 
watt  of  RF  power.  The  power  divider  was  tested  for  loss,  equality 
of  power  split  and  relative  phase.  These  data  are  shown  in  Figures 
3-2  and  3-3. 

4.0  Rotary  Field  Phase  Shifters 

The  eight  phase  shifters  used  in  this  array  are  high  accuracy 
rotary  field  ferrite  phase  shifters.  These  phase  shifters  are  the 
magnetic  analog  of  the  Fox  rotary  vane  phase  shifter.  Typical  per¬ 
formance  data  is  shown  in  Figures  4-1  through  4-3. 

5.0  Beam  Steering  Controller 

The  Beam  Steering  Controller  inputs  are  in  the  form  of  three  analog 
voltages  corresponding  to  target  range,  bearing,  and  altitude  as 
measured  from  a  shipboard  TACAN. 

The  computer  receives  the  analog  signals  through  an  A/D  converter. 
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The  Z80  based  computer  with  CP/M  and  Fortran  software  generates  the 
appropriate  8-bit  sine  and  cosine  data.  The  eight  pairs  of  drivers 
receive  the  data  sequentially  and  latch  it  into  D/A  converters. 

The  power  stage  of  the  driver  aonsists  of  2  parts*  a  precision 
current  sink  and  polarity  switching  network.  The  current  sink 
employs  a  single  power  transistor  and  differential  current  sensing 
to  eliminate  crosstalk.  The  polarity  switching  network  which  re¬ 
ceives  polarity  information  from  the  8-bit  sine  or  cosine  word,  uses 
semiconductor  switches  to  effectively  switch  the  phase  shifter  coil 
orientation  in  the  circuit.  This  design  insures  good  bi-polar 
symmetry  and  eliminates  the  nead  for  a  negative  power  supply. 

The  beam  steering  controller  block  diagram  is  shown  in  Figure  5-1. 

6.0  Array  System  Tests 

System  evaluation  tests  included  radiation  patterns  and  gain 
measurements  throughout  the  angular  coverage  at  frequencies  of  14.0 
GHz  and  14.4  GHz  and  15.2  GHz.  A  composite  of  the  H-plane  scanned 
patterns  is  shown  in  Figure  6-1  and  a  composite  of  the  E-plane 
scanned  patterns  is  shown  in  Figure  6-2. 
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DESIGN  CONSIDERATIONS  FOR  THE  BEAMWAVEGUIDE 
RETROFIT  OF  A  GROUND  ANTENNA  STATION* 

by 

T.  Veruttipong,  J.  Withingtcn,  V.  Galindo-Israel , 
W.  Imbriale,  D.  Bathker 

Jot  Propulsion  Laboratory 
California  Institute  of  Technology 

j 


1.  (  INTRODUCTION 

> 

A  primary  requirement  of  the  NASA  Deep  Space  Network  (DSN) 
is  to  provide  for  optimal  reception  of  very  low  signal  levels. 

This  requirement  necessitates  optimizing  the  antenna  gain  to  the 
total  system  operating  noise  level  quotient.  Low  overall  system 
noise  levels  of  16  to  20  K  are  achieved  by  using  cryogenically 
cooled  preamplifiers  closely  coupled  with  an  appropriately 
balanced  antenna  gain/ spillover  design.  ^Additionally ,  high-power 
transmitters  (up  to  400  kW  CW)  are  required  for  spacecraft 
emergency  command  and  planetary  radar  experiments.  The  frequency 
bands  allocated  for  deep  space  telemetry  are  narrow  bands  near 
2.1  and  2.3  GHz  (S-band) ,  7.1  and  8.4  GHz  (X-band) ,  and  32  and 
34.5  GHz  (Ka-band).  In  addition,  planned  operations  for  the 
Search  for  Extraterrestrial  Intelligence  (SETI)  program  require 
continuous  low-noise  receive  coverage  over  the  1  to  10  GHz  band. 

To  summarize,  DSN  antennas  must  operate  efficiently  with  low 
receive  noise  and  high-power  uplink  over  the  1  to  35  GHz  band,  yL.. 

J 

Feeding  a  large  low-noise,  ground-based  antenna  via  a 
beamwaveguide  system  has  several  advantages  over  directly  placing 
the  feed  at  the  focal  point  of  a  dual-shaped  antenna.  For 
example,  significant  simplifications  are  possible  in  the  design 
of  high-power,  water-cooled  transmitters  and  low-noise  cryogenic 
amplifiers,  since  these  systems  do  not  have  to  rotate  as  in  a 
normally  fed  dual  reflector.  Furthermore,  these  systems  and 
other  components  can  be  placed  in  a  more  accessible  location, 
leading  to  improved  service  and  availability.  Also,  the  losses 
associated  with  rain  on  the  feedhorn  radome  are  eliminated 
because  the  feedhorn  can  be  sheltered  from  weather. 


Many  existing  beamwaveguide  systems  use  a  quasi-optical 
design,  based  on  Gaussian  wave  principles,  which  optimizes 
performance  over  an  intended  operating  frequency  range.  These 
designs  can  be  made  to  work  well  with  relatively  small  reflectors 


*  The  research  in  this  paper  was  carried  out  by  The  Jet  Propul¬ 
sion  Laboratory,  California  Institute  of  Technology,  under 
contract  with  the  National  Aeronautics  and  Space  Administration. 


(a  very  few  tens  of  wavelengths) ,  and  may  te  viewed  as  "band¬ 
pass/*  since  performance  suffers  as  the  wavelength  becomes  very 
short  as  well  as  very  long.  The  long  wavelength  end  is  naturally 
limited  by  the  approaching  small  D/X  of  the  individual  beam 
reflectors  used;  the  short  wavelength  end  does  not  produce  the 
proper  focusing  needed  to  image  the  feed  at  the  dual-reflector 
focus.  In  contrast,  a  purely  geometrical  optics  (G.O.)  design 
has  no  upper  frequency  limit,  but  performance  suffers  at  long 
wavelengths.  These  designs  may  be  viewed  as  "highpass." 
Considering  the  need  for  practically  sized  beam  reflectors  and 
the  high  DSN  frequency  and  performance  requirements,  the  G.O. 
design  is  favored  in  this  application. 

Retrofitting  an  antenna  that  was  originally  designed  without 
a  beamwaveguide  introduces  special  difficulties  because  it  is 
desirable  to  minimize  alteration  of  the  original  structure.  This 
may  preclude  accessing  the  center  region  of  the  reflector 
(typically  used  in  conventional  beamwaveguide  designs) ,  and  may 
require  bypassing  the  center  region.  A  discussion  of  the 
mechanical  tradeoffs  and  constraints  is  given  herein,  along  with 
a  performance  analysis  of  some  typical  designs.  In  the  retrofit 
design,  it  is  also  desirable  to  image  the  original  feed  without 
distortion  at  the  focal  point  of  the  dual-shaped  reflector.  This 
will  minimize  gain  loss,  reflector  design,  and  feed  changes. 

To  obtain  an  acceptable  image,  certain  design  criteria  are 
followed  as  closely  as  possible.  In  1973,  Mizusawa  and 
Kitsuregawa  [1]  introduced  certain  G.O.  criteria  which  guarantee 
a  perfect  image  from  a  reflector  pair  (cell) .  If  more  than  one 
cell  is  used  (where  each  cell  may  or  may  not  satisfy  Mizusawa* s 
criteria) ,  application  of  other  G.O.  symmetry  conditions  can  also 
guarantee  a  perfect  image.  The  problems  and  opportunities 
associated  with  applying  these  conditions  to  a  34-m  dual-shaped 
antenna  are  discussed. 

The  use  of  various  diffraction  analysis  techniques  in  the 
design  process  are  also  discussed.  Gaussian  (Goubau)  modes 
provide  important  insight  to  the  wave  propagation  character¬ 
istics,  but  Geometrical  Theory  of  Diffraction  (GTD) ,  FFT, 
Spherical  Wave  Expansion  (SWE) ,  and  Physical  Optics  (PO)  have 
proven  more  accurate  and  faster.  GTD  and  FFT  algorithms  are 
particularly  necessary  at  the  higher  frequencies.  PO  and  S\E 
have  been  necessary  at  the  lower  frequencies. 


2.  DESIGN  CONSIDERATIONS 
2 • 1  Hiahoass  Design  Feed  Imaging 

For  a  2.4-m  (8-foot)  reflector  beam  waveguide  operated  over 
1  to  35  GHz  with  near  perfect  imaging  at  X-band  and  Ka-band  and 
acceptable  performance  degradation  at  L-band  and  S-band,  a 
highpass  type  beamwaveguide  should  be  used.  This  type  of  design 
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is  based  upon  G.O.  and  Mizusawa's  criteria,  Mizusawa's  criteria 
can  be  briefly  stated  as  follows: 

For  a  circularly  symmetric  input  beam,  the  conditions  on  a 
conic  reflector  pair  necessary  to  produce  an  identical  output 
beam  are: 

(1)  The  four  loci  (two  of  which  may  be  coincident) 
associated  with  the  two  curved  reflectors  must  be 
arranged  on  a  straight  line. 

(2)  The  eccentricity  of  the  second  reflector  must  be  equal 
to  the  eccentricity  or  the  reciprocal  of  the 
eccentricity  of  the  first  reflector. 

Figures  la,  b,  and  c  show  some  of  the  orientations  of  the 
curved  reflector  pair  that  satisfy  Mizusawa's  criteria.  We  term 
a  curved  reflector  pair  as  one  cell. 

For  the  case  of  two  cells  where  at  least  one  cell  does  not 
satisfy  Mizusawa's  criteria,  a  perfect  image  may  still  be 
achieved  by  imposing  some  additional  conditions,  described  below. 

Let  ,  S2,  S3,  and  S4  be  curved  surfaces  of  two  cells  as 
shown  in  Figure  2.  Each  surface  can  be  an  ellipsoid,  hyper¬ 
boloid,  or  paraboloid.  Keeping  the  same  sequence  order,  the 
surfaces  are  divided  into  two  pairs  [first  pair  (S2,  S3)  ,*  second 
pair  (Slf  S4 ) ]  as  shown  in  Figure  2. 

For  a  circularly  symmetric  input  pattern,  an  identical 
output  pattern  (in  the  G.O.  limit)  can  be  obtained  if  Mizusawa's 
criteria  are  satisfied  for  both  pairs  in  the  following  manner: 
First  pair  (S2,  S3)  satisfies  Mizusawa's  criteria;  second  pair 
( ,  S4)  satisfies  Mizusawa's  criteria  after  eliminating  the 
first  pair. 

It  is  noted  that  the  first  pair  can  be  eliminated  because 
the  input  is  identical  to  the  output.  Also,  this  concept  can  be 
applied  to  cases  with  more  than  two  cells.  Examples  of  an 
extension  of  Mizusawa's  criteria  for  a  multiple-reflector 
beamwaveguide  are  given  in  Figures  3  and  4. 

Figure  5  shows  a  geometrical  optics  field  reflected  from 
each  reflector  of  a  beamwaveguide  system  (refer  to  Figure  3a) . 

It  is  clear  from  Figure  5  that  the  distorted  pattern  from  the 
first  cell  is  completely  compensated  for  by  the  second  cell  and 
yields  an  output  pattern  identical  to  the  input  pattern. 

2 . 2  Bandpass  Design  Feed  Imaging 

For  many  systems,  a  single-frequency  or  bandpass  design  can 
be  advantageously  employed.  The  design  considerations  can  best 
be  described  with  reference  to  Figure  6,  where  the  center 
frequency  is  given  as  fQ,  and  L2  is  the  spacing  between  two 


curved  surfaces.  A  bandpass  beamwaveguide  system  is  usually 
composed  of  two  non-confocal  (shallow)  ellipsoids  (eccentricity 
close  to  one).  Again  from  Figure  6,  Fjq  and  FA2  are  G.O.  foci  of 
ellipsoid  A,  while  FAp  is  the  phase  center  of  the  scattered  field 
from  surface  A  (evaluated  at  frequency  =  fQ)  in  the  neighborhood 
of  surface  B.  Similarly,  Fg1#  fB2,  anc*  FBP  are  f°r  ellipsoid  B. 
The  distances  from  FA2  and  fap  from  surface  A  are  very  large 
compared  to  and  L2 .  The  locations  of  FAp  and  Fgp  depend  on 
frequency  as  well  as  surface  curvature,  ,  and  L2.  For  example, 
with  a  2.4-m  reflector  with  eccentricity  =  0.97  at  f  =  2.3  GHz 
and  L2  =  8  m  (26  feet) ,  FAp  is  about  120  m  (400  feet)  to  the  left 
of  ellipsoid  A,  as  shown  in  Figure  6a.  In  the  G.O.  limit,  FAP 
and  Fa2  are  at  the  same  location,  to  the  right  of  ellipsoid  A. 

It  is  desirable  to  have  two  identical  surfaces  for  low 
cross-polarization  and  a  symmetrical  system.  Trial  and  error  is 
needed  in  order  to  determine  surface  parameters  for  the  desired 
operating  frequency  and  bandwidth  within  specified  losses. 

Figure  7  shows  the  input  and  output  patterns  from  a  bandpass 
beamwaveguide  system  where  FAp  and  FBp  are  chosen  to  be  at  the 
same  locations  as  F^2  and  FA2 ,  respectively  (the  choice  may  not 
be  the  optimum  condition) .  The  two  identical  ellipsoids  are 
designed  at  fQ  =  2,3  GHz.  The  results  show  good  agreement 
between  the  input  feed  and  the  imaged  feed.  Bandpass  beamwave¬ 
guide  systems  appear  useful  when  a  limited  band  coverage  is 
required,  using  modestly  sized  (D  =  20  to  30A)  reflectors. 
However,  these  systems  do  not  perform  well  as  the  wavelength 
approaches  either  zero  or  infinity.  In  contrast,  a  highpass 
(G.O.)  design  focuses  perfectly  at  zero  wavelength  and  focuses 
very  well  down  to  D  ~40A.  The  performance  then  decreases 
monotonically  as  D  becomes  smaller  in  wavelengths. 


3.  APPLICATION  CONSIDERATIONS  FOR  THE  DSN 

The  DSN  presently  operates  three  34-m  high  efficiency  (H.E.) 
dual-shaped  reflector  antennas  with  a  dual-band  (2 . 3/8 . 4-GHz) 
feed  having  a  far-field  gain  of  +22.4  dBi  that  is  conventionally 
located  at  the  Cassegrain  focal  point.  The  structures  were 
designed  prior  to  beamwaveguide  requirements,  and  feature  a 
continuous  elevation  axle  and  a  carefully  designed  elevation 
wheel  substructure.  The  elevation  wheel  substructure,  shown  in 
Figure  8,  plays  a  key  role  in  preserving  main  reflector  contour 
integrity  as  the  antenna  rotates  in  elevation.  To  maintain 
contour  integrity  at  8  GHz  and  above  is  of  prime  concern  for  RF 
efficiency  performance  as  well  as  retrofit  costs.  Figure  9  shows 
a  centerline  beamwaveguide  approach  which,  although  it  is  a 
compact  and  straightforward  RF  design,  severely  impacts  retrofit 
costs  and  the  contour  integrity  of  the  main  reflector  (hence  RF 
efficiency  at  8  GHz  and  above) .  Figure  10  shows  an  unconven¬ 
tional  approach  and  represents  attempts  to  reduce  structural 
impacts.  Figure  10  is  a  six-reflector  beamwaveguide  (eight- 
reflector  antenna)  based  on  our  extension  of  Mizusawa's  criteria. 


390- 


Although  several  detailed  options  are  possible,  most  options  use 
two  cells  (four  curved  reflectors)  with  two  flat  reflectors. 

Some  of  the  options  make  use  of  the  flexibility  afforded  by 
allowing  each  cell  to  be  distorting  (of  itself) ,  but  then 
compensated  for  by  the  second  cell  as  described  in  Figures  3,  4, 
and  5. 

Our  goal  is  therefore  to  perfectly  image  a  feed  located 
perhaps  15  to  25  m  (50  to  80  feet)  below  the  main  reflector  to 
the  original  Cassegrain  focus.  This  goal  applies  over  the  1  to 
35-GHz  frequency  range,  using  a  beamwaveguide  housing  limited  to 
about  2.4  m  (8  feet)  in  diameter.  The  image  should  be  a  1:1 
beamwidth  transformation  of  the  original  +22.4  dBi  feed, 
permitting  reuse  of  that  feed  and  no  changes  to  the  subreflector 
or  main  reflector  contour.  Lastly,  the  goal  includes  minimal 
structural  impacts,  particularly  to  the  integrity  of  the  main 
reflector  contour.  The  generalized  solution  to  these  goals  is 
reflected  in  the  approach  shown  in  Figure  10,  termed  the  bypass 
beamwaveguide . 


4.  ANALYTICAL  TECHNIQUES  FOR  DESIGN  AND  ANALYSIS 

The  software  requirements  for  the  study  and  design  of 
beamwaveguides  are  extensive.  These  include  the  capability  for 
G.O.  synthesis,  Gaussian  wave  analysis,  and  high  and  low 
frequency  diffraction  analysis.  These  requirements  are  discussed 
below. 

4 . 1  Geometrical  Optics  Synthesis  Capability 

This  capability  includes  software  which  synthesizes  as  well 
as  analyzes  reflectors  satisfying  the  Mizusawa-Kitsuregawa 
conditions  [l]  for  minimum  cross-polarization  and  best  imaging. 

In  the  high-frequency  domain  (8  to  35  GHz  for  the  designs 
considered  herein) ,  the  focused  system  shown  in  Figure  ll  is 
desired.  Of  course,  two  paraboloids  or  mixtures  of  various 
conic-section  reflectors  can  be  synthesized.  Optimization  at 
lower  frequency  bands  generally  is  accomplished  by  appropriate 
defocusing. 

4 . 2  Gaussian  Wave  Analysis  Capability 

The  required  defocusing  for  the  lower  bands  can  be  deter¬ 
mined  by  using  various  beam  imaging  techniques  [2]  which  are 
based  on  Gaussian  beam  analysis  methods  first  developed  by  Goubau 
and  Schwering  [3],  While  Gaussian  mode  analysis  is  useful  at 
high  as  well  as  low  frequencies  for  "conceptual"  designing, 
conventional  diffraction  analysis  methods  are  found  to  be  more 
suitable . 

One  consideration  is  that  Gaussian  mode  analysis  does  not 
supply  spillover  losses  directly  with  as  great  an  accuracy  as 
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conventional  diffraction  analysis  methods.  The  Gaussian  modes 
supply  the  discrete  spectrum,  whereas  the  spillover  losses  are 
directly  related  to  the  continuous  spectrum  [4].  Spillover  for 
Gaussian  modes  is  computed  by  a  method  more  suitable  in  terms  of 
computational  efficiency  for  a  great  number  of  reflection 
(refraction)  elements  [3].  The  antenna  systems  considered  here 
rarely  contain  more  than  four  curved  reflectors.  Further, 
spillover  losses  must  be  reliably  known  to  much  less  than  a  tenth 
of  a  decibel  for  high  performance  systems. 

4 . 3  Low  Frequency  Diffraction  Analysis  Capability . 

Because  of  the  large  bandwidth  of  operation  (1  to  35  GHz) , 
there  is  no  single  diffraction  analysis  method  which  will  be  both 
accurate  and  efficient  over  the  entire  band. 

Efficiency,  in  the  sense  of  speed  of  computation,  is  criti¬ 
cal,  since  in  a  constrained  design  many  different  configurations 
may  be  analyzed  before  an  optimum  beam  waveguide  configuration  is 
selected. 

In  the  region  of  2.3  GHz  (for  the  34-m  antenna  considered 
herein) ,  the  reflector  diameters  are  generally  about  20A.  Since 
a  very  low  edge  taper  illumination,  at  least  -20  dB,  is  used  to 
reduce  spillover  loss,  the  "effective'1  reflector  diameters  are 
very  small  in  this  frequency  range. 

A  comparison  between  three  diffraction  algorithms:  PO,  GTD, 
'■nd  Jacobi-Bessel  (JB)  leads  to  the  conclusion  that: 

(1)  GTD  is  not  sufficiently  accurate  at  the  low 
frequencies . 

(2)  JB  is  very  slowly  convergent  in  many  cases  and  gives 
only  the  far-field  in  any  case.  We  must  determine 
near-field  patterns. 

(3)  PO  is  both  accurate  and  sufficiently  fast  below  3  GHz. 

'-V.  above  results  are  illustrated  in  Figures  12  and  13  when  an 
ac  .al  +22.4  dBi  corrugated  feedhorn  is  used.  It  can  be  seen 
fr>  Figure  13  that  PO  and  JB  agree  perfectly  to~+20°. 

There  are  two  general  PO  computer  algorithms  useful  at  the 
low  f  aquencies.  One  is  a  straightforward  PO  algorithm  which 
subdi\ ides  the  reflectors  into  small  triangular  facets.  This  is 
essentially  a  trapazoidal  integration  of  the  near-field  radiation 
integral  and  is  a  very  flexible  algorithm. 

A  second  algorithm  is  based  on  a  Spherical  Wave  Expansion 
(SWE)  and  is  also  a  PO  technique.  It  has  two  useful  character¬ 
istics:  (1)  when  a  high  degree  of  circular  symmetry  exists  in 

the  scattered  fields,  then  the  two-dimensional  radiation  integral 
is  reduced  to  a  small  number  of  one-dimensional  integrals  with  a 


resultant  marked  decrease  of  computational  speed;  and  (2)  an  (r) 
interpolation  of  the  scattered  field  (at  different  radial 
distances  from  the  coordinate  origin)  is  done  very  accurately. 

The  PO  (direct-trapezoid)  and  SWE  algorithms  are  useful  for 
cross-checking  results.  Near-field  and  far-field  computations 
and  comparisons  in  both  amplitude  and  phase  are  shown  in  Figures 
14,  15,  and  16.  Figure  14  contains  near-field  PO  and  SWE  results 
for  scattering  from  one  paraboloid  reflector.  Figure  15  contains 
the  results  for  far-field  scattering  from  two  reflectors.  Figure 
16  contains  the  near-field  scattering  for  two  reflectors  and  a 
comparison  of  the  "imaging"  with  the  feed  pattern.  At  higher 
frequencies  (>8  GHz) ,  the  feed  pattern  will  be  virtually 
perfectly  imaged  over  an  angle  of  +15°. 

4 • 4  High  Frequency  Diffraction  Analysis  Capability 

For  reflector  diameters  of  70  or  more  wavelengths  (>8  GHz) , 
including  a  -20  dB  edge  taper,  PO  analysis  methods  become  too 
expensive  and  time  consuming.  (The  SWE  may  still  be  useful  if  a 
high  degree  of  rotational  symmetry  exists.)  An  alternative 
approach  is  to  use  GTD  analysis.  The  GTD  computation  speed  does 
not  increase  with  increasing  reflector  diameters,  but  the 
accuracy  of  the  analysis  does  increase. 

In  order  to  test  the  accuracy  of  GTD  at  8  GHz,  comparisons 
were  made  between  GTD  and  PO.  Results  for  diffraction  of  a 
single  ellipsoid  are  shown  in  Figure  17.  The  results  for  the 
phase  of  the  scattered  field  were  equally  as  good.  GTD  was 
determined  to  be  accurate  at  8  GHz  and  higher. 

Analysis  of  two  or  more  reflectors  by  GTD  involves  some 
manipulation  of  the  fields  scattered  between  any  two  reflectors. 
The  fields  scattered  from  one  reflector  must  be  placed  in  a 
format  suitable  for  GTD  scattering  from  the  next  reflector.  This 
is  accomplished  by  computing  the  vector-scattered  field  in  the 
vicinity  of  the  next  reflector  and  then  interpolating  as  follows: 

(1)  Use  of  an  FFT  for  0-variable  interpolation. 

(2)  Use  of  a  second-order  Lagrangian  local  interpolation 
for  0  interpolation  (for  a  z-axis  along  the  axis  of  the 
reflector,  0  and  9  are  spherical  coordinates) . 

(3)  For  the  r  (of  [r,  9,  0]  spherical  coordinates) 
interpolation,  an  approximation  consistent  with  the  GTD 
approximation  was  to  assume  a  1/r  variation  in  ampli¬ 
tude  and  a  kr  variation  in  phase.  This  approximate 
interpolation  should  be  checked  against  exact  computa¬ 
tions  of  the  near  fields. 

By  the  method  described  above,  multiple  reflector  computa¬ 
tions,  even  with  a  large  number  of  reflectors,  can  be  calculated 
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with  both  great  speed  and  accuracy  at  frequencies  above  ~8  GHz 
for  reflectors  of  >70X  in  diameter. 

A  typical  result  is  shown  in  Figure  18  for  a  pair  of 
ellipsoids  which  satisfy  the  Mizusawa-Kitsuregawa  criteria.  The 
object  is  to  perfectly  image  the  input  feed  over  about  +20°. 

This  is  accomplished  with  great  accuracy  and  virtually  no  loss  at 
X-band  frequencies  and  higher. 


5.  CONCLUSIONS 

A  generalized  solution  has  been  achieved  for  retrofitting  a 
large  dual-shaped  reflector  antenna  for  beamwaveguide.  The 
design  is  termed  bypass  beamwaveguide.  Several  detailed  options 
within  the  bypass  category  remain  to  be  studied,  and  work 
continues. 

With  the  analysis  capability  available,  we  are  gaining  some 
valuable  views  of  the  RF  performance  behavior  of  some  of  the  many 
options.  It  appears  fairly  clear  for  the  1  to  35-GHz  requirement 
that  highpass  (pure  G.O.)  designs  are  necessary  in  contrast  to  a 
bandpass  (non-confocal  ellipsoids)  approach.  It  appears  that 
deep  confocal  ellipsoids  satisfying  the  Mizusawa-Kitsuregawa 
criteria  operate  (focus)  well  with  reflector  diameters  of  about 
7 OX  and  larger  but  may  not  be  tolerable  at  longer  wavelengths. 

As  a  part  of  this  activity,  an  important  extension  of  the 
Mizusawa-Kitsuregawa  criteria  has  been  revealed.  The  principle 
revealed  shows  how  a  two-reflector  (cell) ,  although  in  itself 
distorting,  may  be  combined  with  a  second  cell  which  compensates 
for  the  first  and  delivers  an  output  beam  which  is  a  good  image 
of  the  input  beam. 
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SECOND  PAIR 


Figure  2.  Multiple  Curved  Reflector  Beainwaveguide  System; 
Si  =  Ellipsoid,  Hyperboloid,  or  Paraboloid 
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I  (s2,  S3)  SATISFIES  MIZUSAWA 
V  AND  IS  ELIMINATED 


Figure  3.  Demonstrating  an  Extension  of  Mizusawa's  Criteria 

for  a  Multiple  Reflector  Beamwaveguide  (Paraboloids) 
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AND  IS  ELIMINATED 


(Su  $4)  SATISFIES  MIZUSAWA 


Figure  4.  Demonstrating  an  Extension  of  Mizusawa's  Criteria 
for  a  Multiple  Reflector  Beamwaveguide 
(Ellipsoids/Hyperboloids) 
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Figure  5.  Geometrical  Optics  Field  Reflected  from  Each  Surface 
of  a  Beamwaveguide  System  (Shown  in  Figure  3a) 
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Figure  9.  34-m  H.E.  Centerline  Beamwaveguide 
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Figure  10.  34—  xn  H.E.  Bypass  Beainwaveguide 
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Analysis  Geometry  for  Ellipsoid  Pair  (EL^,  EL2) 
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Figure  12.  Near-Field,  One  2.4-m  Ellipse  Diffraction  Patterns 
(RCP,  2.3-GHz,  Offset  Plane) 
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Figure  16.  PO  and  SWE  Near-Field  Scattering  from  Two  2.4-m 
Paraboloids  at  2 . 3  GHz  (comparison  with  Feed 
Pattern) 
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Figure  17.  Near-Field,  One  2.4-m  Ellipse  Diffraction  Patterns 
(RCP,  8.4  GHz ,  Offset  and  1  Plane  Cuts)  * 
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Figure  18.  GTD  Scattering  for  Two  2.4-m  Ellipses,  RCP  at  8.4  GHz 
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AN  INEXPENSIVE  RELATIVELY  BROADBAND  MULTI-PURPOSE  ANTENNA 

Samuel  C.  Kuo  and  Warren  M.  Shelton 
GTE  Government  Systems,  Western  Division 

'j/ 

The  one-wavelength  loop  antenna  has  been  widely  used  In  the  HF 
and  low  VHF  range  by  amateur  radio  operators  In  the  form  of  the 
quad  antenna  ever  since  It  was  Introduced  In  the  late  1940 's  by 
Moore.  ^Three  types  of  the  one-wavelength  loop  antennas  are  shown 
In  Figure  1-a  through  1-c.  In  (a)  and  (b)  the  sides  of  the  square 
loops  are  equal  to  1/4  wavelength  and  the  sides  of  the  triangular 
loop  In  (c)  are  equal  to  1/3  wavelength.  The  relative  direction 
of  current  flow  is  also  shown  In  the  drawing  to  demonstrate  the 
polarization  of  these  antennas.  When  the  input  terminals  are  on 
the  bottom  and  are  excited  with  equal  amplitude  and  out-of-phase 
currents,  these  antenna  provide  horizontally  polarized  radiation. 
The  original  Quad  antennas  introduced  by  Moore,  shown  In 
Figure  1-d  and  1-e,  consisted  of  two  one-wavelength  loops;  one  as 
the  driven  element,  and  the  other  as  the  reflector.  Subsequent 
development  added  more  loops  as  directors  and  the  antenna  with 
directors  added  was  called  a  Quagi  antenna  which  is  claimed  to 
have  a  gain  of  approximately  2  dB  over  a  Yagi  of  the  same  array 
length. 

Various  types  of  the  Quad  and  Quagi  antennas  have  been  devel¬ 
oped  through  the  years,  and  they  have  always  been  used  as  rela- 
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Figure  2-a  has  better  characteristics  In  terms  of  VSWR,  gain,  and 
bandwidth  than  the  two  other  types.  It  was  selected  for  bandwidth 
comparison  with  a  single  loop. 

A  Bl-loop  antenna,  as  shown  In  Figure  2-a,  was  fabricated  with 
each  side  of  the  loop  antenna  approximately  4  Inches  long  and  was 
mounted  3  Inches  above  a  14  by  18  Inches  ground  plane.  Using  a 
VSWR  of  3.0:1  as  the  criteria,  this  antenna  has  a  bandwidth  of 
approximately  2.4:1,  while  a  single  loop  of  the  same  dimension  and 
configuration  as  one-half  of  the  Bl-loop  has  a  bardwldth  of 
approximately  2.0:1.  The  swept  gains  of  the  two  antennas  from  .55 
to  1.2  GHz,  shown  In  Figure  3,  shows  the  gain  of  the  Bi-loop 
antenna  to  be  approximately  2  dB  higher  than  that  of  a  single  loop 
antenna. 


The  bandwidth  of  a  dipole  or  monopole  antenna  depends  a  great 
deal  on  the  diameter  of  the  radiator.  The  same  Is  also  true  for 
the  loop  antenna.  The  Bi-loop  antenna,  snown  in  Figure  4,  was 
made  of  relatively  large  diameter  tubing.  The  bandwidth  of  this 
particular  antenna  Is,  however,  limited  by  the  balun.  Subsequent 
antennas  built  were  fed  with  coaxial  cables  as  shown  in  Figure  5. 
Since  the  characteristic  impedance  of  this  antenna  is  approxi¬ 
mately  75  ohms,  a  50/75-ohm  transformer  was  included  in  the  feed 
line  for  impedance  transformation.  When  a  Bi-loop  antenna  with 
4-1nch  sides  was  built  with  0.650  diameter  tubing,  and  fed  with 
75-ohm  coaxial  cables  and  a  50/75-ohm  impedance  transformer,  the 
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tlvely  narrow  band  antennas  by  the  amateurs.  For  Instance,  a 
popular  version  of  the  antenna,  the  "Three-band  Quad,"  which  has 
three  sets  of  different  sized  ful 1 -wavelength  loops,  is  used  to 
cover  14.4,  21.1,  and  28.1  MHz. 

In  the  high  VHF,  UHF  and  microwave  frequency  range  where  broad 
frequency  coverage  is  often  desired,  the  original  Quad  antenna 
which  has  limited  broadwldth,  can  not  meet  the  frequency  bandwidth 
requirements.  In  Moore's  design,  one  of  the  two  loop  antennas  is 
used  as  a  reflector,  which  Is  very  sensitive  to  the  operating  fre¬ 
quency  and  consequently  limits  the  bandwidth  of  the  antenna.  The 
bandwidth  of  the  Quad  antenna  can  be  increased  by  replacing  the 
reflector  loop  element  with  a  flat  ground  plane.  The  performance 
of  a  one-wavelength  loop  In  front  of  a  ground  plane  Is  similar  to 
that  of  a  dipole  over  ground.  When  the  spacing  Is  small  In  terms 
of  wavelength,  the  radiation  efficiency  Is  very  low  and  the  VSWR 
is  high.  When  the  spacing  approaches  one-half  wavelength,  a  null 
will  appear  on  the  antenna  axis. 

Further  Increase  of  the  operating  bandwidth  of  this  antenna 
can  be  achieved  by  arraying  two  one-wavelength  loop  antennas.  The 
uniqueness  of  the  one-wavelength  loop  antenna  is  that  the  two 
loops  can  be  arrayed  without  a  RF  combiner,  and  the  second  loop 
antenna  can  simply  be  paralleled  to  the  first  loop  antenna  at  the 
feed  point.  The  new  structures  are  called  Bi-loop  antennas,  three 
forms  of  which  are  shown  In  Figure  2.  Since  the  type  shown  in 
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bandwidth  (3.0:1  VSWR)  approached  3.0:1.  Unfortunately,  the  cost 
of  this  antenna  is  very  high  due  to  the  tubular  construction  and 
the  cost  of  the  transformer. 

An  Inexpensive  version  of  the  Bi-loop  antenna  was  developed 
which  costs  less  than  10X  of  the  cost  of  the  tubular  configura¬ 
tion.  This  antenna  was  fabricated  using  printed-circuit  board 
techniques  where  the  Bi-loop  antenna  element  was  etched  on  one 
side  of  the  PC  board  and  the  Impedance  transformer  etched  on  the 
other.  The  antenna  element  trace  serves  as  the  ground  plane  for 
the  microstrip  transformer  which  provides  the  Impedance  transfor¬ 
mation  from  50  ohms  at  the  Input  connector  to  75  ohms  at  the  Input 
to  the  Bl-loop  antenna.  A  sketch  of  this  Bi-loop  antenna  with  Its 
integral  Impedance  transformer  is  shown  in  Figure  6.  The  details 
of  the  connector  and  the  feed  point  configurations  are  also 
shown.  The  material  used  to  fabricate  this  Bl-loop  antenna  and 
its  ground  plane  consist  of  a  piece  of  G- 10  epoxy  fiberglass  PC 
board,  a  panel  type  connector  and  a  piece  of  aluminum  sheet  for 
ground  plane.  The  total  cost  is  less  than  $30.00  including  etch¬ 
ing  the  antenna.  A  slight  decrease  In  bandwidth  and  gain  are  the 
trade-off  for  the  low  cost  in  comparison  to  the  tubular 
conf iguratlon. 

One  such  antenna  was  fabricated  using  l/8-1nch  G-10  epoxy 
fiberglass  PC  board  with  the  following  parameters: 
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Length  of  each  side  of  the  loop  4.0  Inches 

Width  of  the  loop  conductor  0.75  Inches 

Spacing  from  ground  plane  3.0  Inches 

Size  of  ground  plane  11  by  14  Inches 

The  performance  characteristics  of  this  antenna  are  shown 
below,  and  typical  E-  and  H-plane  patterns  are  shown  In  Figure  7. 

VSWR:  600  -  1300  MHz  2.0  Max. 

Gain:  600  -  1200  MHz  7.5  dBI  Min.  8  dBi  Norn. 

3  -  dB  Beamwldth,  E-plane  65°  Nom. 

H-plane  65°  Nom. 

Since  these  antennas  are  fabricated  using  PC  board  etching 
techniques,  they  are  nearly  Identical  and  are  good  candidates  for 
linear  array  elements.  Figure  8  shows  a  10-element,  linear  array 
using  the  Bi-lcop.  Over  a  1.6:1  bandwidth,  the  array  can  be 
steered  to  +/-32°  with  the  VSWR  of  the  antenna  elements  remaining 
under  2.0:1.  Azimuth  radiation  patterns  of  this  array  measured  on 
axis  and  with  32°  scan,  are  shown  in  Flgu#  9  for  low,  medium,  and 
high  frequencies. 

. 
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Figure  1.  One  Wavelength  Loop  and  Quad  Antennas. 
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Figure  2.  Three  Types  of  Bi-Loop  Antennas. 


Figure  3.  Gain  Comparison  of  A.  Bi-Loop  and  A  single-Loop 
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Figure  4.  Bi-Loop  Antenna 
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Figure  6.  A  Bi-Loop  Antenna  on  PC  Board  with  Its 
Integral  Impedance  Transformer 
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Figure  8.  A  10-Element  Linear  Array  of  Bi-Loop  Antennas 
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ABSTRACT 


In  feeding  a  log-periodic  monopole  or  slot  array,  excess  line  must  be 
used  between  radiating  elements  to  achieve  proper  phasing  for  backfire 
radiation.  Discontinuities  presented  by  the  radiating  elements,  which  are 
approximately  one-half  wavelength  apart  on  a  series  feedline,  cause 
additive  reflections  on  the  line.  There  is,  therefore,  a  region  on  the  log- 
periodic  array  which  behaves  similarly  to  a  periodically  loaded  transmission 
line  operating  in  a  stopband.  Reflections  occuring  in  this  region  inhibit 
propagation  of  energy  down  the  feedline.  To  obtain  a  frequency-stable 
impedance,  these  ''structural  stopband*  regions  should  be  reduced. 

One  technique  that  can  be  used  to  minimize  the  stopband  region  on  a 
series  feedline  is  to  reduce  the  reactive  nature  of  the  non-radiating 
elments  using  the  principle  of  complementarity.  A  two-port  monopole-slot 
antenna  has  been  previously  shown  to  have  a  relatively  constant  image 
impedance  below  its  first  resonance.  -4 — 

In  the  work  presented  here,  a  log-periodically  scaled  array  of 
monopole-slot  elements  was  constructed  and  tested.  The  results 
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demonstrate  the  elimination  of  the  structural  stopband  regions  without  the 
use  of  modulation  of  the  feedline  impedance.  Other  performance! 

r 

characteristics  such  as  impedance  bandwidth,  beamwidth,  and  fraWto- 
back  ratio  of  the  array  are  also  presented. 

t 

1.  INTRODUCTION 

In  the  past  few  years  there  has  been  a  great  deal  of  interest  in 

broadband,  low-profile,  surface-mounted  antennas.  Although  microstrip 

patch  and  cavity-backed  slot  antennas  satisfy  the  requirements  of  low 

profile  and  surface  mountability,  neither  structure  can  be  considered 

broadband.  However,  log-periodic  arrays  of  cavity-backed  slots  are  able  to 

satisfy  the  broadband  requirement.  Theoretically,  the  operating  bandwidth 

of  such  an  array  is  limited  only  by  the  number  of  elements  used  and  the 

size  allowed  for  construction  of  the  array.  In  practice,  however,  there  is  a 

fundamental  problem  in  the  feeding  of  such  an  array. 

In  feeding  a  log-periodic  monopole  or  slot  array,  excess  line  must  be 

used  between  the  radiating  elements  to  achieve  proper  phasing  of  the 

elements  for  backfire  radiation.  Ingerson  and  Mayes  ^  found  that  the 

discontinuities  presented  by  slots,  which  were  approximately  one-half 

wavelength  apart  along  a  series  feedline,  cause  additive  reflections  on  the 
* 

line.^1  "Here  is,  therefore,  a  region  on  the  log-periodic  array  that  behaves 
similarly  to  a  periodically  loaded  transmission  line  operating  in  a  stop 
band.  Reflections  that  occur  in  this  stop  region  inhibit  the  propagation  of 
energy  down  the  feedline.  To  achieve  an  input  impedance  that  changes 

,  little  with  frequency,  "structural  stopband"  regions  between  the  source  and 

» 
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the  active  elements  should  be  minimized.  Ingerson  and  Mayes 
demonstrated  that  modulating  the  impedance  of  the  feedline  could  greatly 
reduce  this  stopband  phenomenon.  The  input  impedance  of  the  antenna 
constructed  by  Ingerson  was  approximately  150  ohms.  This  impedance  is 
too  high  tor  many  applications.  With  this  fact  in  mind,  Ostertag  applied 
the  impedance-modulation  technique  to  a  shallower  cavity  with  a 
microstrip  feedline  and  was  able  to  reduce  the  input  impedance  to  50 
ohms^. 

Another  technique  that  can  be  user)  to  minimize  the  stop  band  for  a 

series  feedline  is  to  reduce  the  reactive  nature  of  the  nonradiating 

elements  using  the  principle  of  complementarity.  A  planar  dipole  and  the 

cor  responding  complementary  slot  are  duals  of  each  other.  Using  Babinet's 

principle  it  can  be  shown  that: 

? 

M  V~ 

*7  T  _  _  _ _ [ _ _ 

^s^d”  4 «  "  4 

(1) 

where  Z$  is  the  slot  impedance  and  is  the  impedance  of  the  oomplentary 
dipole.  Thus  it  is  seen  that  the  input  impedance  of  a  slot  and  a 
complementary  dipole  fed  together  is  independent  of  frequency.  This 

3  >4  3 

principle  was  used  b>  Schroeder  ,  Dunlavy  ,  and  Itoh  and  Chen  .  However, 
all  of  these  structures  required  some  sort  of  external  feed  network  for 
proper  operation. 

In  1970  Mayes  and  Weisenmeyer^  demonstrated  that  a  two-port 
monopole-slot  antenna  can  be  made  to  have  a  relatively  constant  image 
impedance  for  frequencies  below  its  first  resonance,  A  simple  microstrip 
feediine  with  a  cylindrical  monopole  attached  directly  above  the  slot 
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feedpoint  was  used  with  this  antenna.  Stripline  versions  of  the  monopole- 

7  8  9 

slot  were  later  studied  by  Cwik  ,  Halpern  ,  and  Paschen  .  Paschen  also 
investigated  the  properties  of  the  monopole-slot  with  regard  to  using  the 
antenna  as  an  array  element.  He  found  that  in  addition  to  having  a  uniform 
image  impedance  (thus  minimizing  the  effect  of  the  structural  stopbands), 
the  element  also  has  other  characteristics  which  make  it  a  desirable 
element  for  a  log-periodic  array.  He  found  that  the  far-field 
characteristics  of  the  monopole  combined  with  the  slot  result  in  an  element 
with  a  fairly  high  front-to-back  ratio.  This  allows  for  easy  construction  of 
an  array  with  a  high  front-to-back  ratio.  He  also  demonstrated  that  this 
characteristic  of  the  monopole-slot  can  reduce  the  mutual  coupling 
between  the  cascaded  elements  of  a  series-fed  array. 

In  the  project  reported  herein,  a  log-periodical ly  scaled  array  of 
monopole-slot  antenr.a  elements  was  constructed  and  tested.  One  objective 
was  to  demonstrate  the  elimination  of  the  structural  stopbands  without  the 
use  of  modulation  of  the  feedline  impedance.  Other  objectives  were  to 
evaluate  the  performance  characteristics  of  impedance  bandwidth,  pattern 
bandwidth,  beamwidth,  and  front-to-back  ratio. 

2.  MEASUREMENTS  FOR  A  SINGLE  MONOPOLE-SLOT  ELEMENT 

Before  construction  of  the  complete  array  could  begin,  the 
impedance  characteristics  of  a  single  representative  monopole-slot  were 
needed.  The  impedance  characteristics  of  a  previously  constructed  cavity- 
backed  slot  were  determined  using  a  vector  network  analyzer  (Figure  1). 
Two  different  techniques  were  initially  used  to  achieve  an  impedance 
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characteristic  for  the  monopole  which  was  complementary  to  the 
impedance  of  the  slot:  top  loading  of  the  monopole  to  reduce  its  resonant 
impedance;  and  increasing  the  feedline  impedance  to  85  ohms,  a  value 
closer  to  the  300  ohm  resonant  impedance  of  the  slot. 

Since  it  was  easier  to  adjust  the  monopole  length  without  the  top 
load,  the  latter  method  was  used  in  construction  of  the  array.  The 
impedahce  characteristics  of  this  monopole-slot  are  shown  in  Figure  2. 
Note  in  this  figure  that  the  VSWR  is  less  than  1.5  for  frequencies  below 
1000  MHz  and  that  nearly  all  of  the  incident  power  is  radiated  between 
1100  and  1200  MHz. 

3.  DERIVATION  OF  CAVITY  AND  ARRAY  PARAMETERS 

The  basic  cell  size  and  scaling  factor  had  to  be  determined  before 
construction  of  the  log-periodic  array  could  begin.  The  cell  size  of  an 
array  of  cavity -backed  slots  is  determined  by  the  length-to-wid;h  ratio  of 
the  largest  element.  It  has  been  previously  found  by  Mayes1®  that  a  wider 
slot  results  in  a  significantly  lower  resonant  frequency  for  a  given  slot 
length  than  does  a  narrow  slot.  To  minimize  the  energy  coupled  into  the 
non-radiating  modes  it  is  desirable  to  have  the  cavity-backed  slot  antenna 
resonate  below  the  lowest  order  mode  (TEjj)  of  the  closed  cavity.  The 
cavity-backed  slot  previously  mentioned  (cavity  dimensions  10x3x.l25  in., 
slot  length  of  9  in.)  was  found  to  satisfy  this  requirement  (l  180  MHz  for  the 
slot  vs.  1256  MHz  for  the  cavity  mode). 

There  are  two  fundamental  opposing  factors  to  consider  when 
choosing  the  scaling  factor  for  an  array: 
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1.  Physical  limitations  due  to  size  limitations  and  material  cost 
and  availability 

2.  Operational  characteristics  such  as  bandwidth,  sidelobes,  etc. 

To  meet  the  physical  requirements  the  scaling  factor  should  be  as  small  as 
possible.  To  meet  the  electrical  performance  requirements  the  scaling 
factor  should  be  as  large  as  possible  to  increase  the  number  of  elements 
operating  in  the  active  region  of  the  array.  For  the  array  constructed  a 
compromise  scaling  factor  of  .89  was  chosen.  This  allowed  for  an  array 
with  a  "structural  bandwidth"  of  4.05:1  and  for  a  nearly  uniform 
distribution  of  power  between  3  elements  for  most  of  the  operating 
frequency  range. 

4.  CONSTRUCTION  OF  THE  ARRAY 

The  cavity-backed  slot  array  was  constructed  using  the  scaling 
factor  of  .89  and  a  base  cavity  size  of  1  1.5x3.375  in.  on  a  single  .125  inch 
thick  double  copper-clad  Rexolite  2200  substrate  of  dimension  12x24  inch 
(Figure  3).  The  cavity  walls  were  constructed  using  brass  screws  spaced  .5 
inches  apa'-t  for  the  largest  cavity  and  (.5  inch)  (.89)n-^  for  the  remaining 
elements.  It  has  been  previously  shown  that  this  arrangement  of  screws  is 
a  fairly  good  approximation  for  a  PEC  wall  for  the  frequencies  of 
interest.  The  Rexolite  substrate  was  mounted  on  a  .125  inch  thick 
aluminum  plate  for  increased  structural  stability.  The  slots  were  formed 
by  removal  of  the  copper  cladding  in  the  desired  area  on  the  top  surface  of 
the  substrate.  The  above  array  of  cavity-backed  slots  was  constructed  2-3 
elements  at  a  time  starting  with  the  low  frequency  elements  and  moving  to 
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higher  frequencies  with  each  succeeding  cavity.  (Cavity  and  slot 
dimensions  are  given  in  Table  1). 


TABLE  1 

DIMENSIONS  OF  CAVITY-BACKED  SLOTS 
(IN  INCHES) 


CAVITY 

NUMBER 

CAVITY 

LENGTH 

WIDTH 

SLOT 

LENGTH 

WIDTH 

1 

11.500 

3.375 

9.750 

.200 

2 

10.235 

3.003 

8.678 

.178 

3 

9.109 

2.673 

7.723 

.138 

4 

8.107 

2.379 

6.873 

.141 

5 

7.215 

2.118 

6.117 

.125 

6 

6.422 

1.885 

5.444 

.112 

7 

5.175 

1.677 

4.846 

,099 

8 

5.087 

1.498 

4.313 

.088 

Using  a  Hewlett-Packard  8510  vector  network  analyzer  each  cavity 
was  matched  to  a  rnor.opole  of  appropriate  length.  The  technique  used  is  as 
follows: 

1,  The  resonant  frequency  of  the  slot  was  determined  using  the 
network  analyzer  and  one  waveiength  50-85  ohm  microstrip 
tapered  ieedlines. 

2,  A  monopole  was  cut  to  a  quarter  of  the  free  space  wavelength 
corresponding  to  this  frequency. 

3,  T!ie  length  of  the  monopole  was  trimmed  to  minimize  the  input 
reflection  coefficient. 
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S-paraineters  were  taken  for  each  of  the  8  monopole-slot  elements. 
A  typical  set  of  S-parameters  for  a  mcnopole-slot  (the  fourth  in  the  array) 
is  given  in  Figure  5.  It  can  be  seen  from  the  S-parameters  shown  in  Figure 
5  that  the  monopole-slot  is  well  suited  for  use  as  an  array  element.  At 
frequencies  below  resonance,  the  element  is  welt  matched  to  the  feedline 
and  transfers  most  of  its  input  energy  with  a  shift  in  phase  to  the  next 
element  in  the  array.  As  the  frequency  is  increased  to  resonance,  the 
element  radiates  most  of  its  input  energy  (represented  by  the  simultaneous 
low  magnitudes  of  S2[  and  Sj  j). 

The  lengths  of  the  matching  monopoles  (3/32  inch  brass  rod  with  // 0- 
72  threads  on  the  bottom  to  fasten  to  the  nuts  soldered  to  the  feedline)  are 
given  in  Table  2. 


TABLE  2 

LENGTHS  OF  MONOPOLES 


NUMBER  LENGTHS  (cm) 


1 

2 

3 

4 

5 

6 

7 

8 


2.87 

3.46 

3.80 

4.07 

4.44 

4.96 

5.49 

6.38 


The  final  phase  of  array  construction  was  building  the  meandering 
feediines  and  the  two  corresponding  85-50  ohrn  linear  tapers.  Three 
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feedlines  with  slightly  different  interelement  spacing  were  constructed  to 
allow  for  experimental  uncertainties  in  the  individual  element  responses. 
Feedline  1  was  constructed  for  proper  phasing  for  the  backfire  radiating 
condition.  Feedline  2  was  7.5%  longer  between  feedpoints  while  feedline  3 
was  7.5%  shorter  between  feedpoints.  The  required  spacing  for  the 


backfire  condition  is  given  by  Equation  1. 


d  = 


( 


-/ 


(  eef  f  ) 


) 

.  5 


(1) 


where  Xq  =  free  space  wavelength  corresponding  to  the  resonant  frequency 
of  the  larger  element  (cm) 

l =  the  physical  separation  between  the  elements  (cm) 
e ef f  =  effective  permittivity  of  the  microstrip  feedline  (  = 

2.048  for  this  85  line) 

Equation  (1)  was  solved  for  the  section  of  line  between  the  first  and  second 
largest  elements.  For  this  condition  of  line  Aq  =  27.78  cm  and  8.10 
cm.  Using  Equation  (1)  results  in  the  following: 

d  =  1 3.75  cm  =  5.41 3  in  (feedline  1) 

d+  =  d  +  (Q.075)d  =  5.819  in.  (feedline  2) 

d"  =  d  -  (0. 07 5)d  =  5.006  in  (feedline  3) 

The  resulting  lengths  of  feedlines  between  elements  are  shown  in  Table  3. 

The  three  feedlines  were  constructed  on  .125  inch  copper-clad 
Rexolite  2200  substrates  using  the  values  given  in  Table  3.  #0-72  nuts  were 
soldered  to  the  feedline  above  the  point  where  the  feedline  crossed  over 
the  slot  to  allow  for  connection  of  the  monopoles  to  the  feedline.  Two  85- 
50  ohm  tapers  to  match  the  85  ohm  feedlines  to  the  50  ohm  measurement 
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system  were  also  constructed  on  .125  inch  Rexolite  substrate.  Using  a 
time-domain  re  Hectometer,  the  whole  feedline  assembly  was  trimmed  for 
minimum  discontinuities  at  the  connectors  and  feedline  corners  for  each  of 


TABLE  3 

LENGTH  OF  FEEDLINE  BETWEEN  ELEMENTS 
(IN  INCHES) 


FEEDLINE  1 

FEEDLINE  2 

FEEDLINE  3 

SECT. 

d 

EXCES 

3  d+ 

EXCESS 

d- 

EXCESS 

1 

5.413 

2.223 

5.819 

2.629 

5.007 

1.818 

2 

4.817 

1.979 

5.179 

2.340 

4.456 

1.618 

3 

4.287 

1.761 

4.609 

2.083 

3.965 

1.440 

4 

3.816 

1.567 

4.102 

1.854 

3.530 

1.281 

5 

3.396 

1.395 

3.651 

1.650 

2.824 

1.140 

6 

3.022 

1.240 

3.249 

1.468 

2,513 

1.015 

7 

2.690 

1.105 

2.892 

1.307 

2.237 

.903 

three  meandering  feedlines.  After  the  lines  were  trimmed,  the  bottom 
ground  plane  of  the  feedlines  was  removed  in  the  area  of  the  feedlines 
above  the  slots  to  allow  for  coupling  of  energy  to  the  slots.  See  Figure  6 


for  the  arrangement  of  the  array  complete  with  feedlines. 


5.  ELECTRICAL  MEASUREMENTS 

The  two  basic  types  of  measurements  made  on  the  complete  array 
were  S-parameter  measurements  and  E-field  pattern  measurements. 

The  S-parameter  measurements  were  made  with  the  Hewlett- 
Packard  8510  network  analyzer  using  feedline  3.  The  input  VSWR  for  the 
array  is  shown  from  500  MHz  to  4  GHz  in  Figure  7.  Sj  j  is  shown  in  Figure 
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8.  The  magnitude  of  the  transmission  coefficient  is  shown  from  500  M Hi: 
to  2.5  GHz  in  Figure  9. 


It  is  seen  from  Figure  7  that  the  input  VSWR  of  the  antenna  is  less 
than  2  for  all  frequencies  below  1950  MHz  and  is  less  than  3  for  all 
frequencies  below  3  GHz.  Thus  it  is  seen  that  the  structural  stopband  in 
regions  observed  for  other  series  fed  log-periodic  slot  arrays  have  been 
effectively  eliminated  by  using  the  monopole-slot  element.  This  antenna, 
therefore,  represents  a  stable  load  *or  most  sources.  Note  also  in  Figure  8 
the  periodic  nature  of  the  impedance  as  expected  for  a  log-periodic  array. 
Looking  at  Figure  9  it  is  seen  that  above  1  GHz  very  little  energy  reaches 
the  second  port.  This  and  the  low  input  VSWR  of  the  array  below  1980  MHz 
show  that  the  antenna  radiates  nearly  all  of  the  energy  presented  to  its 
input.  Since  very  little  energy  is  left  on  the  feedline  after  passing  through 
the  array  few  if  any  large-end  truncation  effects  should  be  observed  when 
using  this  array  above  1  GHz. 

E-field  measurements  were  made  on  20x20  foot  ground  plane  pattern 
range.  For  all  measurements  the  antenna  under  test  was  used  as  the 
transmitting  antenna.  For  the  azimuth  patterns  a  vertically  polarized  log- 
periodic  monopole  array  was  used  for  the  receiving  antenna.  Elevation 
patterns  were  taken  using  a  dipole  mounted  in  a  corner  reflector  attached 
to  a  fiberglass  boom  and  positioned  so  that  it  was  vertically  polarized.  The 
patterns  measured  for  feedline  3  are  shown  in  Figures  10  and  1 1. 

Half-power  beamwidths  and  front-to-back  ratios  were  found  from 
the  patterns  and  are  given  in  Figures  12  and  13  respectively.  In 
computation  of  the  front-to-back  ratios,  the  size  of  the  backlobe  was 
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assumed  to  be  that  of  the  largest  lobe  in  the  back  plane  defined  by  ISO 
±40°. 

From  the  patterns  it  is  seen  that  the  front-to-back  ratio  never 
deteriorates  significantly  over  the  whole  range  of  frequencies  measured. 
This  is  not  a  surprising  result  since  the  monopole-slot  has  a  fairly  high 
front-to-back  ratio  as  an  individual  element  (see  Figure  14).  Therefore,  a 
significant  degradation  of  the  front-to-back  ratio  cannot  be  used  as  a  valid 
limiting  factor  in  determining  the  pattern  bandwidth  for  this  array. 

The  half-power  beamwidth  for  this  antenna  is  relatively  stable  for 
frequencies  below  2.5  GHz.  Above  2.5  GHz  the  pattern  "balloons"  out  to  be 
a  nearly  180°  sector  pattern.  Thus  it  is  seen  that  the  upper  frequency  limit 
for  this  array  is  approximately  2.5  GHz.  However,  if  a  higher  directivity 
and  main  beam  on  axis  are  desired  the  upper  frequency  limit  is 
approximately  i.7  GHz.  Even  within  this  frequency  range  the  half-power 
beamwidth  does  vary  as  a  function  of  frequency  (Figure  12).  Some  of  this 
variaton  is  due  to  unequal  power  distribution  among  elements.  Figure  15 
and  Table  4  show  the  power  distribution  of  each  of  the  elements  as  a 
function  of  frquency.  From  the  figure  and  the  Table  it  is  seen  that 
generally  the  narrowest  beamwidths  were  measured  at  frequencies  where 
at  least  3  elements  were  each  radiating  at  leat  20%  of  the  total  power  and 
that  the  distribution  between  these  elements  was  fairly  uniform.  This  is 
not  surprising  since  element  power  "hogging"  reduces  the  effectiveness  of 
using  the  elements  in  an  array  environment. 

It  is  informative  to  look  at  the  power  distribution  graph  shown  in 
Figure  15  and  the  data  in  Table  4  to  define  limits  on  the  active  region  of 
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the  array.  The  most  active  element  is  seen  to  shift  as  the  frequency  is 
changed.  As  expected  at  high  frequencies,  only  the  first  two  or  three 
elements  radiate  a  significant  fraction  of  the  power  (notice  the  highly 


peaked  response  as  the  frequency  is  increased).  Looking  at  Table  4  for 

TABLE  4 

POWER  DISTRIBUTION  AMONG  ELEMENTS  OF  ARRAY 


FREQ. 

FRACTION  OF  INPUT  POWER  RADIATED  BY  EACH  ELEMENT 

Mhz 

1 

2 

3 

4 

5 

6 

7 

8 

500.0 

0.000 

.003 

.024 

.007 

.003 

.067 

.014 

.014 

600.0 

0.000 

0.000 

0.000 

0.000 

0.000 

.039 

.039 

.036 

700.0 

.013 

.007 

.005 

.006 

.034 

.039 

.042 

.042 

800.0 

.008 

.008 

.010 

.003 

.018 

.065 

.078 

.111 

900.0 

0.000 

0.000 

0.000 

0.000 

0.000 

.084 

.124 

.221 

1000.0 

0.000 

.007 

.004 

.012 

.055 

.090 

.223 

.351 

1100.0 

.016 

.026 

.040 

.044 

.080 

.155 

.289 

.349 

1200.0 

.010 

.033 

.076 

.173 

.125 

.243 

.331 

.006 

1300.0 

.012 

.043 

.079 

.132 

.227 

.344 

.142 

.009 

1335.0 

.019 

.058 

.085 

.145 

.259 

.362 

.054 

.006 

1415.0 

.056 

.080 

.141 

.198 

.379 

.137 

.004 

.001 

1456.0 

.125 

.091 

.184 

.213 

.348 

.033 

.003 

.001 

1500.0 

.084 

.112 

.230 

.253 

.320 

.001 

.000 

.000 

1544.0 

.097 

.162 

.245 

.286 

.200 

.006 

.002 

.002 

1590.0 

.103 

.215 

.268 

.305 

.092 

.010 

.003 

.001 

1685.0 

.150 

.269 

.355 

.224 

.001 

.000 

.000 

.000 

1700.0 

.188 

.308 

.351 

.153 

.000 

.000 

.000 

.000 

1800.0 

.207 

.461 

.322 

.007 

.001 

.000 

.000 

.000 

1900.0 

.333 

.589 

.071 

.004 

.001 

.001 

.000 

.000 

2000.0 

.475 

.517 

.006 

.001 

.000 

.000 

.000 

.ooo 

2100.0 

.793 

.174 

.018 

.006 

.007 

.001 

.000 

.002 

2200.0 

.992 

.005 

.001 

.001 

.001 

.000 

.000 

.000 

1200.0 

.925 

.037 

.014 

.009 

.004 

.002 

.002 

.001 

2400.0 

.732 

.111 

.059 

.068 

.005 

.002 

.019 

.000 

2500.0 

.652 

.131 

.083 

.083 

.009 

.008 

.005 

.004 

frequencies  above 

1.7  GHz.,  it  is 

seen  that 

one  element 

is  radiat' 

ng  well 

over  40%  of  the  power  incident  on  the  array.  Therefore,  one  would  expect 
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to  see  the  pattern  shapes  degrade  abovel.7  GHz.  Looking  at  the  patterns 
given  in  figure  10,  this  is  seen  to  be  the  case.  At  frequencies  near  the 
middle  of  the  operating  band  (1.2  GHz  -  1.7  GHz),  the  power  distribution  is 
more  uniform  than  at  higher  or  lower  frequencies.  Consequently,  one 
would  expect  to  see  well  defined  patterns  in  this  range.  Looking  at  the 
patterns  given  in  Figure  10,  it  is  seen  that  this  is  also  the  case.  At  low 
frequencies  the  higher  frequency  elements  do  not  radiate  a  significant 
fraction  of  the  power  (as  expected).  Hence,  when  frequencies  below  the 
resonant  frequency  of  the  largest  element  are  approached,  little  power  is 
radiated  by  any  of  the  elements.  Thus,  the  radiation  efficiency  of  the 
whole  array  is  decreased  significantly.  However,  as  previously  mentioned, 
well-formed  beams  are  still  observed  due  to  the  relatively  good  patterns 
for  each  element  when  the  element  is  operated  below  its  resonance.  Thus, 
it  is  seen  that  the  patterns  for  this  array  degrade  gracefully  at  frequencies 
beyond  both  ends  of  its  operating  bandwidth. 


6.0  CONCLUSION 

In  this  project  a  log-periodic  array  of  monopole-slot  elements  was 
constructed  and  tested.  It  was  found  that  using  the  monopole-slot  antenna 
in  the  array  does  indeed  eliminate  the  structural  stopband  regions  normally 
associated  with  a  series  fed  log-periodic  array.  In  fact,  the  antenna  had  an 
input  VSWR  below  2:1  for  most  of  its  operating  range.  Thus,  an  antenna 
constructed  similar  to  the  one  built  in  this  project  should  not  present 
loading  problems  for  most  sources.  The  front-to-back  ratio  was  found  to  be 
fairly  high  for  the  total  operating  range  of  the  antenna.  The  half-power 
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bearnwidth  and  pattern  shape  were  found  to  lx?  the  limiting  factors  in 
determination  of  the  upper  frequency  limit  for  this  array.  Radiation 
efficiency  for  the  array  was  found  to  establish  the  low  frequenry  limit. 
Using  these  two  limitations,  the  array  was  found  to  operate  quite 
effecively  between  1100  MHz  -  1800  MHz.  If  a  wider  half-power 
bearnwidth  is  acceptable,  the  high  frequency  limit  can  be  extended  to 
approximately  2500  MHz. 

Use  was  also  made  of  the  wider  element  bandwidth  inherent  with  the 
monopole-slot  in  lowering  the  value  of  the  scaling  factor  required  for 
proper  operation  of  the  array.  This  lowering  of  the  scaling  factor  allows 
for  use  of  physically  smaller  arrays  to  cover  the  same  frequency  range.  If 
narrower  and  more  uniform  beams  are  desired,  the  scaling  factor  should  be 
increased  from  the  value  of  0.89  used  in  this  array  to  allow  for  a  more 
uniform  distribution  of  power  between  elements.  If  a  wider  operating 
bandwidth  is  desired,  this  can  be  easily  accomplished  by  adding  more 
elements  to  the  array. 
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ABSTRACT 


A  Shallow-Cavity  Unity  Gain  Notch  Radiator 
by:  George  J.  Monser 

This  paper  presents  the  results  of  an  empirical  approach  for 
designing  shallow-cavity  (0.080  wavelengths  deep)  notch  radiating 
elements  for  use  in  greater  than  octave  bandwidth  applications. 
Three  different  size  models  were  built  and  evaluated.  It  was 
concluded  that  nearly  unity  gain  could  be  achieved  over  an  octave 
with  excellent  pattern  characteristics  over  greater  than  a  3/1 
bandwidth. 
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/  A  SHALLOW-CAVITY  UNITY  GAIN  NOTCH  RADIATOR 

By:  George  J.  Monser 

SUMMARY 

This  paper  describes  a  shallow-cavity  (0.080  wavelengths  deep) 
notch  radiating  elemenf^with  a  peak  gain  of  unity,  operable  over 
greater  than  one  octave  in  bandwidth.  Based  upon  three  different 
size  $ntennas  excellent  pattern  characteristics  were  measured 
over  a  3/1  frequency  band  when  flush-mounted  in  a  small  ground 
plane.  _ 

1.0  INTRODUCTION 

The  notch  radiating  element  is  a  wide  bandwidth  efficient 
antenna  which  has  been  used  in  many  applications  since  first 
reported  by  Johnson  [1],  These  applications  have  used  Johnson's 
notch  (or  nitch)  in  generic  form.  That  is,  as  a  simple  radiating 
discontinuity.  This  paper  extends  the  technology  to  radiating 
notches  over  shallow  cavities  where  the  plane  defining  the  notch 
is  parallel  to  the  cavity  base.  Efficient  radiating  structures 
with  cavity  depths  of  a  few  hundred's  of  wavelengths  are  attain¬ 
able  offering  an  alternative  to  patch  antennas  [2]. 

2.0  ANTENNA  CONFIGURATION 

Figure  1  shows  the  basic  elements  of  the  notch/cavity  anten¬ 
na.  In  this  figure  the  input  micro-strip  line  couples  energy  to 

*  Patented  by  author  and  Raytheon  Company. 
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the  slot  which,  in  turn,  excites  the  notch.  Without  the  cavity 
in  place  radiation  would  occur  with  a  maximum  along  the  axis  of 
the  slot.  Adding  the  cavity  changes  the  radiation  patterns, 
yielding  a  maximum  normal  to  the  plane  containing  the  slot  and 
notch.  Two  elements  contribute  to  the  far  field  pattern:  the 
microstrip  line  and  slot/notch. 

3.0  TEST  MODELS 

Three  models  were  built  for  evaluation: 

Model  A  1.3  x  1.5  x  0.29  inches. 

Model  B  2.7  x  3.0  x  0.50  inches. 

Model  C  3.1  x  3.8  x  0.62  inches. 

The  first  two  values  pertain  to  aperture  size.  The  last  value 

pertains  to  the  cavity  depth. 

Model  A  was  tested  and  evaluated  over  the  2.0  to  5.0  GHz 
range.  Model  B  was  evaluated  over  1.0  to  3.0  GHz  range.  Model  C 
was  evaluated  over  the  0.75  to  2.0  GHz  band. 

Because  of  space  constraints  the  cavity  depths  were  not 
scaled  in  proportion  to  notch  sizes. 

A  photograph  of  the  Model  A  antenna  is  shown  in  Figure  2. 

4.0  TEST  RESULTS 

Figures  3  through  6  show  typical  pattern  characteristics. 
Broad,  smooth  patterns  were  measured  in  both  planes  (i.e.,  along 
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and  normal  to  the  mlcrostrip  line)  with  the  broader  patterns  in 
the  plane  of  the  mlcrostrip  line. 

Figure  7  shows  gain  versus  cavity  depth  expressed  in  fractions 
of  a  wavelength.  On  the  order  of  unity  gain  is  displayed  over 
approximately  one  octave  commencing  with  a  cavity  depth  of  approx¬ 
imately-  0.06  wavelengths. 

Figure  8  shows  transmission  loss  versus  cavity  depth.  Again 
good  transmission  efficiency  (50  percent  or  better)  is  shown  over 
nearly  an  octave,  except  for  Model  B  with  the  deeper  cavity. 

5.0  DISCUSSION  OF  RESULTS 

Data  displayed  in  Figures  7  and  8  showed  similar  trends 
almost  independent  of  aperture  (slot/notch)  size.  By  computing 
directive  gain  (i.e.,  4itA/A2)  and  deducting  the  mismatch  loss, 
the  reported  gains  (Figure  7)  were  substantiated. 

One  anomaly  in  data  was  observed  in  Figure  8,  where  Model  B 
showed  a  narrower  transmission  band  for  50  percent  efficiency  than 
the  models.  This  occurrence  is  believed  to  be  attributable  to  the 
deeper  cavity.  Since  the  primary  effort  was  to  achieve  a  shallow 
cavity  design,  the  phenomenon  was  not  investigated. 
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6.0  CONCLUDING  REMARKS 


Based  upon  an  empirical  study,  it  is  felt  that  notches  over 
shallow  cavities  can  be  used  in  applications  requiring  broad  cov¬ 
erage  and  unity  gain  over  octave  bandwldths. 

If  the  restriction  of  unity  gain  is  removed,  then  the  full 
pattern  bandwidth  3/1  can  be  realized. 

As  configured,  this  type  of  antenna  offers  a  viable  solution 
for  flush-mounted  applications  requiring  very  little  depth  for 
installation. 
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Figure  3.  Model  B  Patterns 
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Figure  6.  Model  B  Patterns 
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ULTRA-BROADBAND  IMPEDANCE  MATCHING  USING 
ELECTRICALLY  SMALL  SELF-COMPLEMENTARY  STRUCTURES 


K.  G.  Schroeder 


American  Electronic  Laboratories,  Inc.  Lansdale,  Pa.  19446 


ABSTRACT 


M  A  number  of  self-complementary  structures  are  investigated 
with  the  goal  of  providing  frequency-independent  matching  of  the 
inductive  reactance  of  the  other  element.  Impedance  bandwidths 

of  up  to  50:1  are  found  to  be  achievable  with  very  small  (H  approxi- 

;  / 

mately  0.006  element  heights.  Efficiencies  are  determined  by 
achievable  loss  resistances  in  the  radiating  structure  itself, 
rather  than  the  (limited)  quality  of  an  external  loading  coil, 
and  by  energy  channeled  into  the  difference  port  of  the  feed  hybrid, 
which  is  used  to  match  two  complementary  impedances. 

1.0  INTRODUCTION 

The  ultra-broadband  self-complementary  pair  of  elements 
described  here  is  based  on  the  general  concept  first  formulated 
in  1963.  It  basically  uses  a  180-degree  hybrid  to  match  two 
radiators  which  are  closely  spaced  and  enhance  each  other's  radia¬ 
ting  properties  and  which  have  impedances  which  are  complementary. 

In  the  past,  mutual  coupling  between  closely  spaced  radiators  has 
prevented  the  physical  realization  of  such  a  pair,  except  for 
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th«  vertical  monopole/slot  combination  described  by  P.  Mayes,  which 
did  not  use  a  hybrid  feed.  The  evolution  of  a  monopole/half-loop 
combination  is  now  described,  which  both  overcomes  the  mutual 
coupling  problem  and,  at  the  same  time,  results  in  more  practical 
monopoie  structure  fc ■:  mobile  applications.  The  complementary  pair 
Impedance  matching  concept  is  now  discussed  in  more  detail. 

2.0  GENERAL  IMPEDANCE  MATCHING  PRINCIPLE 

The  complementary  matching  approach  is  based  on  the  impedance 
averaging  properties  of  a  magic-tee  or  180-degree  hybrid.  Fig.  1 
shows  the  generalized  impedance  relations  at  the  outputs  of  a  180- 
degree  hybrid  used  in  the  complementary  pair,  and  Fig.  2  shows  typical 
locations  of  complementary  impedances  on  a  Smith  Chart,  where  Z  and 
Z  Lt,  etc.  are  located  on  the  same  VSWR  circle,  but  on  opposite 
sides  with  respect  to  the  center  of  the  Smith  Chart.  It  can  be 
shown  that  the  sum-port  impedance  is  always  the  center  point  of  a 
straight  line  connecting  two  impedances  on  a  Smith  Chart,  even  if 
this  line  does  not  go  through  the  center,  i.e.  ,  even  if  the  two 
Impedances  are  not  exactly  complementary.  Physical  embodiments 
of  a  self-complementary  pair  of  small  elements  consist  of  loops 
(half-loops)  and  monopoles;  slots  and  monopoles  (dipoles);  dipoles 
and  folded  dipoles,  etc.  The  end-fire  complementary  pair  is 
formed  by  using  two  identical  elements  (e.g.  monopoles,  see  Fig.  3) 
and  externally  complementarizing  one  of  them  by  a  network,  such 
as  a  delay  line.  This  end-fire  pair  is  now  not  only  matched  with 
respect  to  the  self-impedances  of  the  two  elements,  but  also  with 
respect  to  mutual  impedances  in  a  phased  array  environment.  How¬ 
ever,  the  bandwidth  of  the  complementarizing  network  generally 
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limits  the  total  pair  bandwidth  to  about  4:1.*  For  applications 
where  10:1  total  bandwidth  or  more  is  required,  the  self-comple¬ 
mentary  approach  currently  offers  a  better  solution,  and  is 
described  below  in  more  detail. 

3.0  SELF-COMPLEMENTARY  MONOPOLE-PLUS-HALF-LOOP  PAIR 

Fig.  4  shows  a  basic  self-complementary  pair  consisting  of 
a  short  monopole  (whip)  and  a  grounded  half-loop.  One  is  an  open 
circuit  at  very  low  frequencies,  and  the  other  a  short  circuit. 
This  means  that  at  these  frequencies,  when  connected  Lo  a  feed 
hybrid  as  shown  in  Fig.  5,  the  sum-port  impedance  should  be 
matched. 

Fig.  6  shows  the  measured  sum-port  impedance  from  1.2  to 
115  MHz  for  a  whip  and  a  half-loop  with  0.002  wavelength  height 
and  stretched  length  (half-circumference),  respectively,  at 
2  MHz.  The  impedance  is  not  as  good  as  it  should  be  at  the  low 
end,  since  the  hybrid  impedance  deteriorates  below  approximately 
5  MHz,  and  the  ground  plane  was  too  small  for  frequencies  below 
approximately  5  MHz. 

Above  approximately  30  MHz  (model  frequency) ,  the  impedance 
match  was  suspected  not  to  be  as  good  as  possible  because  the  whip 
was  too  thin.  Remembering  that  a  fqt  conical  monopole  has  a  good 
high-frequency  performance,  one  might  try  to  improve  the  high- 
frequency  performance  by  using  a  conical  monopole.  However,  the 
mutual  coupling  between  the  loop  and  the  fat  monopole  will  then  be 


(*)  See  Companion  Paper:  "Feed  Designs  for  Broadband  End-Fire 
Complementary  Pairs  Using  Thin-Wire  Elements" 
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Coo  strong,  upsetting  the  self-complementarity.  Hence,  a  new  radia¬ 
tor  approach  was  evolved,  which  is  described  below,  and  which  is 
capable  of  optimizing  electrically  small,  self-complementary 
structures. 

4.0  EVOLUTION  OF  THE  SELF-COMPLEMENTARY  SCHROEDEk  ANTENNA 

Fig.  7  shows  the  well-known  equivalency  between  a  fat  cylin¬ 
drical  monopole  with  conical  feed  section  and  a  flat  sheet  of  a 
width  equal  to  twice  the  diameter  of  the  cylinder  diameter.  For 
the  frequencies  involved  in  most  HF  and  VHF  systems,  the  dimensions 
of  the  radiator  are  still  too  bulky  and  do  not  lend  themselves  to 
easy  installation  and  deployment.  The  next  step  in  simplifying 
the  structure  is  now  shown  in  Fig.  8.  Instead  of  using  a  continu¬ 
ous  metal  surface  or  sheet,  one  can  simulate  this  sheet  by  building 
an  outline  of  the  actual  radiating  structure.  This  can  be  done, 
because  the  radiating  currents  essentially  travel  up  the  extreme 
edges  of  the  radiator.  The  center  portion  does  have  some  effect, 
particularly  at  the  high  frequency  end,  by  virtue  of  the  capaci^^ve 
field  lines  connecting  to  the  surrounding  plane.  Because  of  this, 
the  optimum  angle  B  is  almost  zero  degrees.  That  is,  the  optimum 
equivalent  sheet  radiator  a  flat  (instead  of  angled)  bottom  is 
most  nearly  equivalent  to  the  optimized  conical  section  of  the 
cylindrical  conical  monopole,  where  the  included  cone  angle  is 
about  60  degrees  and  a  also  equals  60  degrees.  This  angle  can 
be  found  empirically  by  varying  a  until  the  impedance  plot  is 
most  ideally  centered  around  the  50-ohm  point  on  the  Smith  Chart 


for  the  widest  frequency  range.  This  new  equivalent-sheet,  wire- 
outline  monopole  can  now  be  used  in  combination  with  a  half-loop 
to  form  a  self-complementary  pair,  where  the  interaction  (through 
mutual  coupling)  between  the  loop  and  the  monopole  is  minimized 
by  the  symmetry  and  orthogonality  of  the  current-carrying  parts 
of  the  structure.  A  further  attempt  in  improving  symmetry  can  be 
made  by  using  four  whips  instead  of  two,  which  allows  the  use  of 
orthogonal  loops  also.  This  is  shown  in  Fig.  9,  except  that  only 
one  of  the  two  orthogonal  half-loops  is  shown.  The  reason  for 
the  second  loop  will  be  explained  below  when  discussing  patterns. 

5.0  TEST  RESULTS 

A  VSWR  plot  of  the  four-whip  arrangement  is  shown  in  Fig.  10 
for  a  small  pattern  model  built  at  a  50:1  scale  factor  so  that 
good  patterns  could  be  measured  svith  a  limited-size  ground  plane. 
The  impedance  match  is  obviously  excellent  (better  than  2:1)  over 
a  more  than  20:1  bandwidth.  From  100  to  1000  MHz  (full-scale 
equivalent  e.g.  3  to  30  MHz)  the  azimuth  patterns  (Fig.  11  to  16) 
show  astoundingly  good  front-to-back  ratios,  averaging  better  than 
10  dB.  The  explanation  for  this  directivity  is  the  fact  that, 
over  the  entire  band  where  the  naif-loop  is  electrically  shorter 
in  circumference  than  the  resonant  length,  the  electrical  current 
vector  of  the  half-loop  is  pointing  upward  at  the  feed  point,  thus 
enhancing  the  electrical  vector  of  the  monopole  which  is  also 
pointing  up.  The  electrical  vector  of  the  half-loop  points  down¬ 
ward  at  the  grounded  end,  and  therefore  cancels  out  the  monopole 
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contribution.  This  forms  a  directional  pattern  in  the  direction  of 
the  half-loop  feed  point.  This  performance  is  also  present  when 
only  two  whips  are  used,  which  are  orthogonal  in  their  feed  arrange¬ 
ment  (e.g.,  using  a  center  feed  bar  going  underneath  the  center 
of  the  half-loop).  With  four  whips,  two  orthogonal  loops  can  be 
installed  with  four  feed  points  that  are  alternately  switched  to 
the  hybrid  output  or  to  ground,  so  that  the  directional  pattern 
can  be  switched  into  four  major  directions  spaced  90  degrees  apart 
(0,  90,  180,  and  270  degrees). 

6.0  CONCLUSION  AND  SUMMARY 

A  new  embodiment  of  self-complementary  pairs  has  been  des¬ 
cribed  which  shows  very  wide  impedance  and  directional  pattern 
bandwidths.  This  antenna  has  also  a  small  physical  size  in  com¬ 
parison  to  product-line  HF  antennas.  The  "hybrid-fed,  wire-outline, 
equivalent-sheet-monopole  plus  half-loop,  electrically  small, 
self-complementary  pair",  (short:  self-complementary  Schroeder 
antenna)  warrants  further  detailed  investigations  in  the  form  of 
model  and  full-scale  experiments.  Already,  impedance  bandwidths 
in  excess  of  50:1  have  been  measured,  and  ultimate  bandwidths 
could  be  as  high  as  100:1,  with  the  only  limitation  being  given  by 
acceptable  efficiency. 
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Fig.  2:  Complementary  Impedances  on  the  Smith  Chart 
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plementary  Impedance  Matching  System  Using 
Half-Loop  (Element  Layout) 
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Fig.  5:  Self-Complementary  Impedance  Matching  System  Using 
Whip  and  Half-Loop  (Feed  Circuit) 


Fig.  6: 


Self-Complementary  Impedance  Matching  System  Using 
Whip  and  Half-Loop  (sum-Port  Input)  Impedance 
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Fig.  9:  Self-Complementary  Pair  of  Equivalent  Conical 
Monopole  Wire-Outline  Radiator  and  Fat  Strip 
Half-Loop 
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Fig.  10:  Impedance  Match  of  Patter,,  Model  for  Self-Complementary  Pair 


Fig.  11:  Azimuth  Pattern  ot  Self-Complementary  Pair 
of  Monopole  and  Half-l,oop,  100  MHz 


Fig.  14:  Azimuth  Pattern  of  Self-Complementary  Pair 
of  Monopole  and  Half-Loop,  600  MHz 
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Fig.  16:  Azimuth  Pattern  of  Self-Complementary  Pair  of 
Monopole  and  Half-Loop,  1000  MHz 
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IMPEDANCE-INVERTING  FEED  DESIGNS  FOR  BROADBAND 
ENDF1RE  COMPLEMENTARY  PAIRS  USING  THIN  WIRE  ELEMENTS 


K.  G.  Schroeder 


American  Electronic  Laboratories,  Inc.  Lansdale,  Pa.  19446 


ABSTRACT 

Broadband  impedance  matching  of  externally  coroplementarized 
endfire  monopole  pairs  is  described,  where  the  radiating  elements 
consist  of  thin  wires.  Broadband  impedance  loci  are  achieved  for 
the  thin  radiators  by  Impedance  inverters  placed  prior  to  the  com- 
plementarizing  circuit.  Impedance  matching  is  then  accomplished 
in  a  standard  broadband  hybrid  or  *Magic  Tee*. ^ 


A 


r 


y 
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1.0 


INTRODUCTION 


From  Ref.  1,  it  is  known  that  fat  monopoles,  arrayed  in 
endfire,  can  be  Impedance-  matched  over  at  least  a  3:1  frequency 
range.  (Fig.  1).  The  conical  monopoles  required  for  this  per¬ 
formance  typically  have  included  cone  angles  of  60  degrees  and 
length-to-diameter  ratios  of  less  than  approximately  4:1.  Solid 
monopole  complementary  pair  arrays,  however,  are  not  very  practical 
for  HF  and  VHF,  particularly  in  mobile  applications.  A  standard 
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way  to  build  conical  monopoles  at  those  frequencies  is  the  wire-cage 
construction. 

To  implement  a  quasi-short  monopole  (height  between  one- 
sixth  of  a  wavelength  and  a  half  a  wavelength)  in  wire-cage  con¬ 
struction,  typically  eight  to  six  wires  are  required,  which  shape 
the  outline  of  the  cone  and  attached  cylinder.  Keeping  in  mind  that 
this  provides  one  monopole,  a  second  wire-cage  structure  has  to  be 
added  to  the  first  one  to  form  a  pair.  This  then  leads  to  12  wires 
minimally  in  an  endfire  complementary  pair  of  conical  wire-cage  mono¬ 
poles.  (Fig.  2). 

These  wire-cage  structures  now  have  to  be  supported  somehow, 
which  leads  either  to  an  almost  impractical  complex  array  and  support 
structure,  or  the  need  for  12  individual  self-supporting  whips  forming 
the  wire-cage  outlines  in  the  verical  direction. 

The  desire  clearly  exists  to  further  reduce  the  number  of 
vertical  radiators.  It  is  known  that  the  conical  cylindrical  mono¬ 
poles  can  be  replaced  by  their  equivalent  flat  sheet  monopoles  pro¬ 
vided  the  width  of  the  sheet  is  twice  the  diameter  of  the  cylinder. 
(Fig.  3).  Figure  4  shows  an  endfire  pair  using  solid  sheets  of 
conducting  material.  (See  also  companion  paper  entitled:  "Ultra- 
Broadband  Impedance  Matching  Using  Electrically  Small  Self-Com¬ 
plementary  Structures".  ) 

Such  solid-sheet  "equivalent"  arrays  offer  no  advantage, 
of  course,  unless  the  sheets  are  also  configured  in  thin-wire 
construction  technique.  The  flat  sheet  monopoles,  however,  can 
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be  simulated  by  using  just  two  wires  or  two  whips,  thus  reducing 
the  total  number  of  elements  to  four.  Mo  connection  is  required 
between  the  tops  of  the  whips.  (Fig.  5)  Another  way  of  explaining 
the  horizontal  connecting  bar  of  the  wire-outline  equivalent  - 
sheet  radiator  is  that  it  constitutes  an  Impedance  inverter.  The 
thin-wire  monopole  Impedance  is  Inverted  about  the  characteristic 
Impedance  of  the  open-wire  transmission  line  leading  from  the  whip 
endpoint  to  the  center  (feed)  point. 

2.0  MODEL  DESIGN 

In  order  to  be  able  to  investigate  the  impedance  matching 
performance  of  the  equivalent-sheet  whip  array  concept,  an  impedance 
t£st  model  was  constructed  and  various  locations  on  a  ground  plane 
were  tested  using  an  end-fire  beam-forming  network.  The  elements 
were  located  on  a  square  and  on  a  rectangle.  The  full-scale  whips 
are  envisioned  to  be  standard  fiberglass  whips  with  conducting  braid 
embedded  in  fiberglass,  mounted  on  a  spring  and  base  Insulator. 

3.0  IMPEDANCE  MATCHING  TEST  RESULTS 

A  single  thin  monopole  (whip)  was  shown  to  be  poorly  matched, 
and  its  impedance  locus  on  the  Smith  Chart  generally  not  useful  for 
complementary  matching.  Fig.  6  shows  the  measured  impedance  plot 
for  a  thin  wire  with  about  0.12  \  height  at  the  lowest  frequency. 

The  wire  was  mounted  in  a  ground  plane  location  as  shown  in  Fig.  6. 

The  sum-port  impedance  for  a  rectangular  array  phased  into  the  front/ 
aft  direction  of  the  ground  plane  is  shown  in  Fig.  7.  This  was  one 
of  the  earlier  test  results,  and  was  aimed  at  proving  the  general 
concept.  The  basic  idea  is,  of  course,  that  two  of  the  whips  con- 
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nected  together,  as  shown  in  Fig.  5  will  behave  lmpedance-wisc 
like  a  fat  conical  mcnopole. 

Fig.  8  and  9  show  the  Input  Impedances  of  two  whips  v < th 
a  height  of  0.16  X  when  fed  in  parallel  with  impedance-inverting 
horizontal  members.  Fig.  6  Is  the  front  pair  and  Fig.  9  the  back 
pair.  Fig.  8  and  9  show  that  the  horizontal  feed  circuits  Indeed 
transform  the  center  of  the  whip  impedance  locus  towards  the 
center  of  the  Smith  Chart,  which  is  a  prerequisite  for  complementary 
matching.  The  sum-port  input  impedances  for  a  fore/aft  phasing 
condition  are  shown  in  Fig.  10.  It  is  seen  that  the  impedance 
match  at  the  very  low  end  had  deteriorated  somewhat,  the  3:1  VSWR 
circle  being  penetrated  above  a  frequency  of  approximately  220  MHz 
instead  of  200  MHz.  The  spacing  between  whips  was  a  constant 
0.12  X  at  the  lowest  frequency.  The  same  is  true  for  the  left/ 
right  phasing  condition  (Fig.  11),  the  two  impedance  plots  being 
very  similar.  This  implies  that  the  dimensions  between  the  whips 
are  more  critical  than  the  ground  plane  conditions  in  front  of  the 
forward  radie*-ing  element,  i.e„,  the  whole  array  Is  matched 
together  as  an  entity,  and  independent  of  ground  plane  extensions. 
The  patterns,  of  course,  may  show  more  gain  reduction  in  the  left/ 
right  direction. 

4.0  DIFFERENCE  PORT  POWER 

One  of  the  many  design  aspects  of  a  broadband  antenna 
system  is  the  matching  network  efficiency.  In  the  case  of  the 
complementary  pair,  there  is  an  immediate  measure  in  the  power  lost 
in  the  difference  port.  That  power  is  shown  versus  frequency  in 
Figures  12  through  16  •  As  one  would  expect,  a  significant  amount 
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Is  lost  in  the  difference  port  at  the  lowest  frequency*  where  much 
Impedance  matching  has  to  be  achieved  In  the  Inversion  and  comple- 
mentaritatlon  process.  A  secondary  peak  appears  between  700  and 
900  MHz.  In  both  cases*  however,  the  power  Is  down  more  than  5  dB 
from  the  Input*  meaning  that  the  matching  efficiency  Is  better  than 
502,  Obviously,  some  additional  losses  will  occur  and  absolute 
gain  measurements  will  have  to  be  made.  Preliminary  gain  measure¬ 
ments  have  confirmed  a  gain  of  not  less  than  0  dBi,  and  up  to  6  dBi. 

5.0  PATTERN  TEST  RESULTS 

Figures  17  through  20  show  representative  azimuth  patterns 
at  210,  400,  600,  and  800  MHz  model  frequencies  for  an  array  phased 
in  the  left/right  direction.  It  can  be  seen  that  even  at  the  lowest 
frequency  there  Is  a  f ront-to-back  (i.e.,  left  side-to-right  side) 
ratio  of  approximately  3  dB,  which  increases  to  6  dB  at  400  MHz, 

10  dB  at  600  MHz,  and  averages  better  than  10  dB  at  800  MHz.  The 
associated  elevation  patterns.  In  the  plane  of  the  main  beam,  are 
shown  in  Figures  21  through  24.  The  roll-off  varies  from  about 
3  dB  over  most  of  the  band  to  about  7  dB  at  the  highest  frequency. 
(The  peak  of  the  main  beam  ir>  pointing  to  the  right  side  of 
Figures  21  through  24.) 
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Fig.  5:  End-Fire  Complementary  Pair  of  Equivalent-Sheet 

Monopoles  Simulated  by  Dual  Whips  Fed  from  Impedance- 
Inverting  Feed  Bar 


Fig.  6:  Individual  Thin  Whip  Model  Impedance  Plot 
from  190  to  800  MHz 
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Fig.  7:  Matched  Sum-Port  Model 
with  Unequal  Spacing 


Impedance  for  Four-Whip  Array 
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Fig.  8:  Combined  Input  Impedance  of  Two  Whip  Antennas 
with  Horizontal  Impedance  Inverter 
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Fig.  9:  Combined  Trput  Impedance  of  Two  Whip  Antennas 
with  Horizontal  Impedance  Inverter 


Fig.  10:  Sum-Port  Input  Impedance  of  Four-Whip  Array 

System  with  Phasing  in  the  Fore/Aft  Direction 
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Fig.  11:  Sum-Port  Input  Impedance  of  Four-Whip 
Array  System  with  Phasing  in  the  Left/ 
Right  Direction 
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Position  (200  to  400  MHz) 


Fig.  13:  Difference  Port  Relative  Power  for  Broadside  Beam 

Position  (400  to  550  MHt) 


Fig.  14:  Difference  Port  Relative  Power  for  Broadside  Beam  Position 

(550  to  1000  MHz) 


Fig.  15:  Difference  Port  Relative  Power  for  Endfire  Beam  Position 

(200  to  550  MHz) 


Fig.  16:  Difference  Port  Relative  Power  for  Endfire  Beam  Position 

(550  to  1000  MHz) 


Fig.  17:  Azimuth  Pattern  of  Vehicular 


Array  (210  MHz) 
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Fig.  23:  Elevation  Pattern  of  Vehicular  Array  (600  MHz) 
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ABSTRACT 

This  paper  describes  the  design  and  measured  performance  of  very 
thick  cylindrical  monopole  antennas  (length-^to-diameter  ratios 
about  2).  These  very  thick  monopole  antennas  exploit  the 
physical  dimensions  of  common  coaxial  connectors  to  achieve  3:1 
frequency  bandwidths  with  VSWR  less  than  2:1  in  the  2  to  20  GHz 
frequency  range.  The  antenna  patterns  of  the  very  thick  monopole 
antennas  differ  significantly  from  those  of  their  thin 
counterparts.  The  body  of  revolution  moment  method  technique  is 
used  to  predict  the  E-plane  patterns  of  the  very  thick 
cylindrical  monopoles  as  well  as  those  of  conical  monopoles  at 
the  center  of  finite  circular  groundplanes .  Agreement  between 
measured  and  predicted  patterns  is  excellent. 

1  .  INTRODUCTION 

There  is  continuing  interest  in  finding  accurate  methods  for 
predicting  the  patterns  and  impedance  of  monopole  antennas  on 
finite  groundplanes  *  .  Both  the  moment  method-1  and  the 
geometrical  theory  of  diffraction  (GTD)^,1,,-’  have  been  used  to 
predict  the  patterns  of  thin  cylindrical  monopoles  at  the  center 
of  circular  groundplanes.  Generally  the  monopoles  are  about  a 
quarter  wavelength  long  which  considerably  simplifies  the 
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problem.  This  is  especially  true  for  the  GTO  where  the  dipole 
pattern  in  the  absence  of  the  groundplane  (the  geometrical  optics 
field)  must  be  known. 

Thin  monopoles  are  generally  narrow  band  antennas  resonating  at  a 
length  a  little  shorter  than  a  quarter  wavelength.  It  has  been 
known  for  a  very  long  time  that  the  impedance  variation  of  thick 
monopoles  is  significantly  less  than  that  of  thin  monopoles **• 
Thick  cylindrical  monopoles  have  been  studied  in  some  detail®’ ^ 
but  these  thick  mono poles  were  usually  driven  by  low  impedance 
coaxial  cables  (about  10  ohms). 

At  microwave  frequencies,  above  about  2  GHz,  it  becomes  possible 
to  make  very  thick  monopoles  (length-to-diameter  ratio  L/D  -  2) 
by  exploiting  the  dimensions  of  commonly  used  coaxial  connectors 
(N,  TNC ,  PTNC  and  SMA) .  Figures  1(a)  to  1(c)  show  the 
progression  to  a  very  thick  monopole  antenna.  The  diameter  of 
the  monopole  conductor  does  not  extend  beyond  the  dielectric  of 
the  50  ohm  coaxial  line.  This  significantly  reduces  the  excess 
base  capacitance 14  resulting  in  monopole  antennas  with  3:1 
frequency  bandwidth  for  a  VSWR  less  than  2:1.  Monopoles  with 
L/D  -  2  are  much  thicker  than  those  treated  in  11 

At  the  high-frequency  end,  these  very  thick  monopoles  are  more 
than  a  half  wavelength  long  and  their  patterns  differ 
significantly  from  those  of  their  thin  counterparts.  The 
patterns  of  the  conical  monopole  antenna  of  Figure  1(d)  can  be 
quite  complex  and  no  simple  solution  exists  although  extensive 
measured  data10  and  theoretical  results  for  a  60°  conical 
monopole11  exist. 
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This  paper  describes  the  design  and  VSWR  performance  of  very 
thick  monopole  antennas  having  3:1  bandwidth  when  fed  directly  by 
a  50  ohm  coaxial  line.  The  highly  efficient  body  of  revolution 
(BOR)  moment  method  formulation  of  Mautz  and  Harrington  is  used 
to  calculate  the  patterns  of  very  thick  cylindrical  as  well  as 
conical  monopole  antennas  on  finite  groundplanes.  In  most  cases 
the  groundplanes  are  quite  small  (one  to  three  wavelengths)  at 
the  low-frequency  limit  of  the  antennas.  The  predicted  patterns 
are  compared  to  measured  patterns  to  assess  the  accuracy  of  the 
technique . 

2.  DESIGN  AND  PERFORMANCE  OF  VERY  THICK  MONOPCLES 

2.1  Design 

To  exploit  the  reduction  in  impedance  variation  as  the  L/D  ratio 
of  monopoles  decreases,  several  monopoles  were  constructed. 
Figure  2  shows  the  typical  construction  of  the  antennas.  The 
antennas  are  all  based  on  standard  panel  mount  50  ohm  coaxial 
connectors.  The  thick  monopoles  (L/D  -  7)  are  made  by  cutting 
the  connector  center  conductor  to  the  desired  length  while  the 
very  thick  monopoles  (L/D  -  2)  are  made  by  adding  a  brass  bush  to 
the  center  conductor.  The  diameter  of  the  brass  bush  is  made  the 
same  as  that  of  the  dielectric  insulator  of  the  panel  mount 
connector.  All  the  monopoles  are  mounted  on  finite  circular 
groundplanes  whose  diameters  were  chosen  somewhat  arbitrarily  to 
lie  in  the  one  to  three  wavelength  range  at  the  low-frequency  end 
of  the  antenna  bandwidth. 
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Even  though  the  conductor  of  the  very  thick  monopole  does  not 
extend  beyond  the  connoctor  dielectric,  there  is  some  excess  base 
capacitance  which  must  be  compensated  for.  The  compensating  gap 
g  (see  Figure  2)  was  determined  experimentally  to  achieve  the 
best  VSWR  performance  over  a  nominal  3:1  bandwidth.  It  was 
observed  that  the  compensating  gap  increased  linearly  with  the 
diameter  of  the  very  thick  monopole.  Figure  3  shows  the 
compensating  gap  g  as  a  function  of  monopole  diameter  D.  This 
curve  can  be  used  to  estimate  the  gap  for  other  monopole 
diameters  but  it  applies  specifically  to  the  case  where  the 
monopole  and  connector  dielectric  have  the  same  diameter. 

Table  1  shows  the  dimensions  of  four  very  thicx  monopole  antennas 
based  on  the  SMA,  TNC,  PTNC  and  N  connectors.  The  symbols  are 
those  used  in  Figure  2. 

TABLE  1  s  Dimensions  and  frequency  limits  of  very  thick  monopole 
antennas . 


NO. 

TYPE* 

D 

(mm) 

L 

(mm) 

L/D 

d 

(mm) 

g 

(mm) 

W 

(mm) 

(mm) 

f 

rmin 

(GHz) 

f 

‘max 

(GHz) 

1 

SMA 

4.0 

7.6 

1.98 

1  25 

1.6 

129.0 

1.5 

■ 

21 .8 

2 

PTNC 

5.3 

11.2 

2.11 

1.62 

2.1 

76.0 

1.5 

| 

14.2 

3 

TNC 

7.0 

13.8 

1 .97 

2.14 

2.9 

129.0 

1.5 

4.1 

11.9 

4 

N 

10.0 

22. 1 

2.21 

3.00 

4.1 

175.0 

3.0 

n 

* 


Type  here  refers  to  the  standard  coaxial  connector. 


2.2  Measured  reflection  coefficient 


Figures  H  to  6  show  the  measured  reflection  coefficient  (VSWR) 
for  three  of  ’■he  very  thick  monopole  antennas.  Also  shown  is  the 
reflection  coefficient  for  thick  monopole  antennas  with  L/D  *  7. 
The  thick  monopole  antennas  show  the  usual  resonant  behaviour  of 
monopoles  while  the  very  thick  monopole  antennas  all  exhibit 
about  3:1  frequency  bandwidths  at  the  VSWR  -  2:1  points.  The 
frequency  limits  for  each  antenna  are  given  in  Table  1.  The 
superior  bandwidth  capability  of  the  very  thick  monopoles 
(L/D  *>  2)  over  the  thick  monopoles  (L/D  *  7)  is  clearly 
demonstrated  in  Figures  M  to  6.  Only  two  very  thick  monopole 
antennas  are  required  to  cover  the  2.6  to  21.8  GHz  frequency 
range. 


Comments  on  measured  patterns 


It  is  well  known  that  the  antenna  patterns  of  very  thick 
monopoles  do  not  differ  too  significantly  from  those  of  very  thin 
monopoles  when  the  lengths  are  about  0.25  X.  As  the  monopole 
length  increases  to  0.625  A  the  very  thick  monopole  patterns 
start  to  deviate  from  the  thin  patterns  by  showing  filled  in 
nulls  and  higher  secondary  lobes^. 


Figure  7  shows  a  set  of  three  measured  patterns  for  antenna 
number  1  (SMA  monopole  on  129  mm  groundplane)  at  10  GHz  where 
L/A  *  0.26.  The  three  antennas  have  L/D  -  1.95  (very  thick 
monoDole),  L/D  ■=  6.2M  (thick  monopole)  and  L/D  *  33.9  (thin 
monopole).  The  patterns  are  identical  (the  peak  value  at  60°  was 
adjusted  to  be  the  same  between  plots)  which  confirms  the 
observation  abou*  '  a  patterns  of  0.25  A  monopoles.  This  result 
indicates  that  the  assumption  of  using  the  thin  dipole  pattern  in 
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the  GTD  is  valid  even  for  very  thick  dipoles  provided  the 
monopole  height  is  about  X/M .  The  result  in  Figure  ?  should  be 
compared  to  that  in  Balanis^  where  the  GTD  plus  ring  currents  are 
used  to  solve  for  the  patterns  of  a  1/4  monopole  on  a  4  A 
groundplane.  Agreement  between  Figure  7  and  Figure  11.38  of 
reference  5  is  very  good. 

Antenna  number  4  (N-type  monopole  on  175  mm  groundplane)  was  used 
to  examine  the  effect  on  patterns  when  the  very  thick  (L/D  ■  2.2) 
and  the  thick  (L/D  -  7.4)  monopoles  are  0.625  A  long.  Figure  8 
shows  that  there  are  significant  differences  between  the  patterns 
in  the  upper  hemisphere  between  50°  and  90°.  In  the  lower 
hemisphere  the  sidelobe  structure  for  the  two  antennas  is 
identical  but  there  is  about  10  dB  difference  between  the  power 
levels.  The  peak  values  at  38°  were  adjusted  to  zero  dB;  the 
peak  value  of  the  ve^y  thick  monopole  lay  only  0.5  dB  above  that 
of  the  thick  monopole.  At  90°  the  levels  differ  by  6.5  dB.  Other 
pattern  measurements  show  that  differences  are  evident  when  the 
monopole  lengths  exceed  0.37  A. 

Because  the  monopole  antennas  are  very  light,  they  can  be 
supported  directly  on  a  0.141  inch  semi-'rigid  coaxial  cable  with 
a  swept  radius  90°  bend.  This  eliminates  the  need  for  test 
fixtures  and  structures  behind  the  antenna  groundplane.  The 
measured  patterns  are  extremely  symmetrical  and  it  is  felt  that 
the  backlobes  are  at  the  correct  levels  because  there  are  no 
obstructions  (except  for  the  short  horizontal  section  of  the  90° 
bend)  behind  the  antenna.  This  is  a  significant  improvement  over 
the  measurement  set-up  in^  where  the  back  mounting  structure 
influences  the  backlobes. 
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3. 


BOR  MODELLING  PROCEDURE 


A  body  of  revolution  is  a  body  whose  surface  is  generated  by 
rotating  a  plane  curve  around  an  axis.  This  leads  to  a  body  with 
perfect  rotational  symmetry.  Examples  of  bodies  of  revolution 
are  cylinders,  cones  and  spheres.  Mautz  and  Harrington12 
formulated  the  BOR  moment  method  technique.  This  is  a 
mathematical  procedure  which  exploits  the  symmetry  of  bodies  of 
revolution  by  using  modal  expansions  of  the  current  around  the 
BOR.  These  modes  are  chosen  to  be  orthogonal  and  are  thus 
uncoupled.  Each  mode  can  therefore  be  solved  independently  from 
other  modes.  This  leads  to  a  few  small  matrices  rather  than  one 
large  matrix.  Small  matrices  are  easier  and  faster  to  solve  than 
large  ones  and  this  gives  a  significant  saving  in  computer  time. 
For  a  completely  symmetric  aperture  as  considered  here,  only  one 
mode  needs  to  be  solved. 

To  model  an  antenna  using  the  BOR  technique  a  generating  curve 
must  be  defined.  Several  examples  of  such  generating  curves  are 
shown  in  Figure  9.  There  is  a  limitation  of  99  on  the  maximum 
number  of  points  that  can  be  handled  by  the  BOR  moment  method 
program  in  use.  This  limits  the  detail  that  can  be  modelled. 

Because  the  moment  method  calculates  the  current  on  the  body 
under  consideration,  it  is  very  important  to  put  segments  at  the 
correct  places  to  ensure  the  correct  current  flow  on  the  body. 
Where  large  changes  in  current  are  expected,  more  segments  must 
be  used  for  accurate  modelling.  Figure  9(a)  shows  that  more 
segments  are  used  on  the  monopole  than  on  the  groundplane.  On 
the  groundplane  more  segments  are  used  close  to  the  monopole  than 
further  away. 
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Figure  9(b)  shows  an  example  of  a  case  where  there  are  too  few 
segments  available  to  model  ill  the  details  of  the  antenna.  The 
lower  portion  of  the  groundplane  has  been  modelled  with  a  few 
segments  which  were  truncated  near  the  edge.  This  will  account 
for  the  radiation  from  the  finite  edge  of  the  groundplan.e. 

Figures  9(c),  9(d)  and  9(e)  show  three  models  of  the  conical 
monopole  of  Figure  1(d).  These  models  were  used  to  demonstrate 
different  effects  found  in  the  models.  Figure  9(c)  shows  the 
normal  model  of  the  capped  cone  which  models  only  the  upper  side 
of  the  groundplane  and  has  a  maximum  segment  length  of  2  mm.  This 
gives  a  segment  length  of  1/15  X  at  10  GHz.  Figure  9(d)  shows  a 
model  with  the  same  resolution  but  modelling  both  sides  of  the 
groundplane.  Figure  10  shows  the  difference  between  the 
predicted  patterns  for  the  conical  monopole  on  its  129  mm  ground- 
plane  using  these  two  models.  The  bottom  side  of  the  groundplane 
has  very  little  effect  on  the  computed  radiation  pattern  at 
10  GHz.  This  is  also  true  at  other  frequencies.  To  prove 
convergence  of  the  moment  method  solution,  a  high  resolution 
model  as  shown  in  Figure  9(e)  was  used.  This  model  has  a  segment 
length  of  1  mm  (1/30  X  at  10  GHz).  Figure  11  shows  the 
comparison  of  the  radiation  patterns  for  these  two  models.  There 
are  slight  differences  between  the  two  models  at  10  GHz  in  the 
±120°  region.  At  higher  frequencies  (18  GHz),  the  patterns  of  the 
two  models  start  to  differ  by  more  than  a  dB  showing  that  the 
normal  resolution  model  is  not  valid  at  these  frequencies  because 
there  are  too  few  segments. 

The  antenna  feed  network  is  modelled  realistically  by  using  the 
TEM  field  of  a  coaxial  cable.  This  is  done  by  defining  a 
constant  radial  field  from  the  center  conductor  to  the  outer 
conductor  of  the  coaxial  line  feeding  the  monopole. 
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An  attempt  was  made  to  compute  the  input  impedance  of  the 
antennas  at  different  frequencies  using  the  BOR  method.  This  was 
done  by  integrating  over  the  field  and  currents  in  the  feed 
aperture.  The  results  obtained  look  promising  but  more  work  is 
needed  on  this  aspect. 

4.  MEASURED  AND  PREDICTED  RESULTS 

The  VSWR  of  the  70°  capped  cone  of  Figure  1(d)  on  a  129  mm 
groundplane  (as  for  antenna  no.  1)  was  2.2:1  at  2  GHz  and  better 
than  1.5:1  from  2.5  to  above  21  GHz.  All  the  predicted  patterns 
were  computed  using  the  high  resolution  model  of  Figure  9(e). 

Figure  12  shows  the  measured  and  predicted  patterns  at  2  GHz 
where  the  groundplane  is  only  0.86  A  in  diameter.  Above  the 
groundplane  the  agreement  is  excellent,  below  the  groundplane  the 
measured  pattern  lies  about  2.5  dB  below  the  predicted  one. 
Because  the  groundplane  is  so  small,  currents  will  creep  onto  the 
underside  of  the  groundplane  and  possibly  excite  the  horizontal 
section  of  supporting  cable.  It  was  found  that  the  length  of  the 
horizontal  section  of  cable  did  influence  the  measured  patterns 
indicating  that  the  antenna  i3  not  completely  Isolated  from  the 
measurement  environment  at  2  GHz.  There  may  also  be  merit  in 
using  one  of  the  BOR  models  which  includes  the  lower  portion  of 
the  groundplane.  Figures  13  and  14  show  the  patterns  at  4  GHz 
and  10  GHz.  Now  the  agreement  is  excellent  both  as  regards  the 
mainlobes  above  the  groundplane  and  the  si  delobes  below  the 
groundplane. 
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Figure  15  shows  the  measured  pattern  for  antenna  no.  *1  at  2.7  GHz 
where  the  groundplane  is  1.58  \  in  diameter.  The  predictions 
were  made  using  the  model  of  Figure  9(b).  The  agreement  is  good 
except  near  the  horizon  where  the  measured  pattern  is  about 
2.5  dB  below  the  predicted  one.  It  is  felt  that  a  more  detailed 
model  of  the  underside  of  the  groundplane  and  a  means  to  suppress 
stray  currents  on  the  0.141  inch  semi-rigid  cable  could  improve 
this  comparison.  Figure  16  shows  excellent  agreement  at  7.9  GHz. 

Measured  and  predicted  patterns  for  antenna  no.  1  are  shown  for 
10  GHz  in  Figure  17.  Agreement  is  excellent.  The  patterns  for 
the  thick  cylindrical  monopole  in  Figure  17  should  be  compared  to 
those  of  the  70°  conical  monopole  in  Figure  14.  The  backlobe 
structure  of  the  two  sets  of  patterns  is  almost  identical  (as 
expected  since  the  shape  of  the  backlobe  structure  is  determined 
almost  entirely  by  the  diameter  of  the  groundplane).  However, 
there  are  significant  differences  in  the  upper  hemisphere  where 
the  maximum  of  the  conical  monopole  occurs  at  about  25°  while 
that  for  the  very  thick  monopole  occurs  at  60°.  This  emphasizes 

i 

the  fact  that  the  basic  radiator  pattern  can  strongly  influence 
the  pattern  structure  in  the  upper  hemisphere. 

5.  CONCLUSIONS 


Measured  reflection  coefficient  data  has  been  presented  to  show 
that  broadband  very  thick  monopoles  can  be  made  quite  simply  by 
exploiting  the  dimensions  of  common  connector  types.  Frequency 
bandwidths  of  3:1  can  be  achieved  with  VSWR  less  than  2:1. 
Because  the  broadband  monopoles  are  very  thick,  their  basic 
radiation  patterns  differ  significantly  from  those  of  their  thin 
counterparts  for  L/A  greater  than  about  0.4.  The  BOR  moment 
method  is  a  powerful  technique  for  taking  into  account  directly 
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the  radiation  characteristics  of  the  thick  monopoles.  In 
addition,  the  more  complex  conical  monopoles  can  also  be 
modelled.  Agreement  between  measured  and  predicted  patterns  is 
excellent,  the  main  differences  occurring  for  groundplane 
diameters  in  the  1  to  1.5  X  region.  Additional  work  is  required 
to  refine  the  BOR  model  and  the  measurement  technique  for  these 
very  small  groundplanes . 
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FIGURE  1  :  Various  monopole  antennas  (a)  thin  (L/D  large), 
(b)  thick  (L/D  «  7),  (c)  very  thick  (L/D  »  2)  and  70° 
conical  monopole  with  truncated  90°  cap. 


FIGURE  2  :  Schematic  diagram  showing  design  parameters  of  the 
monopole  antenna. 
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REFLECTION  COEFFICIENT  (dB) 


FIGURE 


FIGURE  3  :  Compensating  gap  vs.  monopole  diameter. 


£  :  Measured  reflection  coefficient  of  antenna  no.  i|  with 
L/D  *  2.21  ( - )  and  monopole  with  L/D  =  7.37  (— — — ) 
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REFLECTION  COEFFICIENT  (dB) 


FREQUENCY  (GH*) 


FIGURE  5  :  Measured  reflection  coefficient  of  antenna  no.  3  with 
L/D  -  1.97  ( - )  and  L/D  «  6.45  ( - — ) . 


FIGURE  6  :  Measured  reflection  coefficient  of  antenna  no.  1  with 
L/D  -  1.95  ( - )  and  L/D  -  6.24  (— — . 


FIGURE  7  :  Superposition  of  a  set  of  measured  patterns  at  10  GHz 
for  three  monopole  antennas  with  L/D  ~  1.95, 
L/D  -  6.24  and  L/D  -  33.9  on  a  129  mm  groundplane 
(antenna  no.  1  with  various  radiators). 
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FIGURE  8  :  Measured  patterns  at  8.^8  GHz  for  antenna  no.  ^  with 
monopole  radiator?  having  L/D  -  2.2  and  L/D  »  7.^. 
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FIGURE  9  :  Examples  of  BOR  generating  curves  showing  segmentation 
for  (a)  antenna  no.  1,  (b)  antenna  no.  4  and  the 

conical  monopcle  of  Figure  1(d)  with  (c)  normal  model, 
(d)  double  sided  groundplane  and  (e)  high  resolution 
model . 
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FIGURE  1 0  :  Difference  Letween  the  predicted  patterns  for  the 

normal  model  ( - )  and  the  double  sided  model  ( - ) 

of  the  conical  monopole  of  Figiire  1(d)  at  10  GHz. 
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FIGURE  1 1  :  Difference  between  the  predicted  patterns  for  the 

normal  model  ( - )  and  the  high  resolution  model 

( - )  of  the  conical  monopole  of  Figure  1(d)  at 

10  GHz. 
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FIGURE  12  :  Measured  ( 


)  and  predicted  ( - )  patterns  at 

2  GHz  for  the  70°  conical  monopole. 


0 


180 


FIGURE  13  :  Measured  ( - )  and  predicted  ( - )  patterns  at 

**  GHz  for  70°  conical  monopole  on  129  mm  groundplane. 
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FIGURE  1 4  :  Measured  ( - )  and  predicted  ( - )  patterns  at 

10  GHz  for  70°  conical  monopole  on  129  mm  ground- 
plane  . 


Measured  ( - )  and  predicted  ( - )  patterns  at 

2.7  GHz  for  antenna  no.  4  on  175  mm  groundplane. 
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FIGURE  17  :  Measured  ( - )  and  predicted  ( - )  patterns  at 

10  GHz  for  antenna  no.  1  on  129  mm  groundplane. 
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j  ABSTRACT 

Potential  applications  for  a  precision,  monopulse  polarimeter  exist 
in  many  areas,  from  ELINT  to  communications.  Many  applications 
require  coverage  of  a  broad  range  of  frequencies,  rapid  process¬ 
ing,  and  the  ability  to  present  measurement  results  to  both  human 
operators  and  other  digital  systems.  In  addition,  the  form  in 
which  measurement  results  are  presented  must  often  be  modified 
to  suit  a  particular  application.  The  need  for  a  powerful,  flexible, 
monopulse  polarimeter  prompted  the  development  of  the  Watkins- 
Johnson  Company  polarization  measurement  system.  ~¥ - - 

The  polarization  state  of  an  electromagnetic  wave  may  be  specified 
in  several  ways.  The  interpretation  of  polarization  information  in 
any  form,  however,  requires  an  understanding  of  vector  theory. 
In  order  to  visualize  polarization  relationships  the  cone* .^t  of  the 
Poincarg-  Sphere  can  be  very  useful.  From  the  Poincare  Sphere, 
the  polarization  box  may  be  derived  for  use  in  calculating  polari¬ 
zation  components.  A  discussion  of  these  concepts  lays  the 
groundwork  for  evaluation  of  polarization  measurement  techniques. 


POLARIZATION  MEASUREMENT 

David  L.  Hawthorne 
Watkins-Johnson  Company 
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The  Watkins-Johson  Company  polarization  measurement  system  dis¬ 
cussed  will  allow  monopulse  measurements  of  the  polarization  of 
received  signals  in  the  2  to  18  GHz  region.  The  system  utilizes  a 
spread  spectrum  multiplexing  technique  that  preserves  the  ampli¬ 
tude  and  phase  relationships  of  signals  from  two  orthogonal 
antennas.  The  ability  to  simultaneously  measure  the  relative  amp¬ 
litudes  and  phase  of  orthogonal  signals  without  the  detrimental 
effects  of  multiple  receivers  should  greatly  improve  the  accuracy 
of  polarization  measurements,  while  providing  the  flexibility  of 
broad  frequency  coverage.  In  addition,  the  processing  capability 
incorporated  in  the  system  will  allow  the  calculation  of  polarisation 
parameters  as  required  for  presentation  to  operators  or  reporting 
to  other  systems. 

1.0  INTRODUCTION 

The  area  of  microwave  signal  analysis  has  evolved  to  the  point 
where  systems  capable  of  detailed  characterizations  of  signals  are 
in  relatively  common  use.  Signal  parameters  that  are  routinely 
measured  include  frequency,  amplitude,  phase,  direction  of  prop¬ 
agation,  and  modulation  type  and  rate.  For  pulsed  signals,  prop¬ 
erties  such  as  pulse  width  an'j  PRI  are  frequently  determined.  In 
addition,  the  time  variance  of  signal  parameters  is  measured  as  it 
contains  information  about  scan  patterns,  frequency  agility,  dop- 
pler  and  stagger. 
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Careful  analysis  of  these  properties  using  digital  techniques  has 
produced  powerful  results,  but  no  system  can  exploit  the  full 
vector  nature  of  E.M.  waves  while  remaining  insensitive  to  polar¬ 
ization.  Unfortunately,  information  regarding  the  polarization  state 
of  an  electromagnetic  wave  has  traditionally  been  ignored,  dis¬ 
carded,  or  simply  unavailable  in  conventional  signal  analysis 
systems,  although  the  utilitarian  nature  of  such  information  has 
been  discussed  at  great  length. 

Applications  exist  in  several  fields  for  precision  polarization 
measurement.  The  most  obvious  use  for  such  equipment  is  in  the 
design  and  testing  of  antennas,  yet  polarimetry  has  the  potential 
to  enhance  system  performance  in  many  ways. 

As  the  demand  increases  for  improved,  passive,  surveillance 
capabilities,  systems  will  expand  to  provide  previously  unavailable 
measurement  data.  Since  polarization  diversity  and  agility  charac¬ 
teristics  can  be  used  to  distinguish  emitters,  ELINT  collection 
systems  taking  advantage  of  polarization  information  could  add 
another  dimension  to  the  array  of  parameters  used  for  signal 
classification.  Passive  ESM  systems  could  then  utilize  received 
polarization  data  for  purposes  of  emitter  identification,  or  use 
pulse  by  pulse  polarization  information  as  a  deinterleaving 
parameter. 
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PoelmarT  described  a  method  which  allows  improvement  in  target 

detection  performance  of  active  radars  in  clutter  and  interference 

environments  through  adaptive  antenna  polarizations,  when  return 

polarization  information  is  available.  Another  polarimetry  tech - 
2 

nique  looks  promising  for  improvements  in  target  identification 
using  radar  processing  that  examines  the  depolarization  of  back- 
scattered  radiation  to  obtain  geometrical  type  information  relat¬ 
ing  to  target  symmetry  and  orientation  in  space. 

ECM  systems  can  extend  their  effective  range  by  adapting  the 
polarization  of  transmitted  jamming  signals  to  match  that  of  the 
received  signal.  To  optimize  such  a  technique,  detailed  informa¬ 
tion  about  the  received  polarization  is  necessary. 

Earth-space  communication  links  arc  susceptible  to  cross-polarized 
signal  interference  resulting  from  depolarization  during  propaga¬ 
tion,  and  from  small  changes  in  orientation  of  satellite  antennas 
with  respect  to  ground  stations.  It  has  been  demonstrated  that 

cancellation  of  residuals  through  adaptive  antenna  polarization  is 
3 

desirable  ,  but  is  dependent  on  the  availability  of  polarization 
information  about  the  received  signal. 

The  use  of  polarization  processing  techniques  is  considered  the 
next  logical  step  in  the  enhancement  of  many  microwave  systems. 
However,  taking  the  step  depends  on  the  availability  of  equipment 
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that  is  capable  of  precision  polarization  measurements.  The  need 
for  such  equipment  prompted  the  development  of  Watkins-Johnson 
Company's  polarization  measurement  system. 

In  the  following  pages,  the  basic  equations  for  the  vector  repre¬ 
sentation  of  electromagnetic  waves  shall  be  discussed.  The  con¬ 
cepts  of  polarization  and  the  Poincare  Sphere  shall  be  presented 
in  order  to  understand  the  technique  used  to  measure  received 
polarization.  The  implementation  of  the  measurement  technique 
used  in  the  Watkins-Johnson  Company  polarization  measurement 
system  and  the  limitations  on  the  accuracy  of  such  systems  will  be 
examined. 

2.0  COORDINATE  SYSTEM  AND  VECTOR  REPRESENTATION 

Consider  a  rectangular,  cartesian,  coordinate  system  consisting  of 

three  mutually  orthogonal  axes  whose  variables  are  denoted  x,  y, 

and  z.  Unit  vectors  a  ,  a  ,  and  are  oriented  as  shown  in 

x  y  z 

Figure  1.  For  purposes  of  analysis  we  will  assume  that  we  are 

dealing  with  a  plane  wave  traveling  in  the  z  direction.  Phase 

fronts  of  the  wave  will  be  normal  to  a  . 

z 

A  plane  electromagnetic  wave  traveling  in  the  z  direction  is 

composed  of  electric  and  magnetic  time-varying  fields  that  lie  in 
4 

the  xy  plane.  The  fields  are  perpendicular  to  each  other,  and  at 
a  specific  time  can  be  represented  by  orthogonal  vectors  as  shown 
in  Figure  2. 
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Figure  1.  The  Rectangular  Coordinate  System 


Y 


Figure  2.  Electric  and  Magnetic  Fields  in  a  Plane  Wave 


The  E  and  H  components  of  the  plane  wave  vary  with  time  at  the 
same  frequency  and  in  the  same  phase,  and  the  magnitudes  are 
related  by  a  constant.  The  discussion  will  therefore  deal  only 
with  the  electric  field. 


The  superposition  principle  states  that  the  total  electric  field 

vector  of  a  wave  is  the  sum  of  all  electric  field  vectors  composing 
§ 

the  wave.  This  means  the  total  field  vector  for  a  wave  can  be 
decomposed  into,  or  constructed  from,  two  orthogonal  vector  com¬ 
ponents,  as  shown  in  Figure  3. 
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Figure  3.  Decomposition  of  Electric  Field  Vector  into 
Orthogonal  Components 

Component  vectors  E  and  E  are  related  to  the  total  field  vector 

x  y 

2  by  the  equation 

E(t)  =  Ex  ♦  Ey  =  Ex(t)ax  +  Ey(t)ay  (1) 

If  E(t)  is  a  function  of  a  single  frequency,  the  magnitudes  of  the 
orthogonal  components  can  be  expressed  as 

Ex(t)  =  Ex  cos  (cut) 

and 

E  (t)  =  E  cos  ( cut  +  5 

where  -180  deg  <6  <180  deg. 

3.0  POLARIZATION 
The  polarization  of  a  wave  is  related  to  the  orientation  of  the 
electric  field  vector  with  respect  to  the  coordinate  system  in  use. 

-553- 


(2) 

)  (3) 


aaoMiaamei5gaaMBBBa8Hw«a«aKBBiaoaBaaaHfltMBaBaBiaiiiaiaaiaaaaB5fl-'tf^aaflgai 


If  an  electric  field  vector  always  lies  in  a  given  plane  parallel  to 
the  direction  of  propagation,  the  wave  is  said  to  be  linearly 

Jfc  a 

polarised.  Since  E  and  E  fit  this  description,  the  orthogonal 

x  y 

components  we  have  defined  can  be  thought  of  as  linearly  polar¬ 
ised  fields  which  differ  in  time  phase  by  the  angle  8. 

In  general,  any  electric  field  can  be  resolved  into  orthogonal 

linear  components  of  appropriate  magnitudes  and  phase.  Combin- 

ing  E  and  E  to  form  the  total  field  ’  ector ,  and  plotting  the 
x  y 

locus  described  by  the  tip  of  the  resultant  vector  over  time,  will 

g 

generate  an  ellipse  (Figure  4). 


V 


Figure  4.  Polarization  Ellipse  Generated  from  the  Locus  of 
Points  at  the  Tip  of  E  over  Time 

The  electric  field  represented  in  Figure  4  is  said  to  be  elliptically 
polarized,  due  to  the  way  the  field  changes  with  time.  The  locus 
is  called  the  polarization  ellipse.  If  the  tip  of  the  electric  field 
vector  had  traced  out  a  circle,  the  polarization  would  be  called 
circular,  or,  as  mentioned  above,  if  the  vector  remains  at  a  con¬ 
stant  angle  with  respect  to  the  coordinate  system,  the  field  is 
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considered  linearly  polarized.  In  actuality,  all  polarizations 
can  be  considered  elliptical,  as  the  circular  and  linear  polariza¬ 
tions  are  simply  degenerate  (special)  cases  of  the  polarization 
ellipse. 

Discussions  of  polarization  commonly  refer  to  horizontal,  vertical, 

7 

slant-linear,  and  circular  polarizations.  These  terms  refer  to  the 
fact  that  the  coordinate  system  has  been  defined  in  relation  to  a 
particular  (usually  earth-bound)  frame  of  reference.  However, 
polarized  waves  are  rarely  of  a  purely  linear  or  circular  polari¬ 
zation,  due  to  depolarization  during  propagation,  as  a  result  of 
reflections,  or  due  to  imperfections  in  the  radiating  antenna.  Most 
waves  are  actually  elliptically  polarized.  However,  it  is  instructive 
to  consider  the  degenerate  cases  in  order  to  gain  an  understand¬ 
ing  of  the  relationships  between  the  orthogonal  components  for 
different  polarizations.  The  relative  amplitudes  and  phases  of  E 
and  Ey  for  these  cases  are  given  in  Table  1,  which  assumes  that 
the  coordinate  system  has  been  aligned  such  that  the  x  axis  is 
parallel  to  the  horizon. 

If  the  components  are  in  phase  but  have  unequal  amplitudes  the 
total  field  is  linearly  polarized  in  a  direction  that  makes  an  angle 
with  the  x  axis  that  is  given  by  the  equation 

angle  =  arctan  (|e  I/Ie  |)  (4) 

y  x 
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Table  1.  Amplitude  and  Phase  Relationships  between  Orthogonal 
Linear  Components  of  Degenerate  Polarisations 


Polarization 

Amplitudes 

Phase  (  8  ) 

Horizontal 

!eJ  -  o 

N/A 

Vertical 

IeJ  =  0 

N/A 

+45  Deg  Slant 

Kl  *  ^y1 

0  degrees 

-45  Deg  Slant 

II 

180  degrees 

Arbitrary  Linear 

Polarization 

Don’t  Care 

0  or  180  degrees 

Right  Hand  Circular 

IsU  =  lei 

*  i 

-90  degrees 

Left  Hand  Circular 

1**1  *  lfyl 

+90  degrees 

If  8  equals  ±90° ,  the  ellipse  is  oriented  with  its  major  and  minor 
axes  aligned  with  the  x  and  y  coordinate  axes.  Under  these 
conditions,  if  Ex  >  E^,  the  major  axis  is  2EX  and  the  minor  axis 
is  2E  .  If  E  =  E  the  ellipse  degenerates  into  a  circle.  If  phase 

y  *  y 

is  negative  (positive)  the  sense  of  rotation  is  clockwise  (counter¬ 
clockwise)  in  the  z  -  0  plane,  looking  in  the  direction  of  propa- 

g 

gation,  and  is  called  right-handed  (left-handed). 

4.0  AXIAL  RATIO,  TILT  ANGLE,  ELLIPTICITY  ANGLE 

The  ratio  of  the  major  axis  to  the  minor  axis  is  called  the  axial 

ratio  (r),  which  is  a  frequently  used  parameter  for  characterizing 


polarization.  The  second  parameter  used  in  conjunction  with  the 
axial  ratio  to  characterize  a  polarization  is  the  sense  of  rotation. 
The  final  parameter  needed  to  fully  define  the  polarization  of  a 
plane  wave  is  the  orientation  of  the  major  axis  with  respect  to 
the  coordinate  system.  This  parameter  is  called  the 
tilt  angle  (  r  ) . 

Figure  5  shows  the  major  and  minor  axes  of  an  elliptically  polar¬ 
ized  field.  The  tilt  angle  is  also  indicated.  For  consistency,  the 
value  of  the  tilt  angle  is  limited  to  the  range 

-90  <  r<+90  (6) 


v 


Figure  5.  Relationship  between  Major  and  Minor  Axes,  Tilt  Angle, 

and  Ellipticity  Angle 

Figure  5  defines  an  alternate  parameter  to  the  axial  ratio.  The 
ellipticity  angle  (  *  )  is  oimply 
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f  *  arctan  (1/r) 


(7) 


5.0  ORTHOGONAL  POLARIZATIONS 

Two  polarizations  are  said  to  be  orthogonal  if,  and  only  if,  their 

o 

axial  ratios  are  equal  and  their  tilt  angles  differ  by  90  degrees. 
Any  polarization  can  be  thought  of  as  the  superposition  of  two 
arbitrary,  elliptical  polarizations.  Equations  have  been  derived  to 
describe  polarized  fields  in  terms  of  the  following  orthogonal 
pairs: 

Left-hand/ Right-hand  Circular 

Horizontal  /  Vertical 

Left/Right  Slant  Linear  (±45  slant) 

General  Case  of  Orthogonal  Elliptical  Polarizations 
The  following  discussion  will  focus  on  the  horizontal  and  vertical 
pair. 

6.0  POWER  DENSITY  AND  EFFECTIVE  VALUES 
If  two  fields  are  orthogonal,  the  sum  of  the  powers  contained  in 
the  two  fields  is  equal  to  the  total  power  in  the  field. Given 
orthogonal  linear  polarizations  in  the  x  and  y  directions,  the 
total  power  density  in  the  wave  can  be  expressed  as 

Sw  *  S*  +  Sy  (8) 


where 


(9) 


and 


(10) 


(11) 


E  ,  E  ,  and  E  are  the  effective  values  of  the  fields.  For  linear 
a  y  w 

polarizations  the  effective  value  of  the  field  is  0.707  times  the 
peak  value.  For  circular  polarizations  the  effective  and  peak  field 
values  are  equal. 

7.0  POINCARE  SPHERE 

Poincar€  used  the  relationship  in  equation  (8)  to  construct  the 
unit  sphere  shown  in  Figure  6.  The  Poincari  Sphere11  is  a  very 
useful  graphical  aid  in  visualizing  polarization  relationships. 
Poincare  showed  that  the  polarization  of  a  wave  can  be  repre¬ 
sented  by  a  unique  point  on  the  surface  of  the  sphere.  The 
Poincarl  sphere  can  be  used  to  express  the  results  of  a  polar¬ 
ization  measurement,  rather  than  specifying  the  shape  and 
orientation  of  the  polarization  ellipse. 

Orthogonal  polarizations  are  located  at  opposite  poles  of  the 
sphere.  The  sphere  in  Figure  6  is  oriented  such  that  the  circular 
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Figure  6.  The  PolncarS  Sphere 


polarisations  are  located  at  the  North  and  South  poles.  The  family 
of  linear  polarisations  are  located  on  the  equator.  The  mutually 
orthogonal  axes  of  the  sphere  represent  the  (L,R),  (H,V),  and 
(i  slant)  polarisations. 

Defining  a  point  on  the  surface  of  the  sphere  can  be  done  in 

relation  to  any  of  three  spherical  coordinate  systems,  although 

one  will  be  more  appropriate  for  a  specific  measurement  technique 

than  others.  The  polar  angle  with  respect  to  each  axis  is 

determined  by  the  corresponding  polarisation  ratio.  The  polariza- 
12 

'ion  ratios  are  defined  as 

eR 

P  ■  tan  y  =  (circular  polarisation  ratio),  (12) 

EL 


Ev 

P.  *  tan  a  =  — l  (linear  polarization  ratio),  and  (13) 
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E_ 

PD  *  tan  &  *  (diagonal  polarisation  ratio)  (14) 

where  ER  and  EL  are  the  effective  values  of  the  right-hand  and 
left-hand  circular  components  of  the  electric  field.  Ey  and  E^  are 
the  effective  values  of  the  vertical  and  horixontal  components, 
and  E_  and  E+  are  the  effective  values  of  the  -45  and  +45 
slant-linear  components,  respectively.  The  polarisation  ratio  is  an 
indication  of  the  relative  contributions  of  two  orthogonal 
components  of  a  polarised  wave. 

The  longitude  of  a  point  on  the  sphere  about  a  particular  axis  is 
the  phase  angle  by  which  the  component  in  the  numerator  of  the 
polarization  ratio  leads  the  component  in  the  denominator.  The 
polarization  ratio  together  with  the  corresponding  phase  angle 
uniquely  specifies  a  polarisation. 

8.0  THE  POLARIZATION  BOX 

The  mutual  consistency  of  the  three  spherical  coordinate  systems 
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is  confirmed  by  the  polarisation  box,  which  is  shown  in  Figure 
7.  The  polarization  box  is  inscribed  inside  the  PoincarS  Sphere. 
It  is  a  graphical  aid  to  understanding  the  relationships  between 
the  three  basic  sets  of  polarisation  components. 

The  point  W  on  the  surface  of  the  Poincar$  Sphere  is  located  at 
the  corner  of  the  polarization  box  opposite  the  corner  located  at 
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Figure  7.  The  Polarization  Box 


the  origin.  The  three  axis  form  the  three  edges  of  the  box  that 
intersect  at  the  origin.  Specifying  the  dimensions  of  the  polar* 
isation  box  is  equivalent  to  specifying  the  polarization. 

9.0  POLARIZATION  MEASUREMENT  TECHNIQUES 

The  polarization  box  may  be  described  in  several  ways.  The 
method  of  description  usually  depends  on  the  actual  technique 
used  to  measure  the  polarization.  The  first  manner  in  which  the 
box  can  be  defined  is  by  specifying  its  three  dimensions  (i.e., 
length,  width ,  and  height).  These  dimensions  are  found  directly 
by  measuring  the  normalised,  effective  values  of  the  six  degen¬ 
erately  polarized  components.  In  other  words,  the  lengths  oi  the 
three  sides  are  given  by 
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(15) 


(16) 


E2  -  E2  (1?) 

«  “  •  ' 

To  measure  the  contributions  of  all  six  cardinal  polarisations 
would  require  an  antenna  array  with  all  six  receive  polarizations 
available . 

Another  way  to  define  the  polarization  box  involves  specifying  the 
three  longitudinal  angles  about  the  axes,  which  are  found  directly 
be  measuring  the  phase  angles  between  the  same  three  orthogonal 
pairs  used  above.  This  technique  requires  the  use  of  the  same 
antenna  array  with  multiple  receive  polarisations. 

Finally,  the  three  polar  angles  found  from  the  polarization  ratios 
(equations  12,  IS,  14)  can  be  specified  to  define  the  polarization 
box,  but  the  direct  measurement  of  the  polarisation  ratios  also 
requires  the  redundant  antenna  array  described  above.  Fortu¬ 
nately,  it  is  possible  to  synthesise  the  responses  of  all  six 
cardinal  polarizations  using  only  two  orthogonal  antenn  ts  by  pass¬ 
ing  the  received  signals  through  the  appropriate  phase  shifting 
networks.  Alternatively,  it  is  possible  to  calculate  the  responses 
of  an  orthogonal  pair  of  feeds  from  the  amplitude  and  phase 
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information  extracted  from  a  different  orthogonal  pair.  The  in¬ 
herent  relationship  between  the  basic  sets  of  polarisation  com¬ 
ponents  is  represented  by  the  polarisation  box  itself. 

To  find  the  lengths  of  the  remaining  sides  of  the  polarisation  box 
given  only  the  length  of  one  side  (which  can  be  measured  directly 
using  only  two  orthogonal  antennas)  and  the  longitudinal  angle 
about  that  side  at  which  the  point  W  exists  (which  is  equal  to  the 
phase  angle  between  the  responses  of  the  two  orthogonal  antennas 
used  to  measure  the  length  of  the  side)  requires  simple  geometry. 
In  fact,  any  of  the  parameters  of  the  polarisation  box  can  be 
calculated  from  this  information.  This  allows  simplification  of  the 
antenna  array  for  purposes  of  polarization  analysis,  but  increases 
the  requirements  of  the  data  processing  portion  of  the 
polarization  measurement  system. 


The  relationship  of  the  polarization  box  to  the  tilt  angle,  sense  of 
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rotation,  and  axial  ratio  is  straightforward.  The  tilt  angle  (  r) 
is  proportional  to  \'he  phase  (S')  between  the  left  and  right 
circular  components  of  the  field,  as  given  by  the  equation 


This  phase  angle  can  be  measured  directly  or  can  be  calculated 
using  the  geometry  of  the  polarization  box.  The  sense  of  rotation 
can  be  found  from  the  circular  polarization  ratio,  which  was 
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defined  in  equation  12.  If  the  circular  polarization  ratio  is  greater 
than  unity,  the  sense  of  rotation  is  clockwise  looking  in  the 
direction  of  propagation  (right-hand  sense).  If  p  is  less  than 
one,  the  sense  is  left-hand.  If  1  \e  orthogonal  elements  used  in 
the  antenna  array  are  not  circular,  the  circular  polarization  ratio 
can  be  calculated  from  the  relationships  described  above. 


The  axial  ratio  is  also  related  to  the  circular  polarization  ratio  by 
the  equation 


r 


P  +  J 
"  P-  1 


(19) 


10.0  ACCURACY  OF  POLARIZATION  MEASUREMENTS 
15 

Jensen  derived  a  general  expression  for  estimating  the  errors  in 
polarization  measurements  on  the  basis  of  the  relative  power  of  an 

— k. 

error  signal.  If  the  complex  vector  E  denotes  the  electric  field  of 

-A 

a  received  wave,  and  AE  the  limit  of  error  involved  with  the 
measurement  of  E,  then  E  +  AE  represents  the  actual  measure¬ 
ment  result. 

If  the  polarization  E  is  represented  by  a  point  D  on  the  surface 
of  the  Poincare  Sphere,  and  the  polarization  of  E  +  AE  by  the 
point  T,  then  the  angular  distance  between  D  and  T  reflects  the 
measurement  error.  Assuming  that  the  power  of  the  error  signal 
AE  reletive  to  the  power  of  E  is  known  to  be  less  than  unity  and 
given  by 
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(20) 


it  can  be  shown  that  the  upper  bound  value  of  V,  the  angular 
error,  is  given  by 


Vmax  *  2  arcBin  </?> 


(21) 


If  D  represents  the  polarization  of  the  wave,  then  the  measured 
polarization  will  be  represented  by  a  point  T  on  the  sphere  lo¬ 
cated  within  the  circular  area  Q  ,  which  is  characterized  by  a 


maximum  angular  deviation  from  D  equal  to  vmax>  as  shown  in 
Figure  8. 


Figure  8.  Uncertainty  of  Polarization  Measurements 


The  size  of  the  circular  area  Q  on  the  surface  of  the  Poincare 
Sphere  can  be  specified  for  a  particular  measurement  if  the  re¬ 
ceived  power  level,  and  the  errors  in  measured  parameters,  are 
known.  The  angular  error  in  a  polarization  measurement  can  then 
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be  translated  into  uncertainties  in  the  reported  tilt  angle, 
ellipticity  and  sense  using  the  geometric  relationships  of  the 
polarization  box. 

11.0  POLARIZATION  MEASUREMENT  SYSTEM 
The  polarization  measurement  system  discussed  in  the  following 
pages  is  based  on  Watkins-Johnson  Company's  extensive  experi¬ 
ence  in  the  design,  production,  and  support  of  microwave  receiv¬ 
ing  and  analysis  systems.  It  provides  a  powerful,  reliable  and 
accurate  system  in  a  compact,  cost-effective  package. 

The  block  diagram  of  the  system  is  presented  in  Figure  9.  The 
antenna  is  a  broadband,  dual-polarized  horn  housed  in  a  protec¬ 
tive  radome.  Two  orthogonal  ports  (horizontal  and  vertical 
polarizations)  collect  signals  in  the  2  to  18  GHz  frequency  range. 

The  signals  pass  through  matched  RF  Distribution  paths  to  a 
spread-spectrum  multiplexer,  which  is  essential  to  the  accuracy  of 
the  system.  The  multiplexing  technique  incorporated  in  this 
system  allows  signals  from  both  ports  to  be  received  and  analyzed 
simultaneously,  on  a  pulse-by-pulse  basis  without  introducing  the 
detrimental  effects  of  separate  tuner  channels.  By  encoding  the 
received  signals  at  the  antenna  output  and  downconverting 
through  a  single  tuner,  the  amplitude  and  phase  relationships  of 
the  signals  are  preserved,  allowing  a  more  accurate  measurement 
of  polarization  than  is  possible  with  conventional  systems. 
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Figure  9.  Signal  Polarization  Measurement  System  Block  Diagram 

The  RF  Distribution  components,  multiplexer,  and  wideband 
superheterodyne  tuner  are  contained  in  an  environmentally  pro¬ 
tected  enclosure.  The  proximity  of  the  multiplexer  and  tuner  to 
the  antenna  permits  the  simplification  of  the  RF  Distribution  by 
eliminating  the  need  for  preamplifiers  and  their  inherent  tracking 
problems . 

The  TN-122  wideband  tuner  downconverts  the  multiplexed  signal 
to  a  1  GHz  IF.  The  signal  travels  through  a  single  IF  cable  to 
the  demultiplexer,  thus  avoiding  the  problems  of  matching  a  pair 
of  long  IF  cables  and  maintaining  the  match  over  time.  At  the 
demultiplexer  the  individual  channels  are  recovered  for  processing 
by  the  Polarization  Analysis  Unit,  which  is  interfaced  to  the  host 
system  computer. 
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The  spread  spectrum  multiplexing  technique  was  adapted  from  the 

1 6 

WJ-1988  Monopulse  DF  System.  The  frequency  multiplexing  of 
two  channels  into  one  tuner  is  accomplished  by  applying  a  dif¬ 
ferent  BPSK  modulation  code  to  each  of  the  channels  prior  to 
combining.  The  mixers  in  the  multiplexer  unit  are  driven  with  a 
maximally  linear  sequence  similar  to  that  shown  in  Figure  10.  The 
sequence  simulates  pseudo-random  noise  but  is  repetitive  in 
nature  and  can  be  reconstructed  at  the  demultiplexer  to  recover 
the  information  in  the  individual  channels. 
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Figure  10.  Maximally  Linear  Sequence 


The  31  bit  sequence  used  in  this  system  is  generated  by  a  set  of 


shift  registers  with  the  appropriate  feedback.  The  circuitry  used 


to  generate  the  code  consists  of  high  speed  ECL  logic.  The  code 


rate  is  240  MHz. 


The  autocorrelation  property  of  maximally  linear  sequences  results 


in  noise  and  spurious  signals  unless  the  applied  code  in  the  de¬ 


multiplexer  is  in  time  synchronization  with  the  delayed  signal. 
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The  resultant  channel  to  channel  isolation  is  proportional  to  the 
length  of  the  code,  as  given  by  the  equation 

Isolation  =  20  log  (code  length)  (22) 

For  a  31  bit  code  word,  the  isolation  is  approximately  30  dB. 

The  spectrums  of  the  encoded  signals  are  shown  in  Figure  11. 
The  input  signals  are  displayed  in  the  top  photo.  After  passing 
through  the  multiplexer,  the  bandwidths  of  the  individual  chan¬ 
nels  are  many  times  those  of  the  original  signals.  The  spectral 
distribution  has  a  sin(x)/x  distribution,  with  a  main  lobe  width  of 
twice  the  code  rate,  or  480  MHz.  The  spectral  line  spacing  is 
equal  to  the  code  rate  divided  by  the  code  length.  It  is  interest¬ 
ing  to  note  that  the  balanced  nature  of  the  code  sequence  results 
in  the  suppression  of  the  carrier  frequency. 

At  the  multiplexer,  the  two  modulators  are  driven  with  time- 
shifted  versions  of  the  same  code  sequence.  As  long  as  the  time 
shift  is  greater  than  1  bit,  the  two  channels  are  orthogonal  and 
can  be  combined  without  unwanted  interferometry  effects.  The 
combined  signals  are  shown  in  the  center  photo  of  Figure  11. 

The  500  MHz  bandwidth  of  the  TN-122  tuner  is  required  to  sup¬ 
port  the  combined  signals.  The  requirement  for  a  full  code  word 
to  be  applied  to  a  pulse  for  full  isolation  sets  a  limit  on  the 
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Figure  11.  Bi-Phase  Encoding 
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effective  bandwidth  of  the  system.  The  minimum  pulse  width  that 
the  system  can  process  is  therefore  135  ns. 

The  composite  signal  is  downconverted  and  routed  to  the  demulti¬ 
plexer.  The  signal  is  split  and  routed  to  two  demodulators.  At 
this  point  the  identical  code  sequence  must  be  applied  to  the 
divided  signals  in  the  same  time  relationship  as  in  the  multiplexor 
to  recover  the  individual  channels,  as  shown  in  the  lower  photo¬ 
graphs  in  Figure  11. 

The  outputs  of  the  demultiplexer  are  two  IF  signals  with  the  same 
amplitude  and  phase  relationships  as  the  original  received  signals. 
The  preservation  of  these  relationships  is  critical  to  the  polar¬ 
ization  measurement  process,  as  was  shown  above.  The  observed 
performance  of  the  multiplexing/ demultiplexing  chain  is  quite 
good.  The  gain  tracking  is  on  the  order  of  2  dB,  and  the  phase 
tracking  is  3  degrees  RMS. 

The  Polarization  Analysis  Unit  (Figure  12)  converts  the  IF  signals 
to  log  video  for  amplitude  measurements  and  processes  the  IF  to 
recover  the  phase  information.  The  CPU  accepts  control  informa¬ 
tion  from  the  main  system  as  well  as  information  regarding  the 
amplitude  and  PRI  of  the  signal  of  interest.  The  main  system 
directs  the  TN-122  to  the  signal  frequency  and  directs  the  pol¬ 
arization  measurement  antennas  to  the  correct  azimuth.  The  CPU 
then  adapts  its  processing  threshold  accordingly  and  finds  the 
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pulse  train  of  interest  using  the  PRI  information  from  the  main 
system . 

After  locking  onto  the  correct  pulse  train,  the  Polarization 
Analysis  Unit  begins  measuring  relative  amplitude  and  phase  of 
each  incoming  pulse.  The  measured  data  is  analyzed  by  the 
numerical  data  processor.  The  polarization  data  is  returned  to  the 
main  system  and  presented  to  the  operator  via  the  front  panel 
display . 

The  Polarization  Analysis  Unit  is  capable  of  returning  polarization 
information  in  several  ways.  The  tilt  angle,  axial  ratio,  and  sense 


can  be  returned,  or  the  contributions  of  the  (H,V),  (RHC, 
LHC),  (+45,  -45)  components  can  be  calculated.  The  degree  of 
polarisation  can  also  he  calculated  for  the  primary  polarisation 
component . 

Calculations  of  the  polarization  parameters  are  based  on  the 
relationships  presented  in  the  preceding  paragraphs.  Having 
measured  the  amplitudes  of  the  horizontal  and  vertical  components 
and  the  phase  between  the  two  fields  allows  the  calculation  of  the 
dimensions  of  the  polarization  box.  The  dimensions  of  the  polari¬ 
zation  box  in  terms  of  the  measured  parameters  are  shown  in  Fig¬ 
ure  13. 


L 


Figure  13.  Tue  Polarisation  Bor.  in  Terms  of  Measured  Amplitudes  and 
Phase  of  Horizontal  and  Vertical  Components 


The  tilt  angle  is  related  to  the  polarisation  box  by  the  equation 


T 


6' 

T 


(23) 


Unfortunately,  the  phase  angle  5'  is  not  directly  measured  by 
this  system.  However,  the  angle  6'  can  be  calculated  from  the 
equation 


S/=  arctan 


f2EjjEyCos  8 
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(24) 


The  axial  ratio  was  defined  in  equation  19  in  relation  to  the 

circular  polarization  ratio.  Although  p  is  not  measured  in  this 
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system,  it  is  easily  calculated  using  the  equation 


(25) 


where 


X 


L 


2Pr 


i  +p: 


(26) 


The  sense  of  rotation  depends  on  whether  p  is  greater  than  or 
less  than  unity,  as  discussed  above. 


Once  the  dimensions  of  the  polsrtsstion  box  have  been  found,  it 
is  possible  to  identify  which  of  the  degenerate  polarisations  (H, 
V,  RHC,  LHC,  +4S,  or  -45)  the  received  signal  most  cloaely  re¬ 
sembles.  This  information  can  also  be  presented  as  output  data. 

The  degree  to  which  this  component  dominates  can  be  calculated. 

18 

The  degree  of  circular  polarisation  is  defined  by  the  equation 

CP  *  e£  -  E*  ■  2?hEv  sin  S  (27) 

where  CP  *  +1  for  a  purely  LHC  polarised  wave  and  -1  for  a 
purely  RHC  polarisation.  The  degree  of  linear  polarization  can  be 
defined  in  a  similar  manner  by  the  equation 

LP  =  E*  -  E*  (28) 

where  LP  =  +1  for  a  purely  horisontal  polarisation  and  -1  for  a 
purely  vertical  polarisation.  Finally,  the  degree  of  slant  polarisa¬ 
tion  is  given  by  the  equation 

SP  =  eI  -  E2  =  2EUE„  cos  5  (29) 

•  ▼  ••  •  n  *  v 


where  +1  indicates  a  purely  +45  slant  linear  polarisation  and  -1 
indicates  a  purely  -45  slant  linear  polarisation.  The  astute 
observer  will  realise  that  these  quantities  are  merely  the  lengths 
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of  the  appropriate  aides  of  the  polarisation  box.  The  curves  for 
the  degree  of  polarisation  can  be  plotted  as  shown  in  Figure  14. 
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Figure  14.  Degrees  of  Polarisation  versus  Polarisation  Ratio 
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12.0  SUMMARY 

Precision  polarisation  measurements  would  be  impossible  without 
the  excellent  amplitude  and  phase  tracking  characteristics  of  the 
orthogonal  antennas  employed  in  this  system.  Combining  this  per¬ 
formance  with  careful  matching  of  signal  paths  to  the  multiplexor 
ensures  the  accuracy  of  each  measurement. 

The  availability  of  precision  polarisation  measurement  equipment 
will  allow  systems  designers  to  exploit  all  signal  parameters  in  the 
pursuit  of  system  optimization.  However,  it  is  also  imperative  that 
systems  designers  become  aware  of  the  potential  that  polarimetry 
holds  for  system  enhancement.  Research  into  applications  for  pol¬ 
arization  measurement  equipment  must  continue  if  strides  are  to 
be  made  in  improving  the  capabilities  of  microwave  systems. 
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Abstract 


Over  the  past  decade  spherical  near-field  measurements  have  made  the  transition 
from  research  to  application.  Based  upon  spherical  modal  expansion  and  a  port-to-port 
transmission  equation,  the  spherical  near-field  technique  utilizes  sophisticated 
mathematical  processing  to  solve  the  practical  problem  of  obtaining  antenna  pattern 
measurements  in  a  confined  space.  Conventional  hardware  —  consisting  of  rotary  test 
positioners,  signal  sources  and  receivers  augmented  by  digital  control  electronics  and 
minicomputers  —  now  makes  spherical  near-field  measurements  practical  by  automation. 


Recent  developments  in  spherical  near-field  software  architecture  have  expanded 
the  basic  scanning  feature  to  include  dual-ported  scanning,  probe  pattern  and  probe 
polarization  correction,  three  antenna  polarization  measurements,  gain  comparison 
measurements  and  correction  for  thermal  drift  of  the  system.  J  .  ... 

In  this  presentation  I  review  the  theoretical  basis  for  the  method  and  describe  these 
more  advanced  features. 


The  text  of  this  article  originally  appeared  in  the  March  1985  issue  of  Microwave 
Systems  News  entitled  "Spherical  Near-Field  Antenna  Measurements  Improve  Through 
Expanded  Software  Features"  by  myself  and  co-authors  J.  R.  Jones,  C.  Green,  B.  Melson 
and  J.  Proctor. 
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INTRODUCTION 


Over  the  past  two  years  a  development  program  at  Sclentif ic-Atlanta  has 
supported  work,  on  an  expanded  software  package  for  performing  spherical  near¬ 
field  measurements  with  the  mini-computer  based  Model  2022A  Automatic  Antenna 
Analyzer.  In  this  article  we  review  the  features  of  this  expanded  software. 

The  original  spherical  near-field  software  package  for  the  Model  2022 
Antenna  Analyzer  was  designed  to  support  non-probe-corrected  single  channel 
near-field  measurements.  With  the  introduction  of  the  Model  1780  three 
channel  receiver  and  it’s  inclusion  in  the  Model  2022A  Antenna  Analyzer,  the 
spherical  near-field  package  was  expanded  to  provide  for  dual  ported  near- 
field  probes.  Now,  as  a  result  of  the  latest  software  improvements,  the 
following  features  are  included. 

e  Gain  Comparison  Measurement 

e  Three-Antenna  Polarization  Measurement 

e  Probe-Polarization-Corrected  Measurements 

e  Probe-Pattern-Corrected  Measurements 

e  Dual  Port  Probe  Correction  by  Channel  Balancing 

e  Thermal  Drift  Corrected  Measurements 

e  Polarization  Based  Output  Presentation 

Details  of  these  features  are  described  individually  later. 

An  easy  way  to  understand  spherical  near-field  measurements  are  is  to 
first  think  in  terms  of  conventional  far-field  measurements.  On  a  far-field 
range  a  source  antenna  illuminates  an  antenna  under  test  with  a  spherical  wave 
whose  phase  front  is  nearly  planar  and  whose  amplitude  distribution  is  nearly 
constant  across  the  aperture. 

On  a  far-field  range  the  test  antenna  is  placed  on  a  two-axis  positioner 
and  the  source  antenna  is  placed  on  a  single  polarization  axis.  The  two- 
dimensional  far-field  pattern  is  recorded  as  the  test  antenna  is  turned  in 
step-scan  fashion  through  its  two  axes.  The  angular  readout  values  correspond 
to  the  $  and  8  axes  of  the  standard  spherical  coordinate  system.  As  the  two- 
axis  rotary  motion  occurs  the  range  source  antenna  scans  over  an  imaginary 
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spherical  surface  fixed  to  the  test  antenna.  (The  radius  of  the  imaginary 
surface  corresponds  to  the  length  of  the  range).  Far-field  testing  is 
properly  thought  of  as  spherical  far-field  scanning. 

Spherical  near-field  testing  is  simply  the  far-field  configuration  with 
the  range  length  collapsed  to  a  small  distance,  usually  about  one  test  antenna 
diameter.  With  the  range  length  reduced,  measurements  can  be  made  indoors  in 
an  anechoic  chamber  where  laboratory  conditions  exist.  This  gives  near-field 
testing  a  great  practical  advantage  over  far-field  testing.  The  far-field 
pattern  is  obtained  by  a  mathematical  computation,  the  near-field  to  far-field 
transform. 

Spherical  near-field  testing  differs  from  far-field  testing  in  several 
important  aspects.  First,  phase  measurements  must  be  made;  amplitude  alone  is 
not  adequate.  Second,  solid-angle  scans  must  be  made;  principal  planes  alone 
measured  in  the  near-field  are  not  adequate  to  yield  the  correct  principal 
planes  in  the  far-field.  Third,  significant  signal  processing  calculations 
are  required  to  obtain  the  far-field.  These  features  require  the 
sophistication  of  automatic  data  acquisition  and  data  management  available  in 
the  Scientif ic-Atlanta  2022A  Antenna  Analyzer. 

THE  MATHEMATICAL  BASIS  OF  SPHERICAL  HEAR- FIELD  SCANNING 


The  mathematical  basis  of  spherical  near-field  scanning  is  the  reduced 
transmission  equation,  usually  written  in  the  following  notation: 


W(kr 


$  9  Y  ) 
o  o*o 
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smn 

opv 
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C  (kr  ) 

oiiv.syn  o 


D(n)($  9  x  ) 
pm  To  o*o 


smn 


where ; 


“(kr0i*0’a0,*0)  '  bi/ao 

is  the  relative  received  phasor  at  the  port  of  the  probe  antenna;  aQ  is  the 
input  phasor  at  the  port  of  the  antenna  under  test  and  b^  is  the  output 
phasor  from  the  probe. 
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The  output  of  the  probe  depends  parametrically  on  the  frequency  and  the  range 
length: 


k  ■  2ir/A 

r0  -  range  length 

The  output  of  the  probe  depends  on  the  position  of  the  probe  on  the  scanning 
sphere  and  the  orientation  of  the  probe  in  polarization;  these  physical 
quantities  are: 


azimuthal  angle  for  probe  position 

0  ■  polar  angle  for  probe  position 

XQ*  polarization  angle  for  probe  orientation 

(These  are  in  fact  the  Euler  angles  that  relate  the  pi.obe  ai.u  test  antenna 
coordinate  systems'  rotational  orientation). 

The  two  angles  and  6q  are  the  turning  angles  at  the  test  antenna  on  the 

two-axis  positioner.  To  within  an  offset  they  are  the  values  read  on  the 
synchro  or  encoder  display  of  the  test  positioner.  xo  is  the  polarization 
angle  of  the  range-source-probe  antenna;  to  within  an  offset  it  too  is  the 
value  of  the  synchro  readout.  The  right  hand  side  of  the  transmission 
equation  contains  scattering  matrix  elements  for  both  the  probe  antenna  and 
the  test  antenna. 

Tgan  *  transmitting  characteristics  of  the  test  antenna's  general 
scattering  matrix 

"  rece^v^-n8  characteristic  of  the  probe  antenna's  general 
scattering  matrix 
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These  quantities  can  also  be  thought  of  as  modal  expansion  coefficients.  They 
characterize  the  transfer  of  excitation  to  or  from  the  port  of  their 
respective  antenna  to  or  from  free  space  spherical  modes  labelled  as  follows: 

s,o  ■  TE  or  TM  label 

m, p  -  azimuthal  label 

n, v  *  principal  label 

The  indices  n  and  m  occur  in  the  associated  Legendre  polynomials  which  appear 
in  the  spherical  coordinate  solutions  to  the  Maxwell  equations.  The  TE,  TM 
index  is  required  to  distinguish  different  vector  solutions. 

The  right  hand  side  of  the  transmission  equation  also  contains  geometrical 
functions  which  relate  the  modal  coefficients  or  scattering  matrix  elements  in 
different  coordinate  systems. 

(0  ,6  ,x  )  *  Rotational  Transformation  Matrix 
pm  To  o **o' 

C  (kr  )  ■  Translational  Transformation  Matrix 

opv.sun'  o' 

These  functions  are  calculable  and  are  evaluated  in  the  process  of  the 
spherical  near-field  to  far-field  transform. 

If  one  knows  both  R'  for  the  probe  antenna  and  T _ for  the  test 

oyv  smn 

antenna  then  the  response  of  the  probe  b^  to  an  input  aQ  at  the  test 

antenna  can  be  computed  using  the  transmission  equation.  Unfortunately  this 

can  be  done  for  an  unknown  test  antenna  only  after  its  characteristic  has  been 

determined.  To  get  the  transmitting  characteristic  for  an  unknown  antenna, 

one  measures  b  '  as  a  function  of  sampled  values  of  6  ,  8  ,  and  by 
o  o  o  o 

Inverting  the  transmission  equation  gets  the  T  .  (It  is  assumed  that 

smn 

R'  for  the  probe  are  already  known  from  a  preliminary  measurement  step 

opv 

and  that  C  and  D  are  known  from  calculation.)  Having  obtained  ,  b^  may 

be  finally  calculated  for  a  far-field  value  of  r  with  a  dipole  probe. 
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Thus  the  necessary  ingredients  of  the  spherical  near-field  measurement 
process  are  the  capability  of  measuring  and  storing  the  probe  output  voltage 
bQ'  as  a  function  of  scanning  position  $0»90«X0  an^  of  reducing  the  data 
thus  obtained  to  yield  the  equivalent  far-field.  The  probe's  far-field 
pattern  can  be  measured  in  a  preliminary  step  on  the  same  range  on  which  the 
test  antenna  near-field  data  is  measured.  Its  modal  coefficients  or  receiving 
scattering  matrix  elements  are  obtained  from  its  measured  pattern  and  must  be 
stored  until  utilized  in  a  test  antenna  computation. 

The  features  of  automatically  acquiring  spherical  coordinate  scanning 
data,  of  storing  data,  of  performing  digital  calculations  with  the  data  and  of 
presenting  the  results  in  usable  hardcopy  formats  are  fundamental  to  the 
Scientif ic-Atlanta  2020  Antenna  Analyzer  design.  The  implementation  of  a 
software  system  for  spherical  ncar-field  scanning  was  straightforward  given 
the  pre-existing  architecture  of  the  far-field  measurement  software. 

The  main  addition  to  the  software  architecture  for  the  expanded  spherical 
near-field  system  has  been  the  files  for  storing  spherical  scattering  matrix 
transmitting  and  receiving  coefficients.  Files  have  also  been  added  to  store 
data  from  dual-ported  probes'  channel  balancing  procedures,  from  drift 
correction  procedures,  and  from  gain  transfer  measurement  procedures. 

SPHERICAL  WAR-FIELD  FEATURES  OF  THE  AHTEMHA  AHALYZER 

Spherical  Near-Field  Gain  Comparison  Measurements 

Gain  measurements  in  conjunction  with  near-field  scanning  can  be  carried 
ou  Ny  either  of  two  methods  -  the  gain  standard  substitution  method  or  the 
range  Insertion  loss  method.  The  gain  measurement  technique  Implemented  in 
the  S" '.entific-Atlanta  Antenna  Analyzer  is  the  gain  standard  substitution 
method. 

The  advantages  of  the  gain  standard  substitution  method  are  that  the  gain 
of  the  probe  antenna  does  not  have  to  be  known  and  that  the  range  does  not 
have  to  be  "short  circuited”  by  use  of  a  long  cable  to  determine  its  insertion 
loss.  Furthermore,  the  gain  standard  method  makes  use  of  existing  standard 
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gain  antennas  whose  gains  are  know  to  within  0.25  to  0.50  dB.  In  the  range 
inuertion  loss  method,  the  insertion  loss  of  the  near-field  range  (typically 
20  to  40  dfi)  oust  be  measured  to  within  the  accuracy  tolerated  by  the  error 
budget  (usually  0.25dB).  This  restriction  requires  additional  calibration  of 
a  detection  or  receiving  system.  The  gain  standard  substitution  method, 
although  it  usually  requires  the  scanning  operation  to  be  carried  on  the  gain 
standard  as  well  as  the  test  antenna,  is  readily  implemented  and  bears  a 
strong  resemblance  to  the  technique  used  on  far-field  ranges. 

Recall  that  for  the  far-field  case,  the  operating  equation  for  gain 
measurement  by  comparison  to  a  gain  standard  is: 

p  AUT 

GAUT  "  GSTD  p  STD 
r 

where  GAUT  and  GSTD  are  the  gains  of  the  unknown  test  antenna  and  to  gain 

AUT  STD 

standard,  respectively.  Here  and  Pr  are  the  power  levels 

received  under  polarization  matched  conditions  when  the  two  antennas  are 
illuminated  by  plane  waves  of  equal  power  density.  (Gain,  although  it  is 
defined  upon  the  assumption  that  an  antenna  is  transmitting,  is  usually 
measured  with  the  antenna  receiving.)  To  understand  how  gain  comparison 
measurements  are  done  with  spherical  near-field  scanning,  one  can  simply  take 
the  pure  far-field  case  and  rewrite  it. 

The  operating  equation  for  spherical  near-field  gain  comparison 
measurement,  written  in  terms  of  power  quantities  is  in  simplified  form: 

p  AUT-FF  p  AUT-NF  p  p  STI>-NF 

G  .  G  _£ _  _£ _  I  NT  _r _ 

AUT  STD  p  AUT-NF  p  STD-NF  STD-FF 

r  INT  r  r 

The  right  hand  side  consists  of  the  gain  of  the  standard  and  four,  rather  than 
one,  gain  comparison  factors.  The  second  and  third  of  the  four  factors  can  be 
combined  into  one  factor  if  desired,  which  is  the  comparison  between  the  power 
received  under  identical  near-field  illumination  by  the  test  antenna  and  the 
gain  standard;  at  times  it  is  convenient  to  use  an  intermediate  reference 
level  to  assist  in  the  comparison  measurement  when  one  passes  between  the  test 
antenna  and  gain  standard  level  comparison.  Hence  the  reason  for  two  terms. 
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The  first  and  fourth  comparison  terms  relate  the  power  levels  received 
under  near-field  illumination  to  those  received  under  far-field  Illumination; 
they  are  the  near-field  to  far-field  comparison  terms.  The  equation  used  for 
the  most  general  case  of  polarisation  correction  and  Impedance  mismatch 
correction  is  written  in  terms  of  field  quantities  rather  than  power 
quantities.  In  concept ,  it  is  identical  to  this. 

In  spherical  near-field  scanning,  one  augments  the  scanning  operation  by 
making  a  comparison  of  the  response  of  the  test  antenna  under  Illumination  by 
the  range-source-probe  antenna  to  an  intermediate  power  level  which  may  be  the 
AFC  signal  level  on  the  receiver.  In  the  2022A  Antenna  Analyzer  this 
comparison  is  stored  in  a  computer  disc  file.  Previously  this  intermediate 
level  has  been  calibrated  by  comparing  the  response  of  a  gain  standard  antenna 
under  Identical  illumination  conditions,  which  result  has  also  been  stored  in 
a  file. 

Scanning  operations  are  carried  out  on  both  the  gain  standard  antenna  and 
the  antenna  under  test.  The  Initial  points  of  the  scanning  operations  are 
identical  to  the  point  on  the  patterns  of  the  test  antenna  and  gain  standard 
where  the  comparison  measurements  were  made.  When  identical  transforms  are 
computed  to  both  antennas  the  points  of  Interest  in  the  far-field  radiation 
intensity  functions  can  be  compared  in  level  to  the  initial  near-field 
points.  These  comparisons  yield  the  near-field  to  far-field  gain  comparison 
terms  in  the  equation.  Thus  a  near-field  gain  measurement  consists  of  a  near¬ 
field  scan  plus  a  transfer  measurement  for  both  the  test  antenna  and  the  gain 
standard. 

In  the  2022A  Antenna  Analyzer  the  gain  of  an  unknown  antenna  is  computed 
by  a  program  that  goes  back  to  the  data  files  where  the  requisite  measurement 
data  are  stored.  The  value  of  the  gain  of  the  gain  standard,  the  near-field 
comparison  data,  the  near-field  scanning  data  and  the  far-field  radiation 
intensities  for  both  antennas  are  all  available  once  the  procedures  are 
carried  out.  The  gain  analysis  program  brings  the  data  together  for 
calculation  in  an  equation  that  closely  resembles  the  one  above. 
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Polarization  Measurements 


Because  electric  field  (and  magnetic  field)  it  a  vector  quantity,  no 
antenna  measurement  is  complete  without  a  determination  of  the  polarization 
characteristic  of  the  field.  Many  antenna  measurements  which  do  not  properly 
account  for  polarintion  suffer  a  lack  of  accuracy.  Polarization  is  an 
important  antenna  parameter. 

The  usual  method  of  polarization  measurement  Is  the  polarization  pattern 
method.  Two  antennas  in  their  mutual  far-field  regions  couple  via  the  two 
orthogonal  polarization  components  of  a  plane  wave  propagating  along  their 
line  of  sight.  A  polarization  pattern  is  made  by  rotating  a  linearly 
polarized  sampling  antenna  about  the  line  of  sight  in  the  field  of  an  unknown 
antenna.  The  peaks  and  nulla  of  the  familiar  pattern  determine  the  axial 
ratio  and  tilt  angle  of  the  field  of  the  unknown.  However,  no  sampling 
antenna  is  perfectly  linear  and  the  peaks  and  nulls  of  a  polarization  pattern 
are  in  truth  only  apparent  ones;  they  are  not  the  true  values  because  of  the 
finite  axial  ratio  of  the  sampling  antenna.  To  get  at  the  true  values  from 
the  measured  data,  one  needs  the  polarization  of  the  sampling  antenna  from  a 
previous  measured  result.  This  process  of  correcting  measured  polarization 
data  using  the  polarization  of  one  antenna  to  obtain  the  true  polarization  of 
the  other  is  called  a  two-antenna  polarization  transfer  measurement. 

In  many  cases  where  accurate  polarization  measurements  are  desired,  the 
polarization  of  the  sampling  antenna  needs  to  be  measured  as  well.  In  this 
case,  without  a  sampling  antenna  of  known  polarization,  a  three-antenna 
polarization  measurement  must  be  made. 

The  three  antenna  polarization  measurement  entails  taking  three  unknown 
antennas  and  making  three  polarization  patterns  pair  wise  among  the  set.  The 
three  polarization  patterns  may  then  be  analyzed  for  the  true  polarization 
characteristic  of  each  of  the  three  members  of  the  set.  (The  method  assumes 
that  none  of  the  three  antennas  is  circularly  polarized.)  The  three-antenna 
polarization  measurement  is  a  far-field  measurement  and  is  usually  carried  out 
on  smaller  antennas  (such  as  near-field  probe  antennas).  It  is  applied  one 
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frequency  at  a  time,  although  automatic  systems  make  multiple  frequency 
acquisition  straightforward. 

Once  the  polarization  of  a  probe  antenna  it  determined  it  nay  be  utilized 
to  perform  a  polarization-corrected  pattern  measurement,  in  which  the 
polarization  of  the  sampling  antenna  is  taken  into  account  and  corrected  for. 

The  expanded  spherical  near-field  software  package  in  the  2022A  supports 
both  two  and  three  antenna  polarisation  measurements.  An  acquisition  routine 
for  acquiring  polarisation  patterns  and  a  data  reduction  routine  are  Included 
an  features  of  the  2022A;  these  features  are  necessary  for  accurate  cross- 
polarization  measurements  with  spherical  near-fleld  scanning*  These  software 
routines  permit  users  to  calibrate  probe  antennas  for  later  use  in  near-fleld 
scanning. 

Piobe-Pattern-Corrected  Pattern  Measurements 


A  probe  that  is  used  to  measure  the  near-field  of  an  antenna  possesses  a 
response  which  weights  the  field  of  the  probe  over  Its  volume  to  produce  a  net 
signal.  The  manner  in  which  the  rest  antenna  near  field  is  weighted  can  be 
removed  from  affecting  the  measured  result  by  a  process  known  as  probe 
correction.  Probe  correction  presumes  that  one  knows  the  probe  antenna 
characteristic  in  advance  of  the  test  antenna  measurement  and  that  this 
characteristic  can  be  utilized  to  determine  the  far-ficld  of  the  test  antenna. 

Probe  correction  in  general  can  be  broken  into  three  categories  which 
describe  the  degree  of  correction  —  in  the  order  of  increasing 
sophistication: 

Oth  order  -  No  Probe  Correction 

1st  order  -  Probe  Polarization  Correction 

2nd  order  -  Probe  Pattern  Correction 
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In  the  ease  of  no  proba  correction,  tha  proba  la  asaunad  to  ba  an  idaal 
alanancal  dipole  which  raaponda  Co  Che  field  ac  a  poine.  Hare  no  preliminary 
knowledge  of  Che  proba  paccern  ia  assumed.  For  probea  which  are  highly 
linear,  this  la  ofeen  a  vary  acceptable  approach. 

In  cha  caae  of  probe  polarisation  correction,  cha  probe 'a  polarisation 
characteristic  ia  assumed  known  and  is  modelled  as  a  crossed  dipole,  as 
described  earlier.  Tha  probe's  polarisation  characteristic  may  have  been 
determined  either  by  a  two-antenna  or  a  three-antenna  polarisation 

measurement. 

In  the  case  of  proba  pattern  correction,  tha  complete  pattern  of  the 
probe  is  assumed  to  ba  known  by  a  preliminary  pattern  measurement.  Usually 
the  pattern  will  have  been  obtained  in  a  polarisation  corrected  measurement 
with  a  secondary  probe  whose  polarisation  haa  been  determined  one  step 
earlier. 

In  general,  whether  the  probe  correction  is  polarisation  correction  only 
or  a  pattern  correction,  the  probe  pattern  is  represented  by  an  expansion  of 
spherical  modes.  The  expansion  of  the  probe  pattern  in  modal  coefficients  is 
the  key  element  by  which  pattern  correction  operates.  For  polarisation 
correction  only,  the  expansion  is  truncated  abruptly,  keeping  only  the  dipole 
modes.  For  general  pattern  correction  only  enough  modes  are  utilised  to  give 
the  accuracy  desired;  too  many  modes  appear  to  give  difficulties.  The  modal 
coefficients  are  determined  by  applying  the  near-field  to  far-fleld  transform 
in  a  specialised  way,  reducing  the  measured  patterns  to  obtain  them.  Modal 
coefficients  for  the  crossed  dipole  model  are  obtained  by  reducing 
polarization  pattern  data. 

Probe  pattern  correction  (over  and  above  polarisation  correction)  has 
been  useful  only  for  small  probes  that  require  only  a  few  modes  to  represent 
their  fields.  It  is  Important  only  when  the  probe  is  very  close  to  the  test 
antenna  (within  a  few  wave  lengths). 

In  the  2022A  Analyzer,  probe  patterns  are  acquired  with  the  same  data 
acquisition  routine  that  is  used  to  acquire  test  antenna  data.  A  probe 
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pattern  can  ba  used  to  correct  teet  antenna  measurement*  once  the  pattern  is 
reduced  to  receiving  coefficient* .  This  ia  don*  by  the  same  spherical  near- 

i 

field  subroutine  that  yields  far-field  transformed  data*  j 

Probe-Polar last Ion-Corrected  Pattern  Measurements 


Once  the  polarisation  of  a  probe  antenna  is  known  that  information  may  be 
utilised  to  correct  any  pattern  measurements  made  with  that  probe.  This 
technique  is  termed  probe  polarisation  correction.  It  la  extremely  Important 
in  measurement  of  cross-polarised  patterns  in  order  to  achieve  accurate 
results.  Probe  polarisation  correction  in  spherical  near-field  scanning  has 
application  in  the  measurement  of  probe  patterns  and  in  the  measurement  of 
unknown  test  antennas. 

It  la  well  known  that  the  response  of  an  idealised  elemental  dipole 
antenna  Is  directly  proportional  to  the  linear  component  of  electric  field  at 
the  location  of  the  dipole.  Furthermore ,  one  knows  that  an  elemental  electric 
dipole  produces  a  perfectly  linearly  polarised  plane  wave  In  Its  far-fleld  at 
the  peak  of  Its  pattern. 

Analogously,  the  response  of  a  coupled  pair  of  crossed  dipoles  to  an 
electric  field  will  be  a  superposition  of  the  responses  of  the  crossed  dipole 
elements  to  two  orthogonal  linear  components  of  electric  field.  And,  in  the 
far-field  of  a  crossed  dipole  pair  that  is  transmitting,  a  superposition  of 
two  orthogonal  linearly  polarised  plane  waves  will  be  produced.  Intuitively, 
It  is  clear  that  the  polarisation  characteristic  of  a  probe  antenna  can  be 
related  to  the  equivalent  exitation  amplitudes  of  a  pair  of  crossed  dipoles. 
This  approach  is  extremely  useful  in  fitting  prube  polarisation  correction 
into  the  same  framework  as  probe  pattern  correction  to  be  discussed  later. 

Since  a  nonideal  probe  responds  to  both  linear  components  of  electric 
field  one  must  utilize  polarisation  correction  if  a  pattern  made  with  a 
nonideal  probe  can  be  modelled  as  a  pair  of  crossed  dipoles,  one  weakly 
coupled  in  comparison  to  the  other.  By  using  this  approach,  polarization 
correction  of  pattern  measurements,  either  near-field  or  far-fleld,  can  be 
looked  upon  as  a  special  case  of  probe  correction. 
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Probe  polarization  correction  is  especially  important  in  obtaining  the 
patterns  of  primary  probe  antennas  used  in  probe-pattern-correc.ted 
measurements.  Usually  the  primary  probe  and  a  secondary  probe  whose 
polarization  is  then  known  can  be  utilized  to  make  a  polarization-corrected 
measurement  of  the  pattern  of  the  primary  probe. 

Probe  polarization  correction  can  also  be  applied  to  spherical  near-field 
scanning  measurements  by  making  use  of  the  crossed  dipole  model  of  the  probe 
as  an  approximation  to  its  modal  expansion.  Very  often,  polarization 
correction  alone  is  adequate  to  give  perfectly  acceptable  far-field  patterns 
of  an  unknown  test  antenna. 

Dual  Port  Probe  Correction  by  Channel  Balancing 

The  strict  method  of  near-fleld  scanning  requires  two  scans  of  the 
measurement  surface  (one  for  each  polarization).  This  approach  has  always 
given  good  results  and  is  theoretically  correct.  However,  electromagnetic 
technology  has  provided  methods  of  fabricating  dual-ported  probe  horns  (now 
available  in  round  waveguide)  which  permit  both  polarizations  to  be  measured 
simultaneously  in  one  scan.  Probe  horns  have  been  developed  which  permit 
spherical  near-field  measurements  to  be  made  with  dual-ported  probe  horns  over 
half-octave  bandwidths. 

Each  port  of  an  ideal  dual-ported  probe  couples  to  only  one  linear 
component  of  the  field,  and  responds  equally  to  equal  field  strengths.  In 
reality,  however,  each  port  has  different  polarization  responses  and  field 
strength  responses. 

Not  only  is  the  probe  horn  slightly  imperfect  but  also  the  separate 
channels  of  the  receiving  system  are  imperfect.  Different  amplifier  gains  and 
different  phase  lengths  of  connecting  cable  require  matching  via  a  calibration 
procedure  called  channel  balancing. 
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Channel  balancing  is  carried  out  with  a  single  ported  utility  antenna  or 
a  secondary  probe  whose  polarization  is  known.  The  channel  balancing  antenna 
is  mounted  on  the  test  positioner  with  the  dual-ported  probe  mounted  on  the 
polarization  axis  of  the  range  tower. 

First  the  polarization  characteristic  of  each  port  of  the  dual-ported 
probe  can  be  found  from  a  two  antenna  polarization  transfer  measurement  if 
they  are  not  already  known.  Then,  since  the  polarizations  of  both  the  utility 
antenna  and  the  ports  of  the  dual-ported  probe  are  known,  the  polarization 
efficiency  factors  in  the  far-field  transmission  equation  can  be  evaluated. 

The  Imbalance  between  the  responses  of  the  ports  can  be  evaluated  by  co¬ 
polarizing  the  utility  antenna  to  each  port  separately  and  observing  the 
receiver’s  response  for  each  channel.  The  calculated  polarization  efficiency 
factors  and  the  observed  responses  provide  the  information  necessary  to 
correct  for  the  port  imbalance  and  the  receiving  system  Imbalance. 

In  the  2022A  Antenna  Analyzer's  new  spherical  near-field  software  this 
channel  balancing  procedure  with  polarization  correction  is  carried  out 
automatically  with  the  measured  data  recorded  and  stored  and  the  calibration 
performed  in  the  process  of  the  near-field  to  far-field  transform. 

Thermal  Drift  Correction 


Qjp-?,  of  the  primary  features  of  near-field  scanning  is  the  necessity  of 
making  phase  measurements.  Phase  measurements  on  an  antenna  range  are 
troubled  by  the  necessity  of  using  long  cable  runs,  as  compared  with  bench 
measurements.  RF  cable  tends  to  be  sensitive  to  changes  in  temperature  as 
small  as  several  degrees  Fahrenheit,  which  causes  a  drift  in  the  phase 
measurements.  This  drift  can  be  removed  from  the  measurement  if  it  is 
nnnitored  during  the  scanning  process. 

The  2022A  Antenna  Analyzer  corrects  for  thermally  Induced  drift  by  a 
method  called  the  return-to-peak  method.  In  this  technique  the  probe  is 
returned  periodically  during  scanning  to  the  near-field  peak.  The  same  point 
Is  remeasured  upon  each  return  and  the  variation  in  phase  and  amplitude  of  the 
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peak  value  is  used  to  produce  a  correction  factor  which  is  applied  to  each 
data  point  of  the  near-field  scan. 

The  accuracy  of  the  correction  is  determined  by  the  amount  of  drift  which 
occurs,  the  mathematical  form  of  the  true  drift,  the  frequency  of  the  returns 
to  the  near-field  peak  and  the  accuracy  with  which  the  peak  point  can  be 
measured. 

Typically  one  might  uet  up  an  antenna  measurement  to  record  data  on  a 
near-field  sphere  every  1  degree  in  and  9q  ,  the  polar  and  azimuthal 
spherical  angles  respectively.  Scanning  and  stepping  9q  gives  a  raster 
of  data  consisting  of  successive  scans  each  taking  several  minutes.  Provided 
that  the  thermal  drift  can  be  linearized  appropriately  over  the  time  required 
for  a  few  scans,  one  would  choose  to  return  to  the  peak  every  five  or  ten 
scans.  Experience  with  the  particular  range  often  yields  the  knowledge  of  how 
often  returns  ought  to  be  made. 

Polarization-Based  Output  Presentation 

It  has  long  been  a  frustating  fact  of  life  that  antenna  engineers  must 
deal  with  the  various  polarizations  that  the  electromagnetic  field  offers. 
Quite  often  an  engineer  will  wish  to  make  measurements  on  a  circularly 
polarized  antenna  on  a  range  equipped  with  linearly  polarized  source 
antennas.  Correction  for  imperfectly  polarized  source  antennas  has 
consistently  been  a  problem  for  those  interested  in  optimum  measurement 
accuracy.  Conversion  between  different  polarization  formats  has  also  been  a 
problem. 

The  2022A  Antenna  Analyzer  Spherical  Near-Field  software  permits  the  far- 
fields  to  be  presented  in  linear  and  circular  polarization  formats.  The 
linear  polarization  formats  can  be  referred  either  to  the  9  -  $  grid  system 
of  reference  directions  or  the  Ludwig  system  of  reference  directions. 
Furthermore  the  far-field  radiation  Intensity  function  can  be  presented  and 
plotted  just  as  if  it  had  been  measured  directly.  Also,  Poincare'  sphere 
parameters  can  be  plotted  versus  position  on  the  radiation  sphere. 
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The  output  presentation  features  of  the  spherical  near-field  software 
package  are  sufficiently  general  that  they  will  be  appropriate  for  any 
particular  antenna . 


CONCLUSIOa 

Whereas  the  original  2022  Antenna  Analyzer  Spherical  Near-Field  Software 
package  was  restricted  to  the  non-probe-corrected  case,  the  new  2021k  software 
is  designed  for  correction  of  all  recognizable  errors  for  which  solutions 
exist.  Among  the  errors  not  yet  correctable  are  (1)  errors  due  to 
unrestricted  probe  patterns  which  contain  azimuthal  index  values  other  than  ±1 
(2)  probe  position  errors  (3)  standing  waves  between  the  probe  and  tesr 
antenna  and  (4)  extraneous  reflection  signals.  It  is  likely  that  correction 
procedures  for  these  unsolved  problems  will  also  be  found  and  implemented  in 
future  generations  of  Antenna  Analyzers. 

In  the  meantime  the  2022A  Antenna  Analyzer  and  the  expanded  software 
option  -08A  for  spherical  near-field  measurements  enables  industrial  users  to 
make  use  of  the  advantages  of  near-field  testing. 
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SWEPT  FREQUENCY  TECHNIQUE  FOR  DISPERSION 
MEASUREMENT  OF  MICROSTRIP  LINES 
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1 .  SUMMARY 

i 

Mlcroitrlp  lines  used  In  microwave  Integrated  circuits  are 
dispersive.  Because  a  microstrip  line  Is  an  open  structure,  the 
dispersion  can  not  be  derived  with  pure  TEM,  TE,  or  TM  mode  anal¬ 
ysis.  Dispersion  analysis  has  commonly  been  done  using  a  spec¬ 
tral  domain  approach,  and  dispersion  measurement  has  been  made 
with  high  Q  microstrip  ring  resonators.  Since  the  dispersion  of 
a  microstrip  line  Is  fully  characterized  by  the  frequency  depend¬ 
ent  phase  velocity  of  the  line,  dispersion  measurement  of  micro¬ 
strip  lines  requires  the  measurement  of  the  line  wavelength  as  a 
function  of  frequency.  In  this  paper,  a  swept  frequency  tech¬ 
nique  for  dispersion  measurement  Is  described.  ^The  measurement 
was  made  using  an  automatic  network  analyzer  with  the  microstrip 
line  terminated  In  a  short  circuit.  Experimental  date  for  two 
microstrip  lines  on  a  10  and  30  mil.  Cuflon  substrates  were 
recorded  over  a  frequency  range  of  2  to  20  GHz.  Agreement  with 
theoretical  results  computed  by  the  spectral  domain  approach  Is 
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good.  Possible  sources  of  error  for  the  discrepancy  are 
discussed . 

2.  INTRODUCTION 

Microstrip  lines  used  In  microwave  Integrated  circuit  net¬ 
work  and  device  designs  are  very  attractive  because  of  low  cost 
and  easy  fabrication.  Unlike  ordinary  transmission  lines  of  TEM 
mode,  the  analysis  of  a  microstrip  line  Is  very  complex  and 
requires  the  hybrid  mode  expansion  of  both  the  TE  and  TM  modes  to 
account  for  the  fringing.  Furthermore,  a  microstrip  line  Is  more 
difficult  to  design  because  of  Its  dispersion  characteristics. 

The  dispersion  characteristics  of  microstrip  lines  have  been 

1  2 

Investigated  experimentally  and  theoretically.  •  Dispersion 

analysis  has  commonly  been  done  using  a  spectral  domain  approach, 

and  dispersion  measurements  have  been  made  with  a  high  Q  mlcro- 

2 

strip  ring  resonator.  This  paper  describes  a  swept  frequency 
technique  for  dispersion  measurements.  The  technique  Is  applied 
to  analyze  the  dispersion  characteristics  of  an  open  microstrip 
structure  as  shown  In  Fig.  2.  The  measured  results  are  compared 
to  the  analytical  results  computed  by  a  spectral  domain  method. 

3.  DISPERSION  MEASUREMENT 

Since  the  dispersion  of  an  open  microstrip  line  Is  fully 
characterized  by  Its  frequency  dependent  phase  velocity,  disper¬ 
sion  measurement  of  the  microstrip  line  requires  the  measurement 

2 

of  the  line  wavelength  as  a  function  of  frequency.  Troughton 


first  proposed  a  technique  for  measuring  wavelengths  and  dls- 
perlson  characteristics  using  a  microstrip  ring  resonator.  By 
capacltlvely  coupling  RF  power  Into  and  out  of  the  ring  resonator 
with  50  Q  transmission  test  probes,  the  line  wavelength  Is 
determined  by  measuring  the  transmitted  power  as  a  function  of 
frequency.  The  effective  dielectric  constant  at  each  resonance 
frequency,  f.  Is  obtained  from 

/  nc  V 

*eff  *  \2*rf/ 

where  c  Is  the  velocity  of  light  and  n  •  1,2,3 _  The  mean 

radius,  r,  of  the  ring  Is  approximated  from  Wheeler's  results4 
which  are  valid  for  frequencies  below  6  GHz.  Although  the  ring 
resonator  approach  eliminates  errors  In  determining  the  electri¬ 
cal  length  due  to  fringing,  the  curvature  of  the  ring  and  the 
coupling  probes  may  Introduce  errors  In  the  resonance  frequency 
measurements.  At  high  frequencies,  the  validity  of  Wheeler's 
theory,  for  approximating  the  mean  radius  r  Is  also  question¬ 
able.  In  fact,  large  discrepancies  between  calculated  and 

5 

measured  results  have  been  reported.  The  swept  frequency 
technique  shown  in  Fig.  2  Is  based  on  the  transmission  line 
theory.  The  dispersion  measurement  was  made  using  an  automatic 
network  analyzer  with  the  microstrip  line  terminated  In  a  short 
circuit.  Treating  the  microstrip  line  like  a  quarter-wavelength 
resonant  stub,  the  line  wavelength,  \,  may  be  calculated  by 
measuring  t  input  ance  as  a  function  of  frequency.  The 
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effective  dielectric  constant,  ceff»  at  each  resonance  can  be 
determined  from  the  following  relations: 

t  -  ^  (n  -  1,2,3...) 


where  t  Is  the  electrical  length  of  the  microstrip  line.  All 
of  the  measurements  were  made  on  1  In.  long  microstrip  lines  on  a 
Cufton  substrate  (*f  «  2.1).  The  line  width  Is  chosen  to  provide 
a  low-frequency  characteristic  Impedance  of  approximately  50  Q. 

A  SWA. 'HI C  flange  mount  connector  was  used  to  connect  the  micro¬ 
strip  line  to  the  automatic  network  analyzer  and  the  short  cir¬ 
cuit  termination. 

4.  NUMERICAL  ANALYSIS 

A  spectral  domain  method  as  outlined  In  Ref.  3  Is  utilized 
to  analyze  the  dlsperlson  of  an  open  microstrip  structure.  In 
brief,  the  formulation  of  the  problem  yields  a  hybrid  mode  solu¬ 
tion  of  linear  combinations  of  TE  and  TM  fields.  The  hybrid-mode 
fields  In  the  spectral  domain  are  given  by 

k?  -  k*  -Jk,z 

Ez1(«,y,z)  - — Jk  <*,y)e 

t-  k1  "  kz  -h  ‘Jkzz 

Sz1(«.y.g)  «  yk  1  4^(o,y)e 
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txi(«,y,z)  »  -  -j 


»v0vr  a^(a,y)  \  -Jkzz 


ay 


IVr  *♦?(»•*>  «h  1  -JM 
Hxi(«».y.z)  -  ^  ip  ^(a.y )je 

The  subscripts  1-1,2  designate  the  regions  1  (substrate)  and 

2(a1r).  The  Fourier  transforms  of  scalar  potentials  satisfy  the 

homogeneous  Helmholtz  equation 

2~  2~ 

y  o 


or 


where 


rl 


2~ 
a  a 

*2  2~ 
8?”  ' 


A  *r*V 


2  2.2.2 
y2  -  a  ♦  k2  -  k0 

.2  2 

k0  "  “  Vo 

By  matching  tangential  components  (x  and  z)  of  the  fields  at 
the  air-dielectric  Interface,  and  satisfying  the  radiation  con¬ 
dition  and  appropriate  boundary  conditions  at  the  conducting 
strip,  two  coupled  equations  of  the  unknown  transformed  currents 
(3*  ,T  )  and  the  transformed  longitudinal  fields  W  are 
obtained. 
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+  *  ^*2^a,d^ 

621(ci,kz>  \(a)  ♦  622(a,k2)Jz(a)  -  Hz2(a,d) 

where  the  expressions  for  the  6reen's  functions  G^Ca.k^), 

(a,k2)  etc.  are  given  In  Ref.  3.  To  find  the  unknown  propa¬ 
gation  constant  kz,  Galerkln's  method  In  the  spectral  domain  Is 
applied  to  reduce  the  coupled  equations  to  a  matrix  equation.  By 
expanding  the  transformed  current  3*x  and  7*  In  terms  of  some 
known  basis  functions  J*xn  and  3*zn,  l.e.. 


3x<*>  ■  cnV»> 
n»i 


*»<•>  «  £  dn?7n<o) 


A  matrix  equation  of  the  unknown  constants  cn  and  dn 


Is  obtained. 


where 


£  Kil,1)cn  +  E  *  0  m  «  1,2, ...N 

ntl  ^"n  n  n-i  mn  n 

£  *in’1>Cn  *  Z  Kin*2)dn  *  0  -  1,2,...H 

nTl  ffin  n  n«l  nm  n 


C1'  ■/’3™<*,S ll<-kz>  3*n<*>d« 

C*1  -  [  3»<«>  V-V  V»>d“ 
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The  unknown  propagation  constant,  k^,  Is  calculated  numerically 
by  finding  the  root  to  the  determinant  of  the  matrix  equation 

[K]  .  0 

The  root  to  the  determinant  Is  computed  by  Iteration  technique. 
Romberg's  numerical  Integration  algorithm  has  been  used  to  com¬ 
pute  K  The  numerical  result  Is  strongly  dependent  on  the 
mn 

choice  of  the  basis  functions.  Different  bases  function  for  3^ 
and  have  been  reported  In  the  open  literature.  For  the 
work  here,  the  following  basis  functions  have  been  chosen. 


■  0  otherwise 


V*>  ■  JX0  S^n  0.8W  ,x|  -  0,8  2 

■  Jxo  dx  °'8  *  «  1*1  i  ? 

5.  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES 
Figure  3  shows  the  measured  and  computed  values  of  the 
effective  dielectric  constant,  «  ^,  as  *  function  of  fre¬ 
quency  for  a  microstrip  line  of  0.89  mm  line  width  on  a  10  mil 
Cuflon  substrate  («r  «  2.1).  The  measured  curve  exhibits  the 


usual  oscillatory  characteristics  of  the  coaxial  microstrip 

9 

transition.  If  the  substrate  thickness  Is  Increased,  the 
transition  effect  Is  enhanced  as  clearly  Indicated  In  Fig.  4  for 
the  30  mil  substrate.  The  transition  effect,  If  not  removed, 
will  obscure  the  small  deviations  In  Impedance  and  phase  velocity 
due  to  the  dispersion  of  the  microstrip  line.  To  remove  the 
oscillations  experimentally  requires  a  broadband,  well  matched 
coaxlal-mlcrostrlp  transition.  However,  by  fitting  the  measured 
data  to  an  Interpolating  polynomial  to  smooth  out  the  transition 
characteristics,  the  measured  curves  show  good  agreements  with 
the  computed  curves.  The  maximum  deviation  Is  approximately 
5  percent  over  a  frequency  range  of  2  to  20  GHz.  The  discrepancy 
between  the  measured  and  the  calculated  values  can  be  traced  to 
different  sources  of  error.  Since  the  physical  line  length  Is 
used  to  determine  the  line  wavelength,  the  effective  dielectric 
constant,  *eff  computed  from  the  measured  line  wavelength  gen¬ 
erally  has  a  higher  value.  The  fringing  at  the  end  of  the  line 
tends  to  increase  the  electrical  length  of  the  line,  and  thus, 
lowers  the  value  of  the  measured  By  Including  the  fring¬ 

ing  effect,  a  closer  agreement  between  the  measured  and  the  com¬ 
puted  values  can  be  obtained.  Also,  It  has  been  observed  experi¬ 
mentally  that  the  type  of  short  termination  has  a  profound  effect 
on  *eff.  By  terminating  the  microstrip  line  with  a  SHA  coaxial 
short  and  by  wire  bonding  the  microstrip  line  to  Its  ground 
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plane,  the  measurements  can  be  changed  by  a  few  percent.  For 
the  latter  case,  the  Impact  of  fringing  on  the  electrical  length 
Is  significant.  Better  agreement  between  the  analytical  and 
experimental  results  are  also  possible  If  a  more  accurate  basis 
function  for  the  unknown  current  or  a  better  analytical  model  for 
the  microstrip  line  Is  used.  Finally,  the  Impedance  variation 
along  a  short-circuit  microstrip  line  can  be  obtained  by  measur¬ 
ing  the  Input  Impedance  as  a  function  of  frequency  or  electrical 
lengths.  Figures  5  and  6  show  the  variation  of  the  normalized 
resistance  and  reactance  with  frequency.  The  resistance  curve 
shows  an  Impulse  type  of  response  at  resonant  frequencies,  while 
the  reactance  curve  demonstrates  similar  characteristics  of  a 
lossless  short-circuited  transmissions  line  of  TEH  modes.  It  Is 
Important  to  point  out  that  because  of  the  dispersion  of  the 
microstrip  line,  the  resonant  frequencies  for  the  higher  order 
modes  do  not  always  occur  at  exactly  the  multiple  of  the  funda¬ 
mental  frequency.  As  a  consequence,  the  phase  velocity  and  the 
effective  dielectric  constant  are  frequency-dependent. 

7.  CONCLUSION 

A  swept  frequency  technique  for  dispersion  measurement  of 
microstrip  lines  has  been  described.  The  technique  Is  easy  to 
apply,  and  yet,  yields  fairly  accurate  results  for  both  thin  and 
thick  substrate  over  a  frequency  range  of  2  to  20  GHz.  From  this 
technique,  one  Is  able  to  obtain  the  Impedance  variation  along  a 
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dispersive  line,  and  a  better  physical  understanding  of  the  dis¬ 
persive  characteristics  of  an  open  microstrip  line. 
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/  ABSTRACT 


For  certain  applications  the  antenna  designer  may  require  a 

knowledge  of  the  small  scale,  or  localized,  reflectivity 

properties  of  microwave  absorbers.  Suppliers  of  microwave 

absorbers  generally  quote  the  average  absorption  achieved  over 

fairly  large  panels  (610  mm  x  610  mm)  as  determined  by  one  of 

several  free-space  measurement  techniques.  By  the  nature  of 

their  construction,  small  samples  of  planar  absorber  panels 

cannot  generally  be  evaluated  directly  in  waveguide,  particularly 

y 

at  X-eand  and  above,  -"yc _ _ 


This  paper  describes  a  swept-frequency ,  one-horn  interferometer 
where  the  absorber  panel  is  placed  directly  against  the  aperture 
of  a  well  matched  horn.  By  careful  selection  of  waveguide 
components  and  by  aperture  matching  an  exponential  horn,  residual 
reflections  within  the  system  can  be  reduced  to  acceptably  low 
levels  (typically  better  than  -35  dB).  The  measurement  procedure 
has  the  advantage  that  fairly  small  areas  (about  the  size  of  the 
horn  aperture)  of  planar  absorber  panels  can  be  evaluated  without 
having  to  cut  the  panels  to  fit  pieces  into  a  sample  holder. 


Reflectivity  measurements  on  a  variety  of  planar  absorber  panels 
show  that  multilayer  Jaumann  absorbers  have  the  greatest 
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inhomogenelty  (reflectivity  can  vary  by  10  or  even  20  dB  over  a 
few  centimeters).  Some  multilayer,  graded  dielectric  absorbers 
show  excellent  homogeneity  down  to  the  ~30  dB  level.  Honeycomb 
absorbers  have  excellent  homogeneity  within  a  panel  but 
occasionally  show  variations  from  panel  to  panel  when  the  ' 

honeycomb  is  excessively  loaded.  The  technique  conveniently 
determines  the  center  frequency  of  resonant  absorbers.  ‘ 

1 .  INTRODUCTION 

For  certain  applications  (e.g.  the  microwave  absorber  in  the 

backing  cavity  of  a  planar  spiral  antenna),  the  reflectivity  of 

the  absorber  over  quite  small  areas  (50  mm  x  50  mm)  can  influence 

antenna  performance1.  Small  samples  of  absorber  are  sometimes 

measured  in  waveguide  but  this  technique  i3  not  satisfactory 

because  the  physical  construction  of  most  planar  panels  make  them 

poorly  suited  for  insertion  into  small  waveguides.  Relatively  i 

small  samples  of  absorber  can  be  measured  by  a  horn 

interferometer  where  the  sample  is  placed  a  meter  or  so  from  the 

p  i 

horn  .  In  this  case  a  tuner  is  usually  used  to  reduce  the  | 

residual  reflection  from  the  horn  at  each  frequency  of  interest. 

] 

Evaluation  of  absorbers  over  even  modest  frequency  bands  can  be  j 

an  extremely  time-consuming  and  tedious  process. 

* 

This  paper  describes  a  one-horn  interf erometer  where  the  residual 
reflection  from  the  horn  is  reduced  to  a  low  level  (typically 
better  than  -35  dB)  by  eliminating  reflections  from  the  throat  of 
the  horn  and  significantly  reducing  reflections  from  the  horn 
aperture.  This  is  achieved  by  using  an  aperture-matched 
exponential  hTn.  Because  the  horn  can  be  matched  over  the  full 
bandwidth  of  the  standard  waveguide  banas ,  swept-frequency 
measurements  can  be  taken  over  1.5:1  frequency  bands. 

-0 16- 

j 

_____  _ _  _  j 


3 


In  order  to  examine  relatively  small  samples,  the  planar  absorber 
panel  is  placed  directly  against  the  aperture  of  the  horn  and  its 
reflectivity  compared  to  that  of  a  metallic  shorting  plate  which 
sets  the  zero  dB  reference  level.  Comparative  measurements  over 
various  parts  of  the  absorber  panel  can  be  made  rapidly  without 
having  to  cut  samples  from  the  panel.  Measured  results  are 
presented  for  a  variety  of.  absorbers  to  illustrate  some  of  the 
observed  inhomogeneities  in  reflectivity. 

2.  MEASUREMENT  PROCEDURE 

The  measurement  system  consists  of  an  RF  swept  source,  some  means 
of  isolating  the  source,  10  dB  directional  couplers  with  a 
minimum  of  40  dB  directivity,  aperture-matched  exponential  horn, 
RF  step  attenuator  and  network  analyser  (see  Figure  1).  The 
system  is  configured  out  of  components  generally  available  in  a 
microwave  laboratory  and  lends  itself  to  rapid,  real  time 
evaluation  of  absorber  panels  over  comparatively  broad  (1.5:1) 
frequency  ranges. 

Figure  2  shows  the  aperture-matched  horn.  The  aperture  matching 
is  achieved  by  attaching  a  20  mm  diameter  rolled  edge  to  all  four 
walls  of  a  conventional  exponential  horn  (e.g.  model  640  of  the 
NARDA  Microwave  Corporation).  Theoretical  and  experimental 
aspects  of  aperture  matching  are  described  by  Burnside  and 
Chuang^.  For  reasonable  aperture  matching  the  cylinder  diameter 
should  be  about  a  wavelength  at  the  lowest  operating  frequency. 
Experience  has  shown  that  significantly  improved  VSWR  performance 
can  be  achieved  by  using  even  quite  small  cylinders  (down  to  a 
quarter  wavelength  in  diameter)1*.  By  choosing  20  mm  cylinders 
(0.53  A  at  8  GHz),  the  aperture  of  the  horn  at  the  tangential 
points  to  the  cylinders  is  not  increased  too  much. 
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The  reflection  coefficient  of  the  conventional  NARDA  horn  is 
shown  in  Figure  3  together  with  that  of  the  aperture-matched 
horn.  Without  aperture  matching  the  reflection  coefficient  is 
about  -31  dB  ( VSWR  1.06:1)  and  with  aperture  matching  better  than 
-37  dB  (VSWR  1.03:1)  over  the  8 . h  to  12  GH2  band.  The  dips  at 
about  9.1,  10.3  and  H.3  GHz  result  from  the  finite  directivity 
(about  HM  dB,  in  this  case)  of  the  couplers  used.  This  effect  is 
analogous  to  that  observed  when  measuring  antenna  patterns  in  the 
presence  of  extraneous  reflections  and  can  be  analyzed  in  the 
5 

same  way  . 

The  residual  reflectivity  (-37  dB)  of  the  aperture-matched  horn 
will  introduce  measurement  uncertainty  depending  on  the  amplitude 
and  phase  of  the  reflected  signal  from  the  absorber.  For 

example,  a  true  reflected  level  of  -25  dB  from  the  sample  could 
lie  anywhere  between  -23  and  -27.5  dB  depending  on  whether  it 
adds  in  or  out  of  phase  with  the  residual  reflection.  At  a  true 
-30  dB  level,  the  residual  reflectivity  could  cause  this  to  be 
between  -26.8  and  -35  dB.  Measured  reflectivities  below  the 
-25  dB  level  should  therefore  be  treated  with  caution. 

The  measurement  procedure  is  as  follows:  place  the  metallic  plate 
against  the  horn  aperture  to  obtain  a  zero  dB  reference  level  (a 
better  baseline  is  achieved  by  placing  the  plate  at  the  waveguide 
flange  between  the  coupler  and  the  horn),  check  RF  linearity  and 
dynamic  range  with  the  RF  attenuator  (60  dB  linear  dynamic  range 
was  achieved),  place  the  absorber  panel  against  the  horn  aperture 
with  the  metallic  plate  behind  it  ensuring  that  the  sample  is 
firmly  pressed  against  both  the  aperture  and  the  plate.  The 
reflection  coefficient  with  the  sample  in  place  is  the 
reflectivity  of  the  absorber  panel.  The  aperture  of  the  horn  is 
not  free-space  nor  is  the  incident  wave  plane,  however  for 
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comparative  measurements  this  is  of  no  consequence.  Measurements 
on  known  samples  show  that  the  test  results  for  homogeneous 
panels  do  not  deviate  significantly  from  the  free-space  results. 

3.  MEASURED  RESULTS 

A  wide  variety  of  planar  absorber  panels  were  evaluated  using  the 
procedure  of  Section  2  to  ascertain  how  homogeneous  common 
absorber  panels  are  when  sampled  over  small  areas  (about  60  mm  x 
80  mm  in  X  band)  within  the  panels  (typically  610  mm  x  610  mm  or 
305  mm  x  305  mm).  The  absorbers  tested  fall  into  three  broad 
classes:  (a)  multilayer  dielectric  absorbers,  (b)  carbon  coated 
honeycomb  panels  and  (c)  thin,  resonant  absorbers  (single  layer 
or  double  layer).  Chapter  9  of  reference  2  gives  a  detailed 
discussion  of  absorber  types. 

3. 1  Multilayer  dielectric  absorbers 

Broadband  absorbers  can  be  made  by  using  multiple  layers  of 
resistive  sheets  and  low  dielectric  constant  spacers.  Several  of 
these  Jaumann  absorbers  were  evaluated,  the  absorber  panels 
differing  in  the  number  of  resistive  sheets  and  spacer 
thicknesses  but  all  specified  to  have  -17  dB  reflectivity  in  the 
8  to  18  GHz  frequency  band. 

Figure  4  shows  the  swept-frequenoy  results  for  a  three  resistive 
layer  Jaumann  absorber.  The  two  plots  were  taken  at  positions 
near  the  center  of  a  610  mm  x  610  mm  panel  with  the  panel 
displaced  about  100  mm  between  the  measurements.  There  are 
localized  inhomogeneities  where  the  reflectivity  differs  by 
between  5  and  10  dB  to  the  -25  dB  level.  The  largest  difference 
in  reflectivity  observed  for  any  of  the  Jaumann  absorbers  tested 
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is  shown  in  Figure  5  for  a  six  sheet  absorber  measured  over  12  to 
18  GHz.  A  difference  in  reflectivity  of  about  20  dB  was  observed 
when  the  horn  was  moved  only  50  mm.  The  absorber  was  dissected 
to  see  whether  the  inhomogeneity  could  be  ascribed  to  mechanical 
defects  -  the  spacers  were  all  of  correct  thickness  and  the 
resistive  layers  were  not  damaged.  The  difference  in  performance 
can  only  be  ascribed  to  the  fact  that  the  lossy  layers  are  not 
sufficiently  homogeneous  and  isotropic.  When  the  absorber  panels 
are  to  be  used  to  cover  large  surfaces,  the  localized 
inhonogeneities  in  reflectivity  are  of  no  consequence  since  these 
variations  will  be  averaged  out. 

A  multilayer,  carbon  loaded,  low  density  foam  absorber  panel  was 
also  evaluated  (see  Figure  6).  The  three  curves  show  the  extent 
of  variations  over  the  610  mm  x  610  mm  panel.  Down  to  the  -30  dB 
level,  the  variation  is  only  about  4  dB.  which  is  typical  of  the 
graded  foam  absorbers  tested. 

3 . 2  Honeycomb  absorbers 

Honeycomb  absorbers  are  made  by  coating  the  hexagonal  open  cell 
honeycomb  with  a  resistive  layer.  These  absorbers  are  often  used 
to  line  antenna  cavities1  and  are  thus  required  to  have  good 
homogeneity  and,  the  reflectivity  should  be  the  same  (or  nearly 
so)  from  panel  to  panel  of  the  same  type.  The  honeycomb  itself 
is  not  homogeneous,  there  being  a  series  of  ribbons  where 
adjacent  hexagonal  cells  meet.  Figures  7  and  8  show  the 
reflectivity  from  8  to  1 2  GHz  and  1 2  to  18  GHz,  respectively.  In 
each  case  there  are  two  curves  -  one  for  the  electric  field 
parallel  to  the  ribbon  and  the  other  for  the  electric  field 
perpendicular  to  the  ribbon  direction.  The  two  plots  track  each 
other  1  airly  closely,  the  perpendicular  plot  being  typically  1  or 
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2  dB  better.  Figures  7  and  8  were  measured  using  different 
one-horn  interferometers.  There  is  excellent  tracking  of  the 
reflectivities  at  the  -20  dB  level  at  the  cross-over  frequency  of 
12  GHz. 

In  general,  excellent  homogeneity  in  reflectivity  was  found 
within  a  sheet  and  several  sheets  of  the  same  type  tracked  each 
other  to  within  a  dB  from  8  to  12  GHz.  One  particular  honeycomb 
ansorber  deviated  from  the  above.  The  absorber  sheets  were  quite 
thin  (12.7  mm)  and  had  been  heavily  loaded  tc  achieve  -20  dB  at 
X  band.  Figure  9  shows  the  difference  in  reflectivity  between 
two  sheets  of  the  same  type,  within  each  sheet  the  homogeneity 
was  good  (about  1  dB) .  This  particular  honeycomb  was  odd  in  the 
sense  that  the  absorption  decreases  with  increasing  frequency. 

3 •  3  Thin,  resonant  absorbers 

The  s'.'ept -frequency ,  one  horn-interferometer  can  be  used  to 
determine  the  resonant  frequency  of  thin  narrowband  absorbers. 
Figure  10  shows  the  resonant  frequencies  of  three  narrowband 
absorbers.  These  measurements  were  taken  at  the  center  of  each 
305  mm  x  305  mm  panel.  Within  each  panel  the  resonant  frequency 
shifted  by  about  500  MHz  as  the  horn  was  moved  around.  The 
manufacturer  does  not  ouote  the  center  frequency  for  the  absorber 
but  simply  states  that  the  absorption  will  be  better  than  -17  dB 
at  9.3  GHz.  Since  these  absorbers  generally  have  high  dielectric 
constants,  srn^il  changes  in  thickness  will  cause  a  shift  in  the 
resonant  frequency. 

Figure  11  shows  an  interesting  resonant  absorber  where  two 
resonances  are  achieved  in  a  ^.7  mm  thick  panel  to  increase  the 
bandwidth  over  that  of  the  conventional  resonant  absoroer  of 
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Figure  10.  The  solid  curve  shows  the  typical  reflectivity  over 
about  90%  of  the  sheet  and  the  dasheu  curve  shows  the 
manufacturer's  typical  data  which  was  taken  using  an  arch  far 
field  range.  The  solid  curve  of  Figure  11  shows  good  agreement 
with  the  arch  measurements.  In  part,  this  agreement  results 
because  the  sheet  tested  with  the  interferometer  is  quite 
homogeneous . 

4.  CONCLUSIONS 

Details  of  a  svept-freouency  interferometer  for  making  rapid, 
real  time  assessments  of  localized  inhomogeneities  in  planar 
absorber  panels  have  been  presented.  An  aperture-matched, 
exponential  horn  was  used  to  reduce  residual  reflections  in  the 
system  to  about  -37  dB.  This  residual  reflection  is  adequate  for 
making  comparative  measurements  on  planar  absorber  panels  whose 
reflectivities  usually  fall  in  the  -15  to  -25  dB  range.  Jaumann 
absorbers  have  the  worst  homogeneity,  vhile  the  honeycomb 
absorbers  generally  have  excellent  homogeneity  within  a  sheet  and 
from  sheet  to  sheet.  The  center  frequencies  of  resonant 
absorbers  were  also  measured  using  the  test  set-up.  With 
suitable  directional  couplers  and  aperture-matched  exponential 
horns,  the  technique  can  easily  be  applied  in  the  standard 
waveguide  bands  from  2  to  MO  GHz. 
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FIGURE  1  :  Swept-frequeney t  one-horn  interferometer. 
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REFLECTIVITY  (dB) 


FREQUENCY  (GHz) 


FIGURE  3  :  Measured  reflection  coefficient  of  two  exponential 
horns. 


FIGURE  :  Reflectivity  of  5  layer  Jaumann  absorber  measured  at 
points  separated  by  100  mm. 


REFLECTIVITY  (dB) 
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FIGURE 


FIGURE 


5  :  Reflectivity  of  6  layer  Jaumann  absorber  measured  at 
points  separated  by  50  mm. 


j6  :  Reflectivity  of  carbon  loaded  loam  absorber  showing 
variation  over  610  mm  x  610  mm  sheet. 
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FIGURE  7  :  Reflectivity  of  honeycomb  absorber  with  E  field 
parallel  and  perpendicular  to  ribbon. 


FIGURE  8  :  Reflectivity  of  honeycomb  absorber  with  E  field 
parallel  and  perpendicular  to  ribbon. 
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REFLECTIVITY  (dB) 


FREQUENCY  (GHz) 


FIGURE  9  :  Reflectivity  of  two  honeycomb  absorber  sheets  of  the 
same  type. 


FREQUENCY  (GHz) 

FIGURE  10  :  Measured  resonance  frequencies  of  three  narrowband 


absorber  panels. 


FREQUENCY  (GHz) 


FIGURE  1 1  :  Comparison  between  measured  data  and  the  typical  data 
supplied  by  the  manufacturer  for  a  doubly  resonant 
absorber  panel. 
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Abstract 

~V 

The  method  for  maximization  of  the  directivity  of 
Yagi-Uda  arrays  of  shaped  dipoles  as  outlined  by  Liang  and 
Cheng*''"  uas  employed  to  optimize  a  three-element  array  of 
shaped  dipoles,  each  1,5  wavelengths  in  length  and  0.01 
wavelengths  in  radius.  Liang  and  Cheng’s  optimal  three 
element  array  uas  confirmed  to  yield  a  maximum  directivity 

of  11.8  dB  using  the  moment  method  uith  seco;.d-order 

.  A- 

expansions  as  formulated  by  Harrington*". 

s' 

C 

Two  simplex  minimization  routines  were  alternately 
used  to  independently  determine  optimal  element  shapes  and 
eMminate  saddle  points  for  a  t hree  -e  1  ement  array  of  the 
same  dipole  lengths  and  radii,  1.5  and  0.01  wavelengths, 
respectively.  Using  this  procedure,  a  shaped-dipole 
Yagi-Uda  array  different  from  Liang  and  Cheng’s  example 
was  found  to  yield  a  maximum  directivity  of  12.6  dB. 

X 


-631- 


This  paper  will  also  present  results  for  optimal 
shapes  of  Yagi-Uda  arrays  of  more  than  three  elements. 

1 .Introduction 

The  Yagi-Uda  array  has  important  applications  in 
the  VHF/UHF  frequency  range  due  to  its  simplicity.  A 
conventional  Yagi-Uda  array  consists  of  a  row  of  parallel 
straight  cylindrical  conductors,  of  which  only  one  is 
driven  by  a  source  and  all  others  are  parasitic.  Among 
the  parasitic  elements,  one  serves  as  a  reflector  and  the 
rest  are  directors.  A  conventional  Yagi-Uda  array  of 
three  elements  is  pictured  in  Figure  1. 

Methods  for  increasing  the  directive  gain  of  this 
type  of  antenna  abound.  The  methods  include  changing 
element  spacing,  element  radius,  element  length,  boom 
length  and  number  of  directors.  Element  lengths  are 
usually  on  the  order  of  one-half  wavelength. 

Landst orf er 1  reasoned  that  for  dipole  elements 


longer  than  one  wavelength,  shaped  dipoles  could  yield 


a 


higher  directive  gain  than  linear  ones  of  the  same  length, 
assuming  a  sinusoidal  current  distribution  that  is 
independent  of  the  shape  or  radius  of  the  dipole.  By 
experimentation,  he  found  that  a  three-element  Yagi-Uda 
array  of  shaped  elements,  each  1.5  wavelength  long,  could 
be  adjusted  to  yield  a  directive  gain  of  11.5  dB.  Cheng 
and  Liang*  proposed  a  numerical  optimization  procedu  e  for 
Yagi-Uda  arrays.  Each  dipole  was  described  by  three 
parameters;  dipole  length  and  radius  were  held  constant  at 
1.5  and  0.01  wavelengths,  respectively,  for  the  array. 
Their  procedure  used  the  method  of  moments  and  point 
matching  to  approximate  the  array  as  a  multiport  network, 
and  a  simplex  routine  to  maximize  the  directivity  of  the 
array.  Liang  and  Cheng  achieved  a  maximum  gain  of  11.8dB, 
which  could  only  be  implemented  with  a  conventional 
Yagi-Uda  array  of  nine  or  ten  linear  elements. 

The  increase  in  the  directivity  of  a  shaped  three 
element  array  indicates  the  possibility  of  increasing  the 
gain  or  Yagi-Uda  arrays  of  more  than  three  elements. 
Thus,  a  study  was  undertaken  to  inverstigate  shaped  arrays 
of  more  than  three  elements.  The  invertiga  ion  is 
numerical  in  nature  and  paralleled  that  of  Liang  and 
Cheng.  The  moment  method  numerical  technique  employed  in 


-63  3- 


this  work  uas  verified  by  reproducing  the  results 
published  by  Liang  and  Cheng*  When  an  optimization 
procedure  uas  instituted  uith  their  results  as  a  first 
guess,  the  procedure  did  not  converge;  the  intermediate 
results  indicated  that  there  existed  other  optimal  shapes* 
Based  on  these  intermediate  results  it  uas  decided  to 
examine  a  three  e’ement  array  as  uell* 

In  this  paper,  the  results  of  numerical  analysis 
of  multi-element  Yagi-Uda  arrays  are  presented.  In  the 
foliouing  the  numerical  moment  method  formulation  of  an 
axially  symmetric  array  is  presented.  Also,  a  brief 
description  nf  the  optimization  procedures  used  in  this 
study  is  provided.  The  results  of  the  optimization  study 
of  optimal  shapes,  element  spacings  for  three  element 
array  are  described  in  the  last  section.  In  that  section, 
the  results  of  directivity  and  impedance  for  arrays  of 
more  than  three  elements  are  also  included.  Finally,  some 
general  conclusions  are  offered* 

2.  Formulation 

As  described  by  Liang  and  Cheng*  ,  each  element  is 
described  by  the  foliouing  equation; 
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y  <  x)  =  At  C  1  -  1/(1  +  BixMD  +  C,  (1) 

where  i  =  1,  2,  3»  N,  and  N  is  the  total  number  of 

elements*  The  parameters  At  »  Bt  ,  and  Ct  are  to  be 

optimized.  The  relation  between  the  total  arc  length  L  of 

each  element  and  half  of  its  projected  length  on  the 
x-axis  is  as  follows: 

r  Xo  i 

L  =  2  vr"T"+"T3y73xTr'l  dx  (2) 

O  J 

where  x,t  is  the  x-coordinate  of  the  end-point  of  the 
antenna  element  number  i. 

Each  element  is  divided  into  segments  such  that 
one  segment,  is  centered  on  the  y-axis.  Each  segment  is 
then  numbered  and  considered  to  be  a  port  in  a  multi-port 
linear  network.  Following  Harrington2,  the  self  and 
mutual  impedances  of  each  segment  are  given  by 

2<m,n)  =  jwu  'dl.  +  (1/juE)  [  <y(n*  ,m*  )  -  t<n'  ,m*  ) 

+  »Mn*  ,m'  )  +  t(n'  ,m'  )  3  (3) 


where , 
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4>(m,  n) 


( l/'8Tt«) 


Exp  ( -  j  k  R.  n  ) 


(4) 


-a 


Ra  n 


dz’ 


j  =  square  root  of  -1, 

•j  -  frequency  in  radians* 

U  =  permeability  of  free  space* 
£  =  permitivity  of  free  space, 
dl.  =  length  of  segment  m, 

2  a  =  I  Jl*  I 


r  */-Tpr-”Tz"-2r7rl  ,  m  *  n 
L  //-Car-+~TzrTrl  ,  m  =  n 

uire  radius 


The  second  order  approximations  of  the  scalar 
function  of  equation  (4)  from  Harrington  are  used,  and  the 
geometry  for  determination  of  this  function  is  described 
i n  F i gure  2 . 


The  directive  gain  for  one  polarization  is  as 


foil ous : 


g ( 9 , 0) 


uher  e , 


ok2 

4n 


I  C  V,  (8,H>,  ]  [  Y  ]  [  Vs  3  1 2 
Pin 


(5) 


n  =  intrinsic  impedance  of  free  space, 
k  =  2n/A  =  the  wave  numoer, 

Vr  (8,0)  —  Er  a  *  din  ,  ffl  =  1 , 2 , 3  ,  •  •  •  ,N 
Era  (8,0)  =  Ur  eXP  (  ~  j  k  r  T»  )  . 
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Note  that  V,  and  E.  are  induced  voltages  and  electric 
fields,  respectively,  induced  on  a  unit  dipole  receiving 

antenna  by  the  antenna  segment  m. 

The  input  pouer  of  a  source  feeding  the  nth 
segment  of  an  antenna  is 

Pi«=Re{IV*l2Y(n,n)}  (6) 

uhere,  is  the  source  voltage,  Y(n,n)  is  the  input 

admittance,  n  is  the  feed  port.  The  current  is  given  by 
the  product  of  the  admittance  matrix  and  source  voltage 
vectors. 


CI]  =  LY3CVs  1  (7) 

C  V8  1  =  C  0,  0 . V„  ,  0 . 0  V  . 

The  admittance  matrix  is  the  inverse  of  the  impedance 
matrix: 


CYJ  =  HZ]'  ‘ 


(8) 


The  above  mentioned  expressions  were  modeled  on  a 
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computer  to  determine  the  directive  gain,  current 
distribution,  bandwidth,  and  input  impedance  for  an 
axially  symmetric  array.  This  computer  code  was  verified 
by  checking  with  the  known  results  for:  a  three  element 
Yagi-Uda  array  of  linear  elements;  a  single  linear  dipole 
and  an  array  of  three  V-shaped  dipoles;  and  most 
importantly  the  results  of  Liang  and  Cheng  for  three 
shaped  element  array. 

The  numerical  task  of  solving  the  above  moment 
method  equations  is  well  known  and  will  not  be  discussed. 
Suffice  it  to  say  that  the  code  was  carefully  verified.  In 
the  next  section,  a  brief  discussion  of  the  optimization 
procedure  used  in  this  work  is  provided. 

3.  Procedure 


The  optimal  analysis  of  the  Yagi-Uda  arrays 
carried  out  here  utilizes  tuo  different  simplex  routines. 
A  simplex  routine  is  essentially  a  root-finding  routine 
for  a  function  of  N  variables.  It  will  determine  the  set 
of  N  values  of  a  function  of  N  variables  which  will  yield 
a  minimum  value  of  the  function.  It  scans  the  "space", 
the  range  of  possible  values  for  all  the  function's 
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variables,  evaluating  the  function  for  many  different  sets 
of  variable  values.  When  a  minimum  is  suspected,  it 
contracts  the  space  it  is  searching  and  conducts  a  more 
detailed  search.  This  process  is  reiterated  until  a  local 
minimum  is  found  or  the  search  proves  futile  (in  which 
case  it  scans  a  much  larger  space  and  the  process  starts 
again)  . 

For  a  function  of  N  variables,  if  the  function's 
inverse  is  found  to  have  a  minimum  for  a  set  of  N  values 
of  those  vari  ables,  then  this  same  set  of  variables  is 
that  which  will  yield  a  maximum  for  the  function. 
Usually,  the  minimum  (or  maximum  for  a  function’s  inverse) 
is  only  the  minimum  within  a  certain  range  of  values, 
termed  a  "neighborhood",  of  the  function’s  variables,  or 
constraints.  This  is  termed  a  "local"  minimum.  A  minimum 
value  of  the  function  for  the  entire  range  of  possible 
values  of  the  function  is  termed  a  "global"  minimum. 
Different  sets  of  values  which  yield  different  local 
minima  that  are  not  the  global  minimum  are  termed  "saddle 
points"  (See  Priemer* ). 

Any  effort  to  reduce  the  number  of  variables  will 
greatly  increase  the  speed  of  the  optimization  procedure. 
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Liang  and  Cheng  did  just  that  by  setting  one  of  the 
independent  variables  to  equal  zero  (fixing  the  antenna 
configuration  at  a  reference  point)  and  thus  reducing  the 
number  of  independent  variables  from  nine  to  eight. 

Liang  and  Cheng  used  this  procedure  to  determine 
the  three  sets  of  parameters  necessary  for  maximum 
directive  gain  on  a  three-element,  shaped  Yagi-Uda  array 
of  element  length  1.5  wavelengths  and  radius  0.01 
wavelength.  The  procedure  required  approximately  40 
minutes  of  CPU  time  on  a  DEC  10  computer,  and  yielded  a 
set  of  parameters  describing  an  antenna  uhose  directivity 
is  11.8  dB,  which  could  only  be  implemented  with  a 
conventional  Yagi-Uda  array  of  nine  or  ten  dipoles1. 

Usually,  for  functions  that  require  a  very  short 
amount  of  CPU  time,  one  can  allow  the  simplex  routine  to 
attempt  to  converge  on  the  global  minimum.  However,  for 
long  routines  such  as  moment  method  code  used  here, 
(approximately  two  seconds  in  the  simplest  case:  three 
elements),  one  must  print  intermediat e  results  in  case  it 
fails  to  converge . 

In  the  numerical  analysis  carried  out  hare,  two 
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simplex  routines  were  alternately  used  as  follows: 

"MINMUM"3  t  This  routine  has  constraints  defined  by  how 
much  the  function  variables*  values  are  varied  from 
initial  first  guesses. 

"ZXMUD":  This  routine  is  part  of  the  IflSL  subroutine 
library  generally  available  on  major  mainframe 
computers.  It  uses  ranges  of  the  function  variables* 
values  as  constraints. 

A  first  run  of  the  simplex  using  MINMUM  with  Liang 
and  Cheng’s  parameters  as  a  first  guess  failed  to  converge 
on  a  local  minimum  within  two  hours  of  CPU  time  on  an  IBM 
3081K64  machine.  This  is  due  to  the  fact  that  even  with 
Ct  fixed  at  zero,  the  dimension  of  the  space  was  eight, 
and  the  function  required  approximately  two  seconds  of  CPU 
time  for  each  evaluation.  However,  intermediate  results 
of  the  optimization  were  printed,  and  they  indicated  a 
global  minimum  existed  other  than  the  one  obtaineJ  by 
Liang  and  Cheng.  Numerous  saddle  points  (local  minima) 
were  also  indicated. 

ZXMUD  was  used  to  search  for  a  local  minimum 
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within  the  neighborhoods  of  each  saddle  point*  This  was  a 
very  costly  process,  requiring  approximately  two  hours  of 
CPU  time*  ZXMUD  rise  failed  to  converge  for  the  same 
reason  as  MINMUM,  but  it  was  found  that  the  neighborhood 
of  a  saddle  point  could  be  sufficiently  investigated  for  a 
local  minimum  within  only  two  minutes  of  CPU  time. 

Therefore,  for  the  case  of  the  three-element 
shaped  Yagi-Uda  array,  using  Liang  and  Cheng's  parameters 
as  a  first  guess,  program  MINIMUM  was  used  to  scan  a  large 
neighborhood  around  the  first  guess,  and  ZXMUD  was  used  to 
investigate  the  indicated  saddle  points  within  much 
smaller  neighborhoods. 

Unfortunately,  for  the  cases  of  shaped  Yagi-Uda 
antennas  of  elements  numbering  more  than  three,  the 
benefit  of  such  an  excellent  first  guess  was  not 
available.  Therefore,  it  was  necessary  to  try  and  find 
one  by  evaluating  the  function  at  widely  spaced  variable 
values  throughout  the  entire  range  of  the  functions.  For 
this  purpose,  to  reduce  the  space,  the  A,  variables  were 
fixed  at  a  constant  value  for  the  entire  array;  the  Ct 
values  were  fixed  such  that  the  interelement  spacing  was 
constant,  and  only  the  B,  parameters  were  varied.  Thus, 
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for-  a  four-element  array,  the  search  space  was  reduced 
from  12  independent  variables  to  four.  Similarly,  for  a 
five-element  array,  the  space  uas  reduced  from  15  to 
five.  Also,  to  further  reduce  the  CPU  time  required  for 
each  evaluation  of  the  function  for  a  set  of  parameter 
values,  each  element  uas  divided  into  15  segments  instead 
of  the  21  used  by  Liang  and  Cheng.  Houever,  once  a 
fruitful  minima  uas  achieved,  21  segments  uere  used. 

4.  Results  and  Discussion 


Using  the  above  procedure  each  parameter  uas 
optimized  independently,  fixing  only  C2  at  0.000.  This 
procedure  required  approximately  40  minutes  of  CPU  time  on 
the  IBM  3081K64.  The  following  parameters  uere  found  to 
yield  a  maximum  directive  gain  of  12.6  dB  and  an  input 
impedance  of  6.99  +  j32.5  ohms. 

Table  1 

Parameters  of  a  Three  Element  Array  with  a  Directivity  of 

12.6  dB. 


E 1  ement 

A, 

Bi 

Ct 

Reflector 

0.378 

15.05 

-0.162 

Active 

0.393 

32.11 

0.000 

Director 

0.362 

77.69 

0.151 
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The  parameters  found  in  Table  1  were  found 
separately  by  optimizing  each  parameter  except  C2  . 

More  interesting  results  were  found  uhen  the  space 
uas  reduced  by  setting  the  value  of  At  parameters  equal  to 
0.370  and  the  inter-element  spacing,  D,  equal  to  0.160. 
Then  Ct  is  given  by  (i-l)D.  The  Bi  parameters  were  then 
separately  optimized.  The  space  for  a  three  element 
antenna  is  now  reduced  to  three,  the  number  of  B, 
parameters.  The  results  of  this  procedure  are  found  in 
Table  2. 

Table  2 

Parameters  of  a  Three  Element  Array  with  a  Directivity  of 

13.2  dB. 


Element 

Bi 

Ref 1 ector 

15.00 

Act i ve 

31.00 

Director 

50.00 

For  antennas  of  N  elements,  the  above  procedure 
uas  duplicated.  All  parameters  except  C2  uere  optimized. 

Then,  the  8,  parameters  were  separately  optimized  together 
with  a  constant  At  equal  to  0.370  and  interelement  spacing 
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of  D  =  0.160.  The  optimized  Bi  parameters  are  presented 
in  Table  3  for  N  equals  three  to  nine. 


Table  3 


Optimized  Parameter  Bi  and  Directive  Gain  for  Multi-element  Arrays 


N 

Bi 

B* 

b3 

B« 

Bj 

B» 

B? 

B. 

B, 

GAIN 

3 

15 

31 

50 

13.2 

4 

22 

30 

58 

73 

13.4 

5 

15 

30 

70 

117 

130 

13.9 

6 

12 

24 

75 

113 

160 

72 

14.1 

7 

13 

30 

60 

■S3 

140 

160 

■Ml 

14.2 

8 

12 

38 

65 

91 

116 

144 

166 

181 

14.4 

mm 

12 

37 

90 

115 

140 

166 

182 

211 

14.5 

Note  that  the  values  of  the  Bi  parameters  are  to 
tuo  significant  digits.  For  N  equal  to  five  through  nine, 
A?  is  set  to  0.370  and  interelement  spacing  to  0.160. 


For  the  procedure  used  in  this  study,  it  uas  found 
that  the  best  results  for  directive  gain  is  achieved  uhen 
the  Bi  parameters  are  optimized  aione,  uhen  compared  to 
optimization  over  the  entire  space.  That  is,  A2  is  fixed 
to  0.370  and  the  interelement  spacing  is  set  at  0.160.  It 
should  be  noted  that  this  has  been  confirmed  for  antennas 
of  three  and  four  elements  only. 
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By  reducing  the  space  over  uhich  the  optimization 
was  carried  out,  and  by  reducing  the  number  of  dimensions, 
the  number  of  possible  saddle  points  on  uhich  the 
optimization  will  attempt  to  converge  is  reduced.  This 
also  increases  the  chances  of  having  any  local  minimum  to 
be  a  global  minimum;  thus  increasing  the  efficiency  of 
the  optimization. 


--  mvv  »» 


In  this  paper,  the  results  of  a  numerical  analysis 
of  shaped  dipole  Yagi-Uda  arrays  are  presented.  This 
analysis  uas  performed  using  the  method  of  moment  and  the 
work  presented  by  Liang  and  Cheng*  Their  optimization  gas 
confirmed  using  this  method*  The  optimization  of  the 
array  presented  by  Liang  ar.d  Cheng  uas  refined  for  greater 
directive  gain* 

This  uas  accomplished  by  alternately  using  tuo 
optimization  routines.  It  uas  found  that  by  reducing  the 
space  ever  uhich  the  optimization  uas  performed  greater 
computing  efficiency  uas  achieved  and  most  importantly, 
directive  gain  uas  increased.  Using  Liang 
and  Cheng’s  results  as  a  first  yuess  and  optimizing  ovsr 
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the  uhole  space,  a  directive  gain  of  12.6  dB  uas  obtained. 
But,  by  fixing  the  inter-element  spacing  and  one  of  the 
parameters  for  each  element,  an  increase  in  the  directive 
gain  to  13.2  dB  uas  achieved.  The  procedure  uas  then 
extended  to  an  array  of  four  elements  and  a  directive  gain 
of  13.4  dB  uas  found.  The  results  for  up  to  nine  elements 
have  also  been  presented.  Specifically,  the  appropriate 
parameters  and  the  directive  gain  have  been  tabulated. 

An  interesting  observation  that  shou’d  be  made 
that  analogous  to  the  case  of  a  linear  Yagi-Uda  array,  an 
increase  in  the  number  of  director  beyond  the  first  tuo 
does  not  increase  the  directivity  significantly.  This 
observation  is  based  on  our  restrictive  optimization  of 
more  arrays  uith  more  than  three  directors  (five  elements) 
and  a  full  optimization  may  produce  significant  change  in 
the  directive  gain.  Houever,  the  computational  resources 
needed  to  carry  out  such  an  exhaustive  analysis  are 
significant  and  no  attempt  uas  made  to  carry  out  such  an 
analysis  here. 


It 

i  s 

planned  to 

cont i nue 

this  study . 

In 

part i cul ar , 

arrays 

uith  four 

and 

more 

e 1 ement  s  will 

be 

expl ored 

further 

.  Also 

an 

i nvest i gat i on  into 

the 

bandwidth  would  be  of 

considerable 

interest . 
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ABSTRACT 


This  paper  describes  an  attempt  to  predict  the  received  voltages 
on  a  VHF  direction  finding  (DF)  array  using  the  geometrical  theory 
of  diffraction  and  the  principle  of  reciprocity.  Once  the 
received  voltages  are  known  bearing  errors  can  be  predicted  from 
the  phase  information. 


The  predicted  receive  mode  radiation  patterns  are  compared  with  a 
simple  model  of  a  C130  Hercules  transport  aircraft.  Further  com¬ 
parison  is  made  with  models  representing  a  low-wing  jetliner  with 
underwing  engines.  The  latter  aircraft  models  use  either  a  blunt 
capped  cylinder  or  a  sectional  cylinder  with  a  hemispherical  cap 
and  a  conic  section.  The  comparison  between  predicted  and 
measured  results  are  reasonable. 

It  is  concluded  that  a  simplified  model  with  a  blunt  cylindrical 
model  yields  acceptable  results.  Furthermore  for  underbelly 
mounted  DF  arrays  the  smallest  bearing  errors  are  obtained  for 
high-wing  aircraft  followed  by  low-wing  aircraft  with  tailmounted 
engines.  Low-wing  aircraft  with  underwing  enginepods  will  yield 
the  worst  results. 
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1. 


INTRODUCTION 


It  is  known  from  previous  modelling  of  Direction  Finding 
(DF)  antennas  on  aircraft  and  ships  1  that  reradiation  from 
structures  in  the  vicinity  of  the  antennas  can  be  a  major  source 
of  errors.  In  this  paper  attempts  to  predict  the  magnitude  of 
structural  effects  using  standard  computer  techniques  and  to 
verify  these  predictions  by  model  measurements  are  described. 

The  initial  decision  was  to  use  a  reasonably  simple  physical  model 
of  the  Lockheed  C130  (Type  L100-30)  Hercules  transport  aircraft. 
This  aircraft  is  used  extensively  in  a  variety  of  roles.  In  South 
Africa  it  has  been  used  in  several  longrange  search  and  rescue 
operations  by  SAFAIR  Freighters.  Some  of  these  missions  are 
listed  in  Table  1.  In  view  of  the  fact  that  frequencies  at 
approximately  121.6,  156.8  and  243.0  MHz  are  designated  for  call¬ 
ing  and  distress  use,  it  seemed  to  the  authors  appropriate  that 
this  aircraft  would  be  a  suitable  case  study  for  an  analysis  of 
structural  scattering  effects  on  a  simple  VHF  DF  array.  The  com¬ 
plexity  of  mounting  appropriate  antennas  on  the  underbelly  of  the 
real  aircraft  in  view  of  the  low  ground  clearance  was  not  con¬ 
sidered  important  for  the  purpose  of  this  study. 

In  view  of  the  somewhat  disappointing  results  obtained  during  com¬ 
parison  of  numeric  and  physical  modelling  it  was  decided  to  extend 
the  measurements  to  include  a  model  of  what  could  be  considered  a 
fairly  typical  low-wing  commercial  jetliner  with  two  underwirig 
engine  pods.  In  order  to  reduce  the  complexity  of  the  computer 
model  it  was  also  decided  to  investigate  the  use  of  a  blunt  capped 
cylinder  for  the  fuselage  as  an  alternate  to  a  more  accurate  scale 
model . 

2.  THEORETICAL  CONSIDERATIONS 
2.1.  Introduction 

The  principles  of  the  method  for  computing  the  effect  of  an 
aircraft  structure  on  the  performance  of  an  airborne  DF  array  are 
described  here. 

Since  a  DF  system  is  essentially  a  receiving  device,  the  direct 
approach  would  be  to  consider  the  response  of  the  antennas  for 
plane  waves  incident  from  various  directions  (bearings).  However, 
this  would  require  that  the  fields  of  the  incident  wave  at  the  DF 
system  antennas  be  computed  in  the  presence  of  the  aircraft  struc¬ 
ture.  That  is,  scattering  of  an  incident  plane  wave  by  the  air¬ 
craft  structure  would  have  to  be  accounted  for.  The  latter  is  an 
extremely  difficult  problem.  On  the  other  hand,  very  efficient 
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methods  exist  3*4’5*6  which  enable  one  to  compute  the  fields  of 
the  DF  antennas  if  we  consider  them  to  be  transmitting  in  the 
presence  of  the  aircraft  structure.  Using  the  reciprocity  theorem 
of  electromagnetic  theory,  a  method  can  be  developed  whe:  eby  these 
"transmit  mode"  fields  are  used  to  determine  the  response  (receive 
voltages)  of  the  DF  system  to  incident  waves,  in  the  presence  of 
the  aircraft. 

In  practice  it  is  more  than  reasonable  to  assume  that  the  source 
of  radiation  (beacon)  whose  bearing  the  DF  system  is  to  determine, 
is  in  the  far  zone  of  the  DF  system,  and  vice  versa.  The  beacon 
may  then  be  represented  as  an  infinitesimal  electric  current 
source.  Since  most  beacons  of  the  "distress"  type  will  have 
electrically  small  omnidirectional  antennas,  such  a  representation 
as  an  infinitesimal  source  will  usually  be  valid  irrespective  of 
range. 

2.2  Definitions  and  Terminology 

Consider  the  configuration  shown  in  Fig  1(a)  which  shows 
schematically  the  n-th  antenna  (antenna  "a")  of  a  DF  system 
radiating  in  the  presence  of  an  aircraft  structure. 

At  location  R.  *  Rs  =  (x5,ys,zs)  is  a  beacon  (antenna  "b") 
of  orientation  hs,  which  in  the  light  of  the  preceeding 
discussion  can  be  considered  transverse  to  J*s.  The  beacon 
represents  an  infinitesimal  electric  current  source  of  length  M 
and  strength  I q.  Thus  its  current  density  can  be  written  as. 


Otj(x,y,z)  =  I0  aa  6(x-xs,  y-ys,  z-zs)  ns 


(1) 


The  orientation  hs  selected  determines  the  polarisation  of  the 
incident  field  being  simulated. 

The  surface  Si  encloses  antenna  "a",  while  S2  encloses  both 
antennas  "a"  and  the  aircraft. 

The  fields  (Ea,Ha)  are  those  of  antenna  "a"  (the  n-th  UF 
antenna)  when  It  “is  transmitting  in  the  presence  of  the  aircraft 
structure  and  the  antenna  "b".  Thus  the  scattering  from  both  the 
aircraft  and  that  from  antenna  "b"  is  incluoeo.  Since  "b"  is  the 
infinitesimal  element,  its  effect  on  (Ea,Ha)  wiil  be  neglig¬ 
ible. 


The  fields  Ub,Hb)  are  those  of  antenna  "b"  (the  beacon)  vhen 
it  is  transmitting  in  the  presence  of  the  aircraft  structure  and 
antenna  "a".  Thus  scattering  from  both  the  aircraft  and  antenna 
"a"  is  included. 

(Ep,Hp)  represents  the  fields  Ua»Ha)  as  observed  at 
R  a  Rs,  while  (Ej*^t)  represents  (C^,Ha)  as  observed  on 
"S2.  “Finally,  will  be  f i elds  lEb.Hb)  as  observed 

on  S2. 

2.3  Reciprocity  Theorem 

The  general  reciprocity  theorems  for  an  electromagnetic 
field  states  that  , 


f,:  ( Ha  x  Hb  “  Hb  x  ila )  * 

=  /y(ib  •  ia  -  la  *  Jb  -  «b  •  Ma  +  Ha  *  Mb)  dV  (2) 


where  the  closed  surface  S  encloses  the  volume  V  containing  two 
sets  of  electric  and  magnetic  sources  (Ja>Ma)>  (Hb»Hb) 
which  give  rise  to  electromagnetic  fields  (E.a»lia)  and 
(Eb,Hb)  respectively. 

The  vector  dS  is  normal  to  surface  S  and  is  directed  positively 
out  of  S.  ” ** 

2 . 4  Application  of  the  Reciprocity  Theorem 


The  reciprocity  theorem  given  in  expression  (2)  will  now  be 
applied  in  two  ways,  to  the  configuration  of  Fig  1(a). 

First,  choose  the  surface  S  in  the  theorem  as  Si  +  S2,  so  that  the 
volume  V  enclosed  by  S  is  a  source-free  region,  containing  neither 
antenna  "a"  or  antenna  "b".  In  such  a  case,  the  volume  integral 
on  the  right  hand  side  of  equation  (2)  is  zero.  The  surface 
integral  on  the  left  side  can  then  be  separated  into  two  surface 
integrals  as. 


E.x  H  J  •  dS 
— d  —a  — 


x  h  )  •  dS 

— a 


(3) 
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with  d$  directed  positively  out  of  V.  Using  the  "shielded 
antenna*,  whereby  surface  Si  enclosing  the  antenna  can  be  con¬ 
sidered  (for  the  wire  antenna  case  in  Fig  1(a)  and  the  aperture 
antenna  case  of  Fig  1(b))  to  coincide  with  conducting  portions  of 
the  antenna  except  for  an  "aperture"  portion  Sap,  Collin  8  has 
shown  that  the  received  voltage  at  antenna  "a"  due  to  source  "b" 
is 


VR  ■  Cn  f  (E,  x  Hp  -  Efc  X  Hj)  •  dS 
Si 


(4) 


where  Cn  is  a  complex  constant  dependent  on  the  particular 
antenna  type.  From  (3)  this  voltage  is  then  also  given  by, 


Vk  ■  -Cn*.(Ea  *  i!b  *  b  *  «a>  •  dS 
52 


(5) 


If  ( E j.JHj )  and  (Eft,Hj<)  on  S2  are  substituted  for 
(£a,Ha)  and  (E^.Hp)  respectively,  this  is  just 


VK  =  -Cn^  (It  x  Hr  -  Er  x  Ht)  •  dS 
52 


(6) 


Let  the  reciprocity  theorem  (2)  be  applied  once  more  with  volume  V 
being  the  region  enclosed  by  surface  $2  and  the  sphere  at  infinity 
Sa>.  For  the  volume  integral  term  on  the  right  hand  side  of  (2), 
the  only  source  term  enclosed  is  the  beacon  antenna  "b",  with 
current  density  ^Jb  given  by  (1).  It  therefore  integrates  to 
-Ef(Rs).  hs  Iq  M,  where  Ef  is  the  field  Ia  as  observed 
aIRs,  as  decided  above.  Thus  expression  (2)  becomes. 


i  Ua  x  Hb  -  Eb  x  Ha)  •  dS  =  -EF(RS)  •  ns  1 0  A* 

S2+5m 
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But  that  portion  of  the  surface  integral  in  (7)  over  the  infinite 
sphere  is  zero  .  With  (Ea,Ha)  and  (Ej,,Hjj)  on  S2  replaced 
by  the  symbols  (Ej,Hy)  “an?  FT  respectively,  (7) 
becomes 


■L  (ET  x  Hr  -  1r  *  HT)  •  dS  ■  -Er (R*>  •  ns  I  o  AA  (8) 

S  2 


Comparison  of  (6)  and  (8)  then  shows  that  the  voltage  received  by 
the  n-th  DF  antenna,  when  beacon  "b"  is  transmitting  is. 


Vr  -  Cn  Ep  (Rs)  •  r»s  I0  At 


(9) 


Note  that  (9)  is  valid  for  any  range  and  orientation  of  the 
beacon. 


2 . 5  Summary  of  Theoretical  Result  for  Received  Voltages 

The  voltage  received  by  the  n-th  antenna  of  the  DF  system, 
when  the  beacon  of  orientation  ns  is  transmitting  from  bearing 
R.s»  is  given  by 


vn  s  Ao  ns  *  If  (is) 


where 

(i)  Ao  is  a  complex  constant,  Cn  Io  At,  which  will  be  of  no 
concern  in  determining  the  relative  complex  values  of  the 
receive  voltages  of  the  identical  DF  antennas. 

(ii)  Ep(Rs)  is  the  field  at  Rs  of  the  n-th  antenna  of  the 
DF  system,  when  it  is  transmitting  in  the  presence  of  the 
aircraft.  These  fields  are  computed  using  the  geometrical 
theory  of  diffraction  (GTD)  methods  ’  ’ . 

In  order  to  evaluate  the  DF  system,  the  fields  Ep(Ks)  of  each 
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DF  antenna ,  with  phase  referred  to  a  common  reference  point,  are 
computed  for  a  known  bearing  (ns  location  Rs).  Then  the  vol¬ 
tages  Vfl  are  computed  for  each  antenna  using  (10)  and  used  in 
the  particular  UF  system  algorithm  to  find  the  predicted  bearing. 
The  bearing  error  is  then  the  difference  between  the  known  and 
predicted  bearings.  This  is  repeated  for  each  direction  of  inci¬ 
dence  of  interest. 

3.  NUMERICAL  PREDICTION 

3.1  Introduction 

The  6TD  is  ideal  for  studying  antennas  in  a  complex  envi¬ 
ronment  as  only  the  most  basic  features  of  a  very  complex  struc¬ 
ture  need  to  be  modelled.  The  components  of  the  fields  incident 
on,  reflected  by  and  diffracted  by  the  various  structures  are 
determined  by  using  ray  optical  techniques.  Components  of  the 
fields  are  found  by  summing  the  Geometrical  Optics  (GO)  terms  with 
the  diffracted  terms  in  the  far  field.  The  diffraction  terms  can 
be  found  by  using  the  GTO  solutions  for  the  individual  rays  inci¬ 
dent  on  structural  elements  such  as  edges  and  corners.  The  rays 
from  any  one  structural  element  tend  to  interact  with  other  struc¬ 
tural  elements  and  cause  various  higher  order  terms.  By  tracing 
all  the  possible  rays  and  evaluating  only  the  dominant  terms  an 
accurate  and  efficient  code  can  be  created  3’4’b,b. 

3.2  Implementation 

A  short  description  of  how  the  GTD  is  implemented  in  the 
code  used  is  now  given.  Th.  code  allows  the  moael  used  to  include 
a  finite  elliptical  cylinder  and  a  set  of  flat  plates.  The 
cylincer  may  be  capped  with  a  flat  plate  which  may  be  at  any  angle 
relative  to  the  cylinder  in  the  vertical  plane  4.  The  plates  must 
be  flat  and  each  plate  may  have  up  to  six  corners. 

The  following  terms  are  included  in  the  code: 

direct  field 

field  reflected  from  a  plate 
field  doubly  reflected  by  plates 
field  diffracted  by  plates 

field  reflected  by  a  plate  then  diffracted  by  a  plate 
field  diffracted  by  a  plate  then  reflected  by  a  plate 

field  scattered  by  the  cylinder 
field  reflected  by  an  end  cap 
field  diffracted  by  an  end  cap  rim 
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field  reflected  by  a  plate  then  scattered  by  the  cylinder 
field  scattered  by  the  cylinder  then  reflected  by  a  plate 
field  reflected  by  the  cylinder  then  diffracted  by  a  plate 
field  diffracted  by  a  plate  then  reflected  by  the  cylinder 

Scattered  by  cylinder  means  that  the  fields  scattered  by  the 
cylinder  includes  both  reflected  and  diffracted  (creeping)  rays. 

The  code  uses  the  method  discussed  in  the  previous  section  to 
determine  receive  voltages  of  the  DF  antennas  in  the  presence  of 
the  aircraft.  As  such  the  code  is  a  transmitting  antenna  code  and 
it  computes  the  transmit  radiation  pattern  of  the  antennas  mounted 
on  the  structure. 

The  code  allows  more  than  one  antenna  to  be  mounted  on  the  struc¬ 
ture  at  any  one  time.  The  free  space  radiation  pattern  of  each 
antenna  on  the  structure  must  be  known  a  priori.  Because  of  the 
simplicity  of  expressions  for  the  radiation  patterns  of  dipole 
(monopole)  and  slot  antennas  and  the  fact  that  their  fields  have 
the  required  ray  optical  form,  they  are  used  in  the  code. 

To  compute  the  radiation  pattern  of  a  two  antenna  DF  system,  the 
program  is  used  in  the  transmit  mode  and  both  antennas  are  trans¬ 
mitting  simultaneously  with  a  phase  difference  of  180  degrees. 
The  total  radiation  pattern  of  the  two  antennas  is  thus  found  for 
comparison  with  measurements. 

To  compute  the  errors  introduced  in  a  DF  system  the  transmit 
radiation  pattern  is  found  for  each  antenna.  The  method 
described  above  is  then  used  to  compute  the  receive  voltages  of 
each  antenna  from  its  transmit  fields  with  amplitude  and  phase 
information  preserved.  The  receive  voltages  of  each  of  the 
antennas  can  then  be  fed  into  the  OF  algorithm  used  to  compute  the 
errors  introduced  because  of  the  distortion  of  the  receiving 
patterns,  The  algorithm  of  a  four  antenna  system  described  by 
Gething  }  has  been  included  to  supplement  the  code.  This  enables 
direct  calculation  of  DF  errors  for  an  interferometric  system. 

In  order  to  establish  that  the  GTD  method  described  correctly 
preserves  the  phase  information  so  essential  for  application  to 
interferometric  DF  applications,  a  DF  array  of  two  antennas 
mounted  on  a  box  was  first  modelled.  Rearing  errors  were  then 
calculated  in  two  ways.  First  the  GTD  w<  s  used  to  find  the 
transmit  fields  and  the  bearing  errors  found.  Then  the  method  of 
moments  was  used  to  find  the  transmit  fields  and  the  bearing 
errors  found.  The  results  were,  for  all  intents  and  purposes, 
identical.  In  what  follows  only  the  GTD  has  been  used  in  the 
computations. 


-663- 


3.3  The  Numerical  Aircraft  Models 


The  models  under  consideration  have  been  built-up  by  using 
a  circular  cylinder  with  capped  ends  for  the  fuselage.  The  wings 
were  modelled  with  one  flat  plate  each  and  the  tail  was  modelled 
with  a  set  of  plates.  The  modelling  of  the  engines  was  done 
using  one  of  two  methods.  The  first  engine  model  was  a  vertical 
plate  hanging  from  the  wing.  This  model  will  have  the  correct 
shadowing  properties  but  does  not  have  the  correct  diffraction 
effects  of  a  real  engine.  For  a  more  realistic  model  of  the 
engine  three  plates  were  used.  Because  the  number  of  plates  that 
may  be  used  in  the  code  is  limited  to  14,  not  all  the  engines  of 
the  Hercules  C130  could  be  modelled  in  this  way. 

4.  EXPERIMENTAL  PROCEDURE 

4.1  Description  of  Models 

4.1.1  Lockheed  C130  Type  L1Q0-30 

The  choice  of  scale  was  determined  by  the  availability  of 
aluminium  tubing  which  happened  to  have  a  diameter  of  220  mm. 
Flat  plate  with  a  thickness  of  3  mm  was  used  for  all  wing  sur¬ 
faces.  Small  cylinders  were  mounted  underneath  the  wings  to  model 
the  engine  pods.  For  convenience  a  20:1  scale  factor  was  used. 
Details  of  the  model  are  given  in  Table  2  together  with  the  posi¬ 
tions  cf  the  Adcock  arrays.  A  hemisphere  was  used  for  the  nose 
section.  The  upwards  tapered  tail  section  of  the  (.130  was 
approximated  by  a  conic  section,  the  base  of  which  had  been  cut  at 
a  slight  angle  to  provide  the  necessary  upward  slope.  All  slots 
formed  were  taped  with  aluminium  tape.  Five  possible 
configurations  were  considered.  These  were: 

(11  Fuselage  only 

(2)  Fuselage  +  Tail  plane 

(3)  Fuselage  +  Wings  +  Tailplane 

(4)  fuselage  +  Wings  +  Tailplane  +  Inner  Engine  Pods 

(5)  Full  model  (excluding  rudder) 

These  configurations  are  adequately  described  by  the  stick  model 
shown  on  the  polar  diagrams.  Polar  diagram  measurements  were 
performed  for  two  antenna  positions  (1  or  2)  with  the  Adcock  array 
mounted  circumferentially  (C)  or  axially  (A)  along  the  centre  line 
on  the  underbelly.  The  mounting  designations  are  shown  on  the 
diagrams  as  1A,  1C,  2A  or  2C. 

A  photo  of  the  model  and  one  of  the  authors  (G.W.R.)  is  shown  in 
Figure  2.  For  the  measurements  the  model  was  mounted  upside  down 
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TABLE  2 


Construction  details  of  C130  Type  L100-3Q 


ercules  aircraft  mode 


Scale 


20:1  (approximate) 


Fuselage  :  Total  length  155  cm 

Nose  section  (hemisphere)  11  cm 

Central  section  (cylinder)  94  cm 

Tail  section  (cone  with  base  cut  at  angle)  50  cm 

Di ameter  22  cm 


Wings  :  Total  span  204  cm 

Chord  length  at  tip  12  cm 

Chord  length  at  root  24.5  cm 

Length  central  untapered  section  61  cm 

Distance  of  leading  edge  from  nose  G8  cm 

Tail  planes:  Total  Span  83  cm 

Chord  length  at  tip  9  cm 

Chord  length  at  root  20  cm 

Distance  of  leading  edge  from  nose  131  cm 

Engines(4):  Length  of  cylinders  25  cm 

Diameter  of  cylinders  3.7  cm 

Protrusion  ahead  of  wing  11.5  cm 

Propel  lor  length  (fixed  vertical  ly  for 
maximum  scattering  effect)  20  cm 

Engines  fitted  directly  to  the  wing: 

Distance  between  centres  of  inner 

cylinders  50  cm 

Distance  between  centres  of  outer 

cylinders  100  cm 


Adcock  array: 

Position  1  (from  nose)  33  cm 

Position  2  (from  nose)  55  cm 

Aluminium  used  throughout 


on  the  antenna  positioner.  A  bottom  view  line-drawing  is  shown  in 
Figure  3a. 

4.1.2  The  Low  Wing  Jetliner 

A  2.2  m  length  of  the  aluminium  tubing  was  left  over  after 
fabrication  of  the  C130  model.  This  was  used  to  make  a  basic 
model  of  a  low-wing  jetliner.  The  dimensions  are  shown  in  Table  3 
together  with  the  Adcock  array  positions.  Two  cylinders  were 
attached  directly  to  and  underneath  the  wings  to  model  the  engine 
pods. 

After  completion  of  the  series  of  measurements,  portions  of  the 
fore  and  after  sections  were  cut  off  and  replaced  by  a  hemisphere 
and  a  conic  section  so  that  the  model  fuselage  would  resemble  an 
actual  fuselage  more  closely.  A  new  tailplane  with  the  same 
overall  dimensions  was  made  to  cater  for  the  wing  root  at  the 
conic  section.  The  overall  length  of  the  improved  model  was  the 
same  as  that  of  the  basic  cylindrical  model  and  a  bottom  view  is 
shown  ir.  Figure  3b.  These  models  are  designated  as  I  and  R  for 
improved  and  basic  respectively.  As  before,  the  Adcock  arrays 
were  mounted  in  two  positions  either  circumferentially  or  axially 
along  the  central  line  on  the  underbelly.  This  model  was  also 
mounted  upside  down  for  the  measurements.  During  the  course  of 
measurements  it  was  found  that  the  welded  bracket  supporting  the 
model  was  slightly  off  centre.  This  resulted  in  the  port  and 
starboard  wings  making  dihedral  angles  of  5°  and  7°  respectively. 
As  for  the  C130  model  a  number  of  configurations  were  possible. 

4.2  The  Adcock  Array 

The  basic  Adcock  array  consisted  of  two  short  monopoles  fed 
via  6  dR  attenuators  to  the  inputs  of  a  180°  hybrid  to  form  sum  or 
difference  patterns. 

The  probes  were  two  pins  soldered  into  the  central  conduc¬ 
tors  of  bulkhead  mounted  SMA  connectors.  The  pins  projected 
through  holes  in  the  fuselage  approximately  2.5  mm  in  diameter. 
For  the  frequency  range  above  1,6  GHz  the  projecting  length  was 
1  cm.  Relow  this  frequency  the  projecting  length  was  2  cm.  The 
spacing  between  the  probes  was  kept  fixed  at  3  cm.  This  was  the 
minimum  practical  value  in  view  of  the  fact  that  the  SMA 
connectors  and  other  hardware  were  fitted  inside  the  cylinder. 
Allowance  had  to  be  made  for  connecting  the  180°  hybrids  with  very 
short  lengths  of  relatively  inflexible  coaxial  cable.  The  gain 
of  the  individual  elements  was  not  expected  to  affect  the  overall 
results  provided  sufficient  signal  to  noise  ratios  could  be 
obtained.  For  future  reference  the  gains  of  the  elements  were 
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TABLE  3 


Construction  details  of  model  to  represent  a  low-wing  jetliner 


Scale:  20:1  (approximately) 

Fuselage:  Basic  model  (designation  B) 

Cylinder  length  (Ends  blocked  with  metal 
disks) 

Fuselage:  Improved  model  (designation  I) 

Total  length 

Nose  section  (hemisphere) 

Tail  section  (cone) 

Central  cylindrical  section 


Diameter  of  both  models: 


kings:  Length  each  wing 

Chord  length  at  tip 
Chord  length  at  root 
fitted  to  give  total  span  of 
Sweep  at  leading  edge 
Dihedral  angle 

Distance  of  leading  edge  from  nose 

Tail  planes-  Length 

Chord  length  at  tip 

Chord  length  at  root 

Sweep  at  leading  edge 

Distance  of  leading  edge  from  nose 


Engines(2):  Length  of  cylinders 

Diameter  of  cylinders 
Distance  between  centres  of  cylinders 
Protrusion  ahead  of  wing 
Fitted  directly  to  the  wing 


Adcock  array: 

Position  1  (from  nose) 
Position  2  (from  nose) 

Aluminium  used  throughout 


220  cm 


220  cm 
11  cm 
50  cm 
159  cm 

22  cm 

107  cm 
14.2  cm 
47.5  cm 
222  cm 
37° 

6° 

79  cm 

24  cm 
10  cm 

25  cm 
43° 

195  cm 

26  cm 
3.5  cm 
100  cm 

17  cm 


60  cm 
45  cm 


I 
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measured  at  1,6  GHz  with  the  6  dR  attenuator  pads  in  place.  The 
gains  were  found  to  be  -20  dRi  and  -24  dRi  for  the  2  cm  and  1  cm 
probes  respectively. 

All  interconnections  were  made  with  RG  142  R/U  cable  terminated  in 
SMA  connectors.  In  order  to  minimise  phase  errors  the  electrical 
lengths  were  checked  using  a  Hewlett  Packard  2-18  GHz  Network  Ana¬ 
lyser.  Appropriate  cable  pairs  were  selected  for  the  closest 
possible  match. 

The  180°  hybrid  couplers  were  a  MARDA  Model  4343  covering  the 
range  2-8  GHz  and  an  ANZAC  Model  H-183-4  covering  the  range  30  to 
3000  MHz.  Measurements  were  made  between  2  and  3  GHz  using  both 
couplers  to  obtain  an  overlap  of  results.  The  amplitude  and 
phase  imbalance  for  the  NARDA  Hybrid  are  given  as  ±0.5  dR  max  and 
±7  degrees  max  respectively  with  an  isolation  of  17  dR.  The 
amplitude  and  phase  imbalance  of  the  ANZAC  hybrid  are  given  as 
±0.4  dB  and  7.5  degrees  max  from  100  to  1500  MHz  rising  to  15 
degrees  max  from  1500  to  3000  MHz  respectively.  The  isolation  was 
given  as  20  dR  min  for  the  ANZAC  hybrid.  Roth  hybrids  were 
checked  using  the  HP  network  analyser  and  found  to  be  well  within 
specification.  The  worst  phase  imbalance  measured  did  not  exceed 
approximately  7.5  degrees.  This  implies  that  broadside  null 
depths  in  excess  of  17  to  20  dR  should  be  treated  with  caution. 
In  general  the  null  depths  exceeded  20  dR. 

The  6  dR  attenuators  ensured  that  the  mismatch  between  the  inputs 
of  the  hybrids  and  the  probes  would  not  exceed  an  SWR  of  2:1. 

The  probes  were  mounted  on  the  rounded  cylindrical  surface  as  is. 
No  attempt  was  made  to  provide  a  flat  groundplane  in  the  vicinity 
of  the  array.  In  the  circumferential  position  this  meant  that  the 
maximum  surface  height  between  the  probes  due  to  base  curvature 
for  the  array  was  1  mm  higher  than  the  line  interconnecting  the 
holes  through  which  the  probes  project.  This  was  not  believed  to 
be  detrimental  to  the  measurements  and  simplifies  the  numerical 
model . 

Access  to  the  bulkhead  mounted  SMA  connectors,  attenuators  and 
hardware  was  through  a  hatch  measuring  25  cm  axially  and  12  cm 
circumferentially  on  the  opposite  side  of  the  cylinder  from  the 
array  itself.  A  single  cable  carrying  the  difference  signal  was 
led  through  a  small  cutout  in  the  hatch. 

4.3  The  Antenna  Test  Facility 

All  measurements  were  performed  at  the  National  Antenna 
Test  Facility  operated  by  the  National  Institute  for  Telecommuni- 
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cations  Research  of  the  South  African  Council  for  Scientific  and 
Industrial  Research.  This  facility  operates  in  the  ground  reflec¬ 
tion  mode  and  is  situated  approximately  30  km  Northeast  of 
Pretoria.  It  is  described  fully  in  the  reports  by  Raker  10  . 

The  antenna  positioner  is  mounted  on  a  concrete  mast  at  the  con¬ 
trol  room.  The  total  height  from  the  ground  to  the  model  was 
approximately  12  m.  Down  range  slope  for  the  ground  reflection 
mode  was  1.4  degrees.  The  source  tower  is  500  m  away.  Source 
antenna  heights  can  be  adjusted  from  0,3  m  to  a  maximum  cf  11  m 
depending  on  the  frequency.  The  source  positioning  system  allows 
adjustment  of  polarisation,  azimuth  and  elevation  pointing  angle. 
A  range  of  dish  antennas  is  used  to  cover  various  frequency 
ranges.  Range  instrumentation  allows  operation  from  500  MHz  to  18 
GHz. 

For  frequencies  above  1  GHz  a  wideband  mixer  must  be  used.  Sig¬ 
nals  with  frequencies  less  than  this  couple  directly  to  the  recei¬ 
ver. 


5.  RESULTS 


All  measured  results  discussed  here  were  digitised  in  one 
degree  increments  in  order  to  display  them  on  the  same  polar  plots 
as  the  predicted  values. 

Figures  4  and  5  show  a  sequence  of  measured  polar  diagrams  with 
increasing  model  complexity  for  the  C130  Hercules  aircraft  at  2 
and  6  GHz.  The  configurations  are  adequately  described  by  the 
line  drawings.  Only  results  for  a  basic  Adcock  array  mounted 
axially  along  the  fuselage  are  shown.  The  nulls  are  thus 
broadside  to  the  fuselage.  Also  plotted  are  the  results  predicted 
by  the  modified  GTD  code. 

At  2  GHz  there  is  an  obvious  bulge  in  the  polar  diagram  in  the 
direction  of  the  tail.  This  is  ascribed  to  diffraction  effects 
and  creeping  waves  along  the  conical  section  of  the  fuselage.  The 
fact  that  the  Adcock  array  is  positioned  near  the  nose  and  the 
fact  that  this  effect  was  observed  at  other  frequencies  suggests 
that  the  bulge  cannot  be  ascribed  to  constructive  interference  of 
diffracted  and  direct  fields  alone.  A  more  pronounced  effect  in 
the  direction  of  the  nose  would  haVe  been  expected.  The  effect  of 
increasing  model  complexity  is  surprisingly  little  as  can  be  seen 
by  comparing  Figure  4(a)  and  4(d).  Some  increase  in  pattern 
ripple  is  noticeable. 

The  predicted  results  show  an  anomaly  towards  the  nose  from  +15  to 
-15  degrees  as  well  as  one  in  the  direction  of  the  tail.  It  is 
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expected  that  the  addition  of  one  further  GTD  ray  term  to  the 
summation  could  correct  this.  Of  interest  is  the  prediction  of 
the  increasing  bulge  in  the  polar  diagram  in  the  direction  of  the 
tail.  The  numerical  model  predicts  somewhat  greater  effects  for 
the  presence  of  wings  and  engine  pods.  These  effects  are  not 
obvious  in  the  measured  results.  The  agreement  is  regarded  as 
fair. 

At  6  GHz  the  measured  pattern  is  remarkably  clean  with  very  little 
ripple  effect  as  can  be  seen  in  Figures  5(a)  to  (d).  Agreement 
between  the  patterns  for  the  fuselage  alone  and  the  full  model  is 
generally  within  1  dB  or  better.  The  ripple  effect  is  consider¬ 
ably  less  pronounced  than  for  the  2  GHz  case.  The  effects  of 
wings  and  engine  pods  should  have  appeared  between  110°  and  140° 
as  measured . from  the  nose  but  aro  not  evident. 

The  agreement  between  predicted  and  measured  results  is  satisfac¬ 
tory  except  for  the  same  anomalous  behaviour  towards  the  nose  and 
tail  as  was  the  case  at  2  GHz. 

From  the  point  of  view  of  having  scattering  or  diffraction  effects 
for  purposes  of  checking  the  code  these  results  were  somewhat  dis¬ 
appointing  and  gave  rise  to  the  next  series  of  measurements  on  a 
model  representing  a  low-wing  jetliner  with  two  underwing  engine 
pods. 

For  purposes  of  illustration  the  sequence  of  measurements  with  the 
basic  (blunt  cylinder)  model  at  7  GHz  are  shown  in  Figure  6  for  an 
Adcock  array  mounted  circumferenti al ly  in  position  1.  From  other 
measurements  it  was  clear  that  the  tailplanes  had  minimal  effect 
and  this  sequence  illustrates  the  effect  of  wings  and  engine  pods 
only.  A  slight  wing  effect  appears  to  be  evident  in  the  top  half 
of  the  polar  diagram  in  Figure  6(c)  about  130°  from  the  nose 
direction.  Reference  was  made  earlier  to  the  unequal  dihedral 
angles.  This  difference  was  found  to  account  for  considerable 
asymmetry  in  many  of  the  measurements.  In  Figure  6(d)  the  effect 
of  engine  pod  shadowing  and  diffraction  can  be  clearly  seen  in  the 
same  direction  noted  previously  and  to  a  lesser  extent  in  the 
lower  half.  It  was  generally  found  that  the  amount  of  ripple  in 
the  polarplot  was  substantially  higher  than  in  the  case  of  the 
C130  model . 

Finally  Figure  6(f)  shows  a  comparison  of  measured  results  for  the 
complete  model  with  the  blunt  cylinder  and  the  improved  cylindri¬ 
cal  body  as  fuselage.  Apart  from  what  appears  to  be  a  2  degree 
offset  (due  to  misalignment  of  the  models)  and  a  slippage  in  the 
positions  of  ripples  the  overall  agreement  is  regarded  as  good. 
Note  that  the  points  of  normalisation  for  the  two  curves  are  not 
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the  same.  Asymmetry  in  the  individual  curves  can  be  ascribed  to 
the  model  asymmetr y  and  imbalance  in  the  180  degree  hybrids. 

There  is  some  discrepancy  between  the  measured  and  predicted 
results  as  shown  in  Figures  6(a)  and  6(b)  and  agreement  is 
regarded  as  fair.  There  is  far  worse  agreement  between  the  curves 
shown  in  Figure  6(c)  and  considerable  asymmetry  in  the  predicted 
results.  Time  did  not  permit  a  check  of  the  numeric  model  used 
although  it  is  suspected  that  there  may  be  an  error  concerning  one 
of  the  wing's  data.  Figure  6(d)  shows  that  the  numerical  model 
predicts  the  effects  of  the  enginepods  in  approximately  the 
correct  position. 

It  was  found  that  using  a  flat  vertical  plate  or  a  combination  of 
three  plates  attached  to  the  wing  to  form  a  box  section  made  very 
little  real  difference  to  the  predicted  results. 

On  the  strength  of  these  results  the  expected  DF  errors  for  the 
configurations  shown  were  calculated  assuming  that  the  DF  system 
was  of  the  interferometric  type. 

The  predicted  DF  errors  are  shown  in  Figure  7  for  the  C130  model 
and  Figure  8  for  the  low-wing  jetliner  model.  Many  of  the  spikes 
can  be  attributed  to  irregularities  which  could  arise  from  the 
lack  of  a  sufficient  number  of  GTD  terms  rays  in  the  numeric 
solution.  It  would  appear  though  that  the  high-wing  C130  would 
make  a  better  DF  platform  than  a  low-wing  aircraft.  Further  work 
is  however  needed  torefine  the  numerical  model  for  the  analysis  of 
DF  systems  on  aircraft. 

6.  CONCLUSIONS 


DF  arrays  mounted  on  top  of  aircraft  fuselages  suffer 
greater  induced  errors  from  stabiliser  effects  than  those  mounted 
on  the  underbelly.  For  the  latter  case  the  DF  errors  are  least 
for  aircraft  with  high  wing  mountings  and  the  error  appears  to 
decrease  with  frequency  due  to  decreased  distortion  of  the  basic 
Adcock  array  patterns  by  scattered  fields.  For  low-wing  aircraft 
the  DF  errors  will  be  less  for  those  with  the  engines  mounted  in 
the  tail  than  for  those  with  underwing  pods. 

The  results  obtained  from  numerical  modelling  are  encouraging.  It 
is  believed  that  further  work  will  yield  in  predictions  which  are 
sufficiently  accurate  for  the  evaluation  of  the  performance  of  DF 
arrays  on  various  aircraft  as  well  as  the  optimum  positioning  of 
such  arrays  on  specific  aircraft. 

The  relatively  good  agreement  between  the  basic  model  using  a 


-671- 


I 


cylindrical  fuselage  and  the  improved  model  using  a  hemispherical 
nosecap  and  a  conical  tailsection  suggests  considerable  simplifi¬ 
cation  in  the  physical  and  numerical  model  of  an  aircraft  in  the 
VHF  range. 
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for  receive  voltages  on  a  l)F  array.  (Wire 
antenna  case) 

(b)  Surface  of  integration  for  aperture  type  ancenr.a. 


are  approximately  tc  scale.  The  basic  Adcock  arrays  were  fixed  at 
positions  1  or  2  either  axially  (A)  along  the  fuselage  or  circumferen¬ 
tially  (C)  across  it. 
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1  on  the  underbelly  of  a  20:1  scale  model  of  a  Cl 30  Hercules  transport 
aircraft.  The  configurations  are  represented  by  the  stick  drawings. 
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neas'  ^ea  ^sona  line;  ana  predicted  iDroken  line;  receive  mode  polar 
diagrams  at  6  GHz  of  a  simple  Adcock  array  mounted  axially  at  position  1 
on  the  underbelly  of  a  20:1  scale  model  of  a  Cl 30  Hercules  transport 
aircraft.  The  configurations  are  represented  by  the  stick  drawings. 
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polar  diagrams  at  7  GHz  of  a  simple  Adcock  array  mounted  circum¬ 
ferentially  in  position  1  on  the  underbelly  of  a  20:1  scale  model 
representation  of  a  low-wing  jetliner  with  two  underwing  engines. 
The  configurations  are  represented  by  the  stock  drawings.  A 
blunt  cylinder  was  used  for  the  fuselage. 
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the  fuselage. 


ERROR  ANGLE  (DEGREES. )  v  ERROR  ANGLE  (DEGREES. ) 


(b)  ANGLE  (DEGREES) 

FIGURE  7.  Predicted  UF  errors  for  the  Hercules  transport  model 

aircraft  at  2  (a)  and  6  (b)  assuming  an  interferometric 
UF  system. 
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ERROR  ANGLE  (DEGREES) 


FIGURE  8.  Predicted  JF  errors  for  the  giodel  low-wing  jetliner 
at  7  GHz. 
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MATRIX  FORMULATION  OF  VECTOR  OPERATIONS 
IN  ELECTROMAGNETICS  ANALYSIS 
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Feng  Ch^ng  Chang 
Antenna  Systems  Laboratory 
TRW  Electronic  Systems  Group 
Redondo  Beach#  California 
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ABSTRACT 

The  manipulation  of  vector  quantities  is  unavoidable  in 
the  study  of  electromagnetic  fields.  An  efficient  technique 
is  developed  to  simplify  the  evaluation  of  vector  algebraic 
and  differential  operations  by  simply  converting  the  vector 
operations  into  matrix  operations.  This  matrix  formulation 
is  especially  useful  when  there  are  mixed  coordinate  basis 
involved  in  the  vector  operations.  Applications  to  some 
antenna  analysis  problems  will  demonstrate  the  merit  of  the 
technique  presented.  ^ 


MATRIX  FORMULATION 

In  analysis  it  is  usually  more  efficient  to  work  with 
matrix  operations  than  to  deal  with  vector  operations. 
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In  vector  operation,  components  and  basis  of  a  vector  are 
inseparate  and  must  work  together;  whereas  in  matrix 
operation,  components  and  basis  can  be  dealt  with  separately 
and  usually  only  the  vector  components  are  involved. 
Therefore,  most  vector  operations  can  be  handled  more 
efficiently  by  simply  converting  them  into  matrix  operations. 

In  matrix  formulation,  a  vector  A  is  expressed  as  a 
column  matrix  with  components  (  A^,  A2,  A3  )  where  it  is 
understood  that  it  has  (  Uj,  u2,  u3  )  as  its  basis.  Four 
basic  conversions  are  required  to  relate  vector  to  matrix 
expressions: 

(1)  a  scalar  quantity  as  a  diagonal  matrix, 

(2)  a  vector  or  gradient  operator  as  a  column  matrix, 

(3)  a  dot-product  or  divergent  operator  as  a  row  matrix, 

(4)  a  cross-product  or  curl  operator  as  a  square  matrix. 

The  notations  \  \  ,  (  ) ,  {  } ,  and  [  ] ,  will  be 

used  to  denote,  respectively,  the  diagonal  matrix,  column 
matrix,  row  matrix,  and  square  matrix.  Explicit  expressions 
of  matrix  formulations  for  the  basic  vector  algebraic  and 
differential  operations  in  the  three  coordinate  systems  of 
interest  are  given  in  Table  1,  where  the  superscripts  c,  d, 
and  s  are  used  to  indicate  the  Cartesian,  cylindrical,  and 
spherical  coordinates. 


TABLE  1 


MATRIX  FORMULATION  OF  VECTOR  ALGEBRAIC 


AND  DIFFERENTIAL  OPERATORS 


Gradient 


Divarganct 


Curl 


7  x 


Scalar  Laplacian  V2  -  7*7 


•  *X  ■ 

Vactor  X 

(  A®  )  - 

»y 

.  A*  . 

Dot-Product  X » 

{Xs  \  - 

[  A* 

Ay  At  ] 

‘  0 

-At  Ay 

Cross-Product  X  x 

[ACJ  - 

A* 

0  -AX 

“Ay 

Ax  0  . 

(  7C)  - 


{  VC} 


tvoj 


6i  1 

6 

oy 

L  6 

oT 


[t 


6 

~Si 

o 

6y 


6 

~S7 


It 

o 


"5a 


6 

TS» 


6 

W 

6 

‘^x 


{ v°  >  ( ve )  -  r  ac  ]  , 


AC  _  6^  .  6*1  ,  d^ 

A  *  'di5  +  S&  +  ~3i2 


7actor  Laplaclan  7*  •  7  7*  -  7*  7x 


(  vo){ vc}  -  t 70] [vc] 

r  Ac  0 

0  A° 


0 
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TABLE  1.  (CONTINUED) 


CYLINDRICAL  COORDINATES 

SPHERICAL  COORDINATES 

Ap  ‘ 

*  a 

(Ad)  * 

Ap 

<  A»  >  - 

A* 

.  **  . 

•  *r  ■ 

<Xd>  - 

[  Ap 

*»  Aj 

i  A*  >  - 

a9 

Ap 

Ar  . 

•  0 

H 

A  tp 

0 

“Ar 

At  1 

[Ad]  - 

)  -Ap 

t  A®  ]  - 

Ar 

0 

“Afl 

-  -A, 

Ap  0 

-Ap 

a9 

0  . 

■  lip 

1 

1  6 
-  T-5? 

(  yd)  ■ 

p  0* 

X 

6 

(  Vs)  » 

r«ino  6? 

J 

£_ 

6r 

J 

<Vd>  - 

[ 

p-tr 

ai 

0 

"1? 

X  d 
p  T5 
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ii 

0 

d 

~-5p 

.14 
p  ”5p 

-p'lpp 
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i  Vd  H  Vd) 

M 

* 

Ad 

1 

*  P 

a  „a 
"5pP?p 

.  X  d* 

*■  p >  6p> 

a1 

+  “5T> 

(  7d ) <  ?d )  - 

CvdJCVd] 

r  * 

kd.l_ 

P* 

_2_  d 

p>  *5p 

0 

9* 

2  d 
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0 

L 

0 

0 

Ad  ■ 

<  v“> 
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T  __x  a. 1  6  i  6  _>  i 
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t  He 
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Also  included  in  Table  1  are  Laplacian  operators  defined  as 


V2f  .  V  •  7  f  and  V2  F  ■  VV*?-7#VxF 
where  f  ■  f{r)  is  a  scalar  function  and  F  -  F(r)  is  a 

vector  function.  Except  for  the  Cartesian  coordinates,  the 
derivation  of  V2  is  very  involved  through  the  conventional 
formulas  [1] .  Applying  the  matrix  formulation  technique  will 
greatly  reduce  the  derivation  to  matrix  multiplication  and 
partial  differentiation. 

Matrix  formulation  is  especially  useful  when  there  are 
mixed  coordinate  basis  involved  in  the  vector  operation. 

The  basic  transformations  in  changing  from  one  basis  to 
another  are  given  in  Table  2.  It  should  be  noted  that 
Vv  =  Vu  in  the  vector  operation  does  not  give  (  Vv)  »  (  Vu) 
in  the  matrix  operation,  unless  v  »  u  is  in  the  same  basis. 
It  does  imply  (  Vv)  »  [VTU]  (  7U)  ,  where  [VTU]  is  a 

coordinate  transformation  matrix  relating  the  u-coordinates 
to  the  v-coordinates .  Employing  the  technique  of  coordinate 
transformation  flowchart  [2],  any  transformation  matrix  among 
the  coordinate  systems  can  readily  be  determined. 

The  following  applications  to  the  reflector  antenna 
analysis  will  demonstrate  the  merit  of  matrix  formulation 
of  vector  operations  in  the  mixed  coordinate  basis. 
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TABLE  2 


TRANSFORMATIONS  OF  COORDINATE  BASIS 


VSCTOlt  ALOtMUUC  OPUATOftS 

Vector  X1  -  1“ 

t  A*>  *  tV)  (  Au> 

Dot-Product  7C  •  ■  T1 . 

l  ifl  .  1 1»>  [Vi 

Croee-Product  7?  *  -  7?  x 

1  A*)  •  IV)  I  A«i  l«TV] 

VSCTO*  DIPPBaiHTZAL  OPERATORS 

Gradient  vv  -  v“ 

(  •  |V)  {  v“) 

Divergence  Vv  •  *  v“  . 

(  -  l  T“)  (“t*) 

Curl  Vvx  »  v“x 

I  »*)  -  |V|  I  »“i  |V| 

z* 


Figure  1 


Reflector  (Non-Star red)  and  Feed  (Starred)  Systems 


< 
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REPLICATIONS 


In  constructing  a  physical  optics  radiation  integral  to 
determine  the  vector  potential*  the  induced  current  on  the 
reflector  surface  illuminated  by  a  positioned  feed  source 
must  be  first  determined.  As  shown  in  Figure  1*  the  induced 
surface  current  J  (r)  at  a  point  r  (  0,  <P,  r  )  on  the 
reflector  surface  F  (?)  «  0  can  be  derived  by 

J  (?)  *  1  ii  (?)  x  r'#  x  E  (?'*) 

v 

The  feed  system  (starred  system)  and  the  reflector  system 
(non-starred  system)  are  related  through  both  the  origin 
translation  along  the  feed  location  vector  ?d  (  0d*  $>d ,  rd  ) 
and  the  three  axsis  rotations  about  the  feed's  positioned 
angles  (  0W*  0W»  )  •  The  radiated  field  E  ( r'")  is 

generally  given  in  the  spherical  coordinates  (  ©'**  4>'\  r '*)  of 
the  feed  system.  If  the  normal  vector  n  (?)  =  V  F  /  |  V  F | 

to  the  reflector  surface  F  (?)  *  0  is  expressed  in  the 
spherical  coordinates  (  0,  <p,  r  )  of  the  reflector  system, 
and  the  desired  surface  current  distribution  J  (r)  in  the 
Cartesian  coordinates  (  x,  y,  z  )*  then 


All  the  quatities  in  the  right-hand  side  of  the  equation  are 
either  known  or  can  be  derived  from  the  given  quantities  [2]. 
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However,  the  direct  execution  of  this  mixed  basis  cross- 
product  equation  is  usually  not  so  easy.  The  computation  can 
be  greatly  simplified,  if  it  is  converted  into  matrix 
expression: 


i  jc) 


O  n  q  a  q  c  ct*  A  c  /  ^ 

--  [CTS]  f  ns]  [STS  ]  t  r"s  3  (  Es  ) 


The  coordinate  transformation  matrices  [  T  ]  and  [STS  ]  can 
be  readily  obtained  from  the  coordinate  transformation  flow 
chart  [2] ,  therefore,  the  evaluation  of  the  induced  current 
distribution  becomes  mere  a  direct  and  straightforward  matrix 
multiplications . 


When  a  vector  potential  is  obtained,  the  magnetic  field 
can  be  determined  by  H  =  V  x  A.  If  a  vector  potential 

A  (  0,  <p ,  r  )  is  given  as  a  function  of  variables  in  the 
spherical  coordinates  but  with  three  components  in  the 
Cartesian  coordinates,  and  the  desired  field  H  (  6,  ,  r  ) 

is  to  be  determined  with  three  components  in  the  spherical 
coordinates,  then 


Converting  it  into  matrix  form,  we  have 


Or  explicitly# 


He 

0  --±-2_r  A  0..  1 

r  dr  rsin8  d£ 

i 

H,  i  « 

J  0  -  Jl 

i  «  r  be 

-Hr  J 

id  l  b  .  . 

•  rsinGd<*>  rsineSeslne  0  J 

•f 


■  COS0COS<p 

cos9sin$> 

-sin© 

' 

-sin$> 

COB<p 

0 

sinGcos^ 

sin8sin4> 

cose 

In  the  similar  manner#  the  electric  field  E  (  6#  $>,  r  ) 
is  determined  by 

es  =  -j«/i  (  1  +  -Jr  Vs  Vs.  )  Ac 

Or  in  the  matrix  form 

1  (  Es)  =  l  rs]  [STC]  (  Ac) 

f  where 

l  rs]  =  \  1  \  +  -^r  (  Vs)  {  Vs} 

Both  the  magnetic  field  and  electric  field  can  therefore  be 
obtained  simply  by  direct  matrix  multiplication  and  partial 
differentiation. 
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Quite  often  a  vector  potential  is  found  to  be  in  the  form  of 


A  (  0#  r  )  =  G  (r)  AQ(  0,  *  ) 

where 


G  (r) 


e~3kr 
4  n  r 


then  the  differential  operators  at  the  far  field  become 


jk 


[  V8]  G  (r) 


r  »  \ 


G  (r) 


0  -1 

1  0 

0  0 


0  " 
0 
0 


[  rs]  G  (r)  G  (r) 


'  1  0 

0  1 

.  0  0 


0  ■ 
0 

0  . 


CONCLUSION 

Derivation  of  compact  matrix  formulas  for  both  vector 
algebraic  operations  and  vector  differential  operations  has 
been  presented  in  detail.  The  simple  technique  is  especially 
useful  when  there  are  mixed  coordinate  basis  involved  in  the 
vector  operations. 

No  claim  is  made  that  the  matrix  formulations  presented 
are  superior  to  those  conventional  formulas  appeared  in  most 
of  vector  analysis  textbooks.  Rather  it  is  hoped  that  these 
matrix  formulas  will  become  more  popular  and  may  be  used  as 
alternatives  for  those  conventional  formulas. 


For  any  curvilinaer  coordinate  system  (  Uj ,  u2,  u3  )  with 
scale  factors  (  hj ,  h2,  hg) ,  the  following  table  may  be  useful 

TABLE  3.  MATRIX  FORMULATIONS  OF  CURVILINEAR  COORDINATES 


VECTOR  ALGEBRAIC  OPERATORS 


Scalar 


Vector 


Dot-Product  A  • 


Cross-Product  A  * 


Gradient 


Divergence  V* 


\  •  \ 


(  A  )  - 


{  A  ) 


I  A  ] 


[  '  . 


VECTOR  DIFFERENTIAL  OPERATORS 


(  V  )  - 


I  V  } 


IV]  - 


1  d 

TT*5! 

i  a 

9u2 
l  d 
du3 


[i7t2h“5u7h2b3  E^EJB^5u^b3hl  EJE^E^5u^blb2  ] 


1  *  K 

hjhjdujbl 

1  d  u 


*  v  L  A  U  t. 


1  0  K  1  0  i, 

b3^20u3°2  b2b30u23 

n  1  d  u 

°  ”  E^E^5u^h3 

1  d  „ 
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